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For the 108th Meeting of the American Astronomical 
Society, at the Maria Mitchell Observatory, the subject 
of comets was selected as the symposium topic. Reasons 
for this choice included the historic discovery of a 
comet which first brought Maria Mitchell to the at- 
tention of the scientific world, the current energetic 
research on comets which has not recently been given 
prominent presentation, and the growing prospect of 
direct space probes for which the instrumentation must 
be guided by present knowledge. 

Arrangements for the Symposium were largely made 
by Dr. Dorrit Hoffleit, Nantucket hostess, with the 
active cooperation of Dr. Herzberg, who also acted as 
moderator. To the originally scheduled set of five 
subjects, the Council at Nantucket added a brief 
presentation by Dr. Rosen of some very recent work. 
Topics and participants were: 


Presiding, Gerhard Herzberg, National Research 
Council, Ottawa, Canada 

“Structure and Development of the Gas Tails of 
Comets,” K. Wurm, Hamburg-Bergedorf Observatory 
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“Astrometric Observations and Orbits of Comets,” 
Elizabeth Roemer, U.S. Naval Observatory, Flag- 
staff Station 

“Photoelectric Photometry of Comets,” William Liller, 
Harvard College Observatory 

“Problems of the Cometary Nucleus,’ Fred L. 
Whipple, Smithsonian Astrophysical Observatory 
and Harvard College Observatory 

“The Survival of Comets and Cometary Material,” 
Ernst J. Opik, University of Maryland 

“New Investigations on the C; Spectrum at 4050 A,” 
B. Rosen, Astrophysical Institute of the University 
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Discussion following the papers was not entirely 
audible on the tape records; B. Marsden kindly helped 
to recover and edit these remarks. 
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Structure and Development of the Gas Tails of Comets 


K. Wurm 
Hamburg-Bergedorf Observatory 


I, INTRODUCTION 


T is well known that in treating atmospheric con- 
stituents of comets we have to distinguish between 
two different types of tails: those consisting of ionized 
gases and those consisting of dust. When we are able 
to look nearly perpendicular to the orbit, we see that 
the two types of tails, when existing, are well separated 
by a different curvature. The dust tail is strongly curved, 
while the gas tail is curved barely at all. A strong curva- 
ture indicates that the particles suffer a small repulsive 
force; straight tails originate from high repulsive forces, 
high with respect to the value of solar gravitation at 
the position of the comet. The main constituent of the 
gas tails is always the singly ionized carbon monoxide 
molecule. 

Twenty to thirty years ago, most of us believed that 
no serious difficulty remained in understanding the ap- 
pearance of such tails, although we realized, of course, 
that all details of their structure had not yet been ex- 
plained. It had been established that the ions of the 
tails shine by fluorescence, excited by solar light. With 
the absorption of solar radiation the particles must suffer 
repulsion by light pressure. It seemed logical, then, to 
assume that light pressure was responsible for the ac- 
celeration of the tail ions. We know now that this ex- 
planation is impossible because of the low f value in the 
tail band system of the carbon monoxide ion. 

Since we may be certain that light pressure plays no 
role, we must look for another force. In addition, the 
tails exhibit a number of peculiar phenomena which have 
yet to find an explanation. We are certainly still far 
from a theory of the tails, and no attempt will be made 
in the present paper to offer such a theory. Rather I 
have tried to prepare a critical representation of the 
main observational facts and in this way to lay the 
groundwork for the planning of new observations. 

From many recent publications on comets one might 
get the impression that a satisfactory explanation for 
the most striking observational facts is easily obtained 
by postulating high-density solar corpuscle streams and 
by considering the collisional interaction of these 


Fic. 1. Schematic dis- 
tribution of particle 
emission from a comet. 


streams with the cometary particles. Although ther 
seems to exist a rather general inclination to accep 
this view, I believe this situation has arisen only be 
cause the current conceptions concerning the observa 
tional facts are too indefinite. 

I have selected for a closer discussion three question: 
connected with observation which appear to me to b 
of particular importance with regard to a foundation 
and test of a physical theory. (% 


(1) The rate of ionization within the heads of thi 
comets. 

(2) The orbits of the tail ions with respect to thi 
nucleus and the prolonged radius vector. 

(3) The range of the repulsive forces exerted on th 
ions. 


I plan to touch only briefly on the third question 
since I have nothing to add beyond what I stated in; 
recent paper, published in the last Liége Symposium 
Proceedings (1961). In my opinion, at the present phas 
of cometary physics questions 1 and 2 deserve particula 
attention. 


II. IONIZATION WITHIN THE HEADS 


Certain explanations must be made concerning th 
monochromatic shapes of comets since these shape 
have an important bearing on question 1. Monochre 
matic refers here to the band emission of the differen 
molecules observed. Our information about monc 
chromatic shapes is gleaned by spectrograms, in pai 
ticular those with objective prisms. The situation + 
rather simple (see Fig. 1). One group of molecules fill 
only the heads of comets; they do not appear withi 
the tails proper. Most representative of this first grou 
are the diatomic cyanogen molecule and the diatomi 
carbon molecule, both of which generally appear wit 
a high emission intensity. A second group comprise 
exclusively ions; their zone of emission is indicated 1 
the drawing of Fig. 1 by the shadowed area. There ° 
some overlapping of the regions of group'I and group 
Some uncertainty exists with regard to classifying th 
two ions near the end of group I: CH*t and OH 
Because of their faintness they are observable with the 
bands only in the vicinity of the nucleus, but very prol 
ably they should be moved into group II. We know thé 
all these particles, ion as well as neutral, start the 
expansion into space from the vicinity of the nucleu 
But how close to the nucleus the ions first appear is n¢ 
precisely known. From pictures of the recent Com 
Burnham 1959 k, which I note later in another conne: 
tion, I found that the COt ions were already certainl 
present at a distance of 5000 to 6000 km from th 
nucleus. 
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(The expansion of the neutral molecules proceeds from 
je nucleus region more or less symmetrically in all 
‘ections, filling their space continuously. If we identify 
je head of a comet, as I will do later, with the extension 
4d space of the neutral particles, we may say that the 
jad is approximately a sphere with the nucleus in the 
Jater. An informative paper on the expansion veloci- 
i's of the head particles was published more than 30 
Jars ago by Bobrovnikoff (1931) in his important study 

i alley’ s comet. For a time after perihelion this comet 

owed brief discontinuous outbreaks of matter from the 

icleus which resulted in expanding halo formations. 
y comparing plates from different observatories over 
e world, Bobrovnikoff was able to study the expansion 

id lifetimes of these halos. He found the expansion 

locities rather small: They averaged below 1 km 

'r second. The average lifetime was about two days, 

ter which they lost definite boundaries because of their 

intness, having reached radii of 130000 to 150 000 

n. The comet’s position at this time was about 1 a.u. 

om the sun. 

The radii of the cyanogen and carbon heads of comets 

‘e generally given from photographs as 150000 to 

0.000 km, We are well aware that no sharply defined 

idius exists. What we should investigate are the in- 

sity gradients within the monochromatic images for 
larger number of comets, and derive from them the 

w of density distribution for a regular expulsion of 

latter from the nucleus. This density distribution 

aould depend on two factors. If p designates the dis- 
ance from the nucleus, the first factor is p-? because 

f the expansion of the gases from the nucleus in all 
‘irections; the second factor e~°/*° enters since there 
aust be a decay of the particles. po is defined by po=27, 

vis the expansion velocity and 7 the average lifetime 
_ f the molecule. 

A careful study of the density distribution should 
llow us to find the average lifetimes of the molecules, 

_ iving the rates for ionization plus dissociation. Studies 
-f monochromatic intensity gradients have not yet 

‘een numerous enough to provide accurate values of 
_ he average lifetimes 7. Isophotic contours or intensity 
_ radients for CN and C2 have been published by Yoss 
1953), Miller and Liller (1956), Miller (1957), and 
/orontsov-Velyaminov (1961). From the material 
_vhich is available at present, one can conclude at least 
hat up to a distance p=200 000 km the influence of 
he decay factor is not yet very strong. It lowers the 
lensity perhaps by a factor of 2 or 3, but certainly not 
nore. 

Hence, we may conclude that the constant po > 200 000 
cm. If we adopt an expansion velocity of 1 km/sec, a 
ittle larger than the average measured by. Bobrovni- 
soff, we find for the lifetime 7> 10° sec. This figure is 
yppropriate for the molecules CN and C, and for helio- 
rentric distances equal to about 1 a.u. From the in- 
tensity gradients along the band lines of the molecules 
9H and NH it follows that the same lifetime is also 
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valid for these molecules. This figure, 10°, is not far 
from the one which is given by Nicolet (1961) as the 
average lifetime of the neutral oxygen molecules in the 
uppermost levels of the earth’s atmosphere, Nicolet’s 
value being 10° sec. Thus the behavior of the observed 
neutral particles in the atmospheres of comets points 
to a rather low ionization rate for particles with ioniza- 
tion energies of, say, 14 ev. 

A striking contradiction will appear when we analyze 
the behavior of the ions and derive from them, insofar as 
possible, the lifetimes of their parent molecules which 
are, of course, identical with the neutral molecules. To 
be sure, the comparison of the behavior of the ions and 
that of the neutral molecules is somewhat complicated 
because we are not in a position to observe the emission 
of the same molecule in its neutral and ionized state. 
Although the band emission of the neutral particles 
mentioned earlier is strong, the particles do not appear 
with a sufficiently strong emission in their ionized 
state, and those with a strong emission in their 
ionized state do not show a sufficiently strong emission 
in their neutral state. However, no particular difficulty 
exists in analyzing the behavior of the ions with respect 
to the ionization of their parent molecules. We do limit 
our considerations to the CO molecule, since the general 
situation with regard to it is clearest. 

Some years ago I gave (Wurm 1947) a figure for 
the ionization rate of the parent molecule of the COT. 
The average lifetime of the neutral CO was estimated 
at rco= 10? sec. This estimate was derived from two 
different observations. First, since we observe COt 
emission already close to the nucleus, ionization must 
take place before the neutral molecules are able to travel 
far from the nucleus. According to this, the lifetime can 
be at most a few hours. The second and even more con- 
clusive observation concerns the fluctuations in the 
stream of CO+ ions moving from the region near the 
nucleus into the tail. Suppose we have a cloud of CO, 
in which by any process ionization is occuring with a 
probability \=1/7co0; suppose further we are observing 
the created ions streaming away from the cloud. A 
stationary state can be reached if first there exists a 
source which keeps the total mass of the CO cloud 
constant, and second if this constancy has persisted for 
a time (>>rco. Fluctuations in the streams may appear 
either by fluctuations in the available mass of the CO 
or by fluctuations in X. It can easily be shown from the 
differential equations of the process that, if we observe 
a fluctuation in the ion stream within a time interval 
At, then rco< At. It has long been known that fluctua- 
tions in the CO* stream of the tails become perceptible 
in less than‘an hour. This means that the average life- 
time of the parent molecule of the CO* must be smaller 
than 10*-* sec. 

When we compare this figure with that which we de- 
rived for the typical head particles, there seems to exist 
an obvious contradiction. The head molecules CV, C2, 
OH, NH have approximately the same ionization energy 
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(b) 


Fic. 2. (a) Comet Burnham 1959, photographed April 25-26 
with 48-inch Asiago reflector. Exposure 15 minutes. (b) Same 
comet, 70 minutes later. 


as CO, but we find a difference for the ionization rates 
which amounts to two powers of 10. Since it had long 
been the prevalent belief that the high rate factor found 
from the carbon monoxide ion was the true one, not 
much weight was given to the existence of neutral 
particles; they were regarded simply as a residue which 


had not yet suffered ionization. According to what 
have explained, this conception is certainly not tenabk 
For a time I thought that a misinterpretation of observe 
tion was possibly connected with the low rate factor 
however, I no longer hold this opimion. The situatio 
now appears to me perfectly clear with the neutr 
particles, far more so than with the ions, for which reg 
son I believe we shall have to look to them for the fe 
moval of the discrepancy. Since either derivation seem 
reliable enough (excluding a drastic change of th 
figures), I suspect we must face the conclusion that th 
discrepancy arises from two different processes of ion 
zation, one for the comet at large originating from a 
external radiation (light and corpuscles), which thu 
far has been the only one considered, second an intern 
process of ionization, limited to the vicinity of th 
nucleus, with a high rate but with an incomplet 
efficiency. I can only, considering the whole situatiot 
make such a suggestion; I am not in a position to di 
velop a picture of this internal source of ionizatior 


III. ORBITS OF THE TAIL IONS 


Coming now to the tails proper, I should summariz 
briefly the most striking differences in the behavic 
of the ions and the neutral particles. In the first plac 
are the much higher repulsive forces exerted on th 
ions. The real range of the repulsive forces is only aj 
proximately known; certaily values 50 to 100 time 
the gravitational attraction of the sun are reachet 
Much higher values mentioned often in the literatur 
such as a thousand to ten thousand times the solar a 
traction, remain very questionable. At least statist 
cally they cannot be important (Wurm 1961). Secondh 
regarding the peculiar behavior of the ions, I woul 
mention their strong tendency to form rays an 
streamers. Normally they do not fill their space cor 
tinuously; they arrange themselves in rays with dian 
eters, as far as can be measured, equal to or below 20€ 
km. When beginning to form, the rays have an ang 
exceeding 45 deg to the tail axis and are convex again: 
the sun. While growing in length [see Fig. 3(c) ] th 
angle decreases; the ray lies finally along the tail ax 
and dies out. An essential difference between the ior 
and the neutral particles concerns the initial velocitie 
which in the case of the ions are much higher than 1 
the case of the neutral molecules. But we must make 
restriction here with regard to the phase in the formatio 
of the rays. The initial velocities are high at the b 
ginning of the formation but small at the end of it. 

Figure 2 shows particularly clearly the turning of ta 
rays to the tail axis. The photographs are of the recer 
Comet Burnham 1959 obtained on the night Apr 
25-26 by Dr. Maffei and myself with the 48-inch Asiag 
reflector. The comet was difficult to guide since it ha 
no clear nucleus. The two outermost rays of the bund. 
formed, on the first exposure, an angle of 19° 40’; 7 
minutes later they had decreased to 11° 30’. Furthe: 


re, on the last picture the two rays mentioned have 
jome somewhat less diffuse than on the first one. It 
ms rather evident that the cause is to be found in a 
brease of the angular velocity of the turning. The 
‘ee or four innermost rays are all less diffuse and 
irper than the outer ones. They are either already 
irly at rest or are moving in a plane which contains 
» line of sight. 
We can conclude directly from certain tail configura- 
‘ns that the initial velocities of the ions within the 
‘| rays are large at the beginning of the ray formation 
F d small] at the end of it. The configurations in question 
‘nsist of very long rays which have accumulated along 
| » tail axis (or radius vector, which is approximately 
-}‘ntical with the tail axis) having lengths of 10 to 20 
t llion km, and very short rays with lengths of the order 
é the bead radius (100000 to 200000 km). A good 
jae of this kind of tail was found in Comet Hassel 
139 at the end of April 1939. For heliocentric distances 
Jtween 0.5 and 1.0 a.u., and assuming reasonable 
‘oulsive forces [(1—y)=50 to 100] we can easily 
«mpute that the ions at the ends of these long rays 
‘It the nucleus 50 to 100 hr earlier than when these 
otographs were taken. From the widths of such long 
yy complexes it follows that the ions have distances 
jm the tail axis less than the radius of the head (say 
+0000 km). Hence, their velocity normal to the radius 
» ctor must have been on the average less than 1 km/sec. 
1 the other hand, the ions at the tips of the short rays 
\ve already reached this distance (about 200000 
~n) within time intervals of the order of one hour! This 
ugh picture suggests initial velocities normal to the 
il axis of 50 to 100 km/sec at the beginning of the ray 
‘rmation and of 1 km/sec or less at the end of the 
rmation. Whether the initial velocity has only a com- 
ment normal to the tail axis is not yet clear. From 
_ is explanation it follows conclusively that the turning 
_ the tail rays to the tail axis is kinematically connected 
. ith and dependent on a decrease in time of the velocity 
_mponent normal to the tail axis. Within each ray 
lis component is, on the average, smaller the later an 
} leaves the region of the nucleus. It remains uncertain 
/hether this velocity can always be identified, over the 
“ital orbit of an ion, with its initial velocity len leav- 
g the nucleus region. 
Stumpff (1959), in treating the kinematics and dy- 
amics of tail ions, concluded that the turning of the 
lil rays can be understood quite simply by the action 
' repulsive forces alone. His idea becomes clear with 
jie aid of the diagram of Fig. 3(a). The straight line 
represents a ray at a certain time /, and the straight 
ne B the same ray at a later time /). The author con- 
ders essentially only the problem of transferring A 
ito B.by remaining in accord with the time scale of 
ve whole turning of the ray. He neglects as outside 
is problem the early phases of the rays when their 
jimensions are comparable with the diameter of the 
ead, and he also ignores the final fate of the rays. For 
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(1-u)=f(y) 
y 


(a) (b) (c) 


Fic. 3. Turning of the tail rays of a comet. 


his restricted problem he finds a solution if he postulates 
a repulsive force which is a function of the distance y 
from the tail axis. The numerical values are about 
(1—p)=10 along the tail axis and (1—) between 10° 
and 104 at distances y>200 000 km. I believe it is un- 
necessary to explain that Stumpff’s conception is not 
tenable. 

In a short note several years ago (Wurm 1953) I 
pointed out that the turning of the rays can be under- 
stood kinematically, at least to a large extent, by a 
constant repulsive force including a steady decrease of 
the velocity of expulsion. Figure 3(b) shows. qualita- 
tively for these conditions the movement of the indi- 
vidual particle: Point a goes over into point a’, b into 
b’, and c into c’. All points move away from the tail 
axis, but each tail ray as a whole turns toward the radius 
vector. This suggestion was only intended to emphasize 
that the turning of the rays depends on a decrease of 
the velocity component vy. 

At first sight, it might seem easy to investigate the 
orbits of the tail particles—one could simply observe 
shifts in time of well-defined condensations within the 
tail rays. The obstacle is that no well-suited material 
exists for this purpose. The time interval between two 
photographs which one can compare is for the most part 
too short to reveal a clear shift of a condensation. I 
have studied carefully any suitable pairs of photographs 
which I found, but very few of them gave any informa- 
tion. While the available material is insufficient to give 
a clear picture of the movement of individual particles, 
there is some indication that the velocity component 
normal to the tail axis is not identical over the whole 
orbit of the ion with its initial velocity v,. It is rather 
probable that a new force appears which accelerates 
the ion in the direction toward the radius vector, di- 
minishing in this way the initial velocity v,. At first 
sight, the appearance of such a force would seem to cast 
doubt on our reasoning concerning the high initial 
velocities of the ions at the beginning of the ray forma- 
tion and their low velocities at the end of the formation. 
However, it can be shown that the ions within the long 
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tail rays have low velocities from the beginning. The 
situation is not so simple that I can explain the reasons 
for this opinion in a few phrases. I intend soon to publish 
a note on this point. 


IV. CRITICISM OF THE CORPUSCULAR CAUSE 
OF THE TAILS 


The subjects touched in this paper have not, of course, 
exhausted our total knowledge about the tails, but the 
points mentioned are certainly of prime importance with 
regard to the foundation of a physical theory. I do not 
intend to deal more closely with any other observations, 
but now wish merely to summarize the arguments which, 
in my opinion, speak against the corpuscular-stream 
theory of the tails as suggested by Biermann (1953) 
and advocated by him and his associates at the Max 
Planck Institute. 

It is claimed that the theory accounts as well for 
the ionization within the heads as for the acceleration 
of the ions. The ionization is ascribed to a charge-trans- 
fer process between the solar protons and the cometary 
neutral particles. In this connection we need not discuss 
whether the proton densities are high enough to explain 
the high ionization rate of the carbon monoxide or 
carbon dioxide. The real problem is to reconcile the high 
rate near the nucleus with the low at larger distances, 
and it is this problem for which the theory can present 
no answer. I am inclined, as indicated already, to con- 
sider this contradiction between the ionization rates as a 
hint not to connect the observed ionization near the 
nucleus with the direct influence of an outer force or 
radiation, but to search for its energy source within 
the comet. Support is lent this belief if we consider 
the two other phenomena connected with the appearance 
of the ions: the formation of the ray structure and the 
high and regularly decreasing initial velocities. We 
know nothing about their causes but it seems to me 
impossible to imagine that the regular development of 
the tail rays, each ray developing largely independently 
of the others, is a phenomenon which can be controlled 
directly by an outer influence. 

The theory runs into difficulties also with respect to 
the repulsive forces, this despite the fact that we need 
not consider at all the extreme values of 1000 to 10 000 
times the solar gravitation which were envisaged by 
Stumpff to explain the turning of the tail rays. Ac- 
cording to Biermann, repulsive forces of the order 100 
times the solar gravitation at 1 a.u. from the sun can 
be explained by the collision theory with proton densi- 
ties of 10° particles per cm’. But we must remember 
that this result depends on an assumed density of the 
cometary particles of the order one to ten. The density 
in the accelerating stream must be by about three 
powers of 10 higher than that of the accelerated stream. 
According to newer estimates the ion density within 
the tails of brighter comets amounts to 10? to 10%, and 
we would correspondingly need proton-stream densities 
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of 105 and more in order to account for the observ) 


tt 


accelerations. Very probably, solar corpuscle streai) 
of such densities can, if they exist at all, be only ol 
localized and Se orlic nature. But because of th} 
distinctly permanent nature it is not possible to conn(| 
the cometary phenomena which I have described ¥ Wi] 
such occasional streams. ) 

One can also find strong arguments against the hig 
density corpuscle theory in a consideration of the exci} 
tion conditions of the emitting molecules. I will menti} 
only one argument, connected with the oscillation a! 
rotation distribution with the C2. Swan bands. The | 
molecule cannot get rid of its rotational and oscillation 
energies by way of radiation since it has no permane 
dipole movement. Collisions with fast protons show 
of course, increase the internal energies of the molecul 
The C2 molecules at the border of the head have exist 
about two days within the atmosphere. There is not t} 
slightest indication that the molecules near the bord 
have gained rotational or oscillational energy compar 
with the molecules near the center. If one assumes ¢ 
effective cross section of 510-16 cm? for the ener; 
transfer, one estimates that a corpuscle stream with 
density howe 100 protons per cm’ cannot exist. 

I believe that the arguments mentioned in this pap 
together have sufficient weight to show that we mu 
give up the collisional theory of cometary tails, at 
search for another one. 


DISCUSSION 


Oprx: Concerning the need for high particle densi 
in the corpuscular streams, I think the observations” 
the spectrum of the zodiacal light by Blackwell she 
that the electron density is certainly less than 100/cn 
To be sure, if the stream velocity is smaller, the inte 
action rises since the cross section is inversely pro 
tional to the 4th power of the velocity, but then 
momentum is insufficient. Perhaps a greater difficu 
is the motion of the particles in the cometary ta 
Repulsive forces of 100 or 1000 times solar were report 
first by Bobrovnikoff with Comet Halley. From A 
photographs made in late May and early June 1910 
assumed that every day there were separate lems 
But if all the observations are joined together they sh¢ 
a beautiful single trajectory, with average radial 1 
pulsive force of from 60 to 70—quite normal—still t 
motion is not that of a central force, the transvet 
component coming to a full stop. And this is similar 
what Dr. Wurm has told us. 

WuippLe: Is there any well-substantiated eviden 
of correlation between tail activity and solar activit 
At Liége there have been several papers on this. 

Wurm: There have been many papers, but I don 
believe the evidence is good. 

: Are effects of rotation small compared to t 
forces you have described? 

Worm: We don’t know what the forces actually ai 


. | otation of the tail rays is certainly not observed. 
“jiprk: In the case of Halley’s comet the force was 
4\ctly in the direction of the radius vector. 
\ferzBerc: You said that there was an apparent 
“tradiction between the two lifetimes that have been 
aimed? If I understand correctly one is for CO, and 
lother is for CV, Cy and some others. Is it not con- 
jrable that these should be different? 

Vurm: No, I don’t think so. Observed also are CHt 
#| CH neutral. 

TerzBerc: Yes. Do I understand correctly that you 
7 ibute the large difference by a factor of 100 to some 
7-hanism inside the comet? 

Vurm: There seem to be two mechanisms for ioniza- 
/). The internal one somehow causes the lifetime of 
y molecules in the nucleus to be very short, yet 
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this mechanism is not fully efficient and is restricted in 
space. 

HERzBERG: Perhaps some kind of shock wave that is 
produced by the breaking up of the core? 

Worm: This would be difficult since the free path for 
particles is of the order of the atmosphere. 
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Astrometric Observations and Orbits of Comets 


ELIZABETH ROEMER 
U.S. Naval Observatory, Flagstaff Station 


NY study of comets begins with observations of 
position. Crude observations may suffice for 
finding purposes for spectroscopic or photometric ob- 
servations of brighter objects. But a large amount of 
information can be gained, for faint objects as well as 
for bright ones, by studying the motions, or orbits. 
Precise study of the motions of comets has two major 
aspects: observation and computation. Observations of 
position are desired, of considerable precision, and ex- 
tending over as long an interval of time as possible. 
To achieve either of these ends, the application of fairly 
large instruments is needed. For astrometric purposes 
comets are not to be regarded as extended objects; the 


advantage in scale and limiting magnitude of large’ 


telescopes is valuable. 


OBSERVING PROGRAM AT FLAGSTAFF 


At the Flagstaff Station.of the U. S. Naval Observa- 
tory an effort is made to obtain one pair of plates during 
each dark of the moon of every comet within reach of 
the 40-inch Ritchey-Chrétien reflector. Additional ob- 
servations may be made if they are likely to be useful, 
as, for example, in the case of newly discovered objects. 
This program is a modification and expansion of the 
type of program carried on for some time by H. M. 
Jeffers of the Lick Observatory, and especially by G. 
Van Biesbroeck at the Yerkes and McDonald Observa- 
tories. Moving objects are recorded regularly to magni- 
tude (pg) 20.5. Since the dividing line between short- 
period comets and minor planets is by no means clear, 
certain minor planets are observed whenever the oppor- 
tunity arises. Some of these minor planets may indeed 
be old comet nuclei. Included on the program are all 
the minor planets having a perihelion distance less than 
1 a.u. and many of those that have a perihelion distance 
only slightly greater than 1 a.u., and for which ephemer- 
ides are available. 

To demonstrate the extent of the program, it may be 
noted that 19 different comets and nine minor planets 
were successfully observed during the 12 months ending 
in June 1961. Eight returning periodic comets were re- 
observed first at Flagstaff, including four within 16 
days in August 1960. 

In making the observations, the calculated motion 
of the object is offset:as a general practice, the plate 
being displaced relative to the crosshairs of the guiding 
eyepiece. Guiding is done by means of a star located 
near the edge of the field being photographed. The tele- 
scope operates at //6.8, and exposures up to two hours 
in duration are useful. Appropriate exposure time is 
expecially important in astrometric investigations. 
Since estimates of the brightness of comets vary greatly, 


depending on the instrument and methods by which {j 
observations are made, we have found it necessary: 
calculate our own version of predicted magnitudes, (|) 
pending upon earlier observations of our own, or th 
of other observers using similar equipment and metho 
The results of Cunningham (1951, 1952), obseryi} 
with the Mount Wilson 60-inch and 100-inch telescop 
1950 to 1953, have been especially valuable. il 
Bright comets, as well as faint ones, are observ} 
regularly. Especially in bright objects is there pos 
bility of a systematic displacement of the center | 
light from the center of gravity. Such a displaceme| 
leads to systematic errors in the orbit, especially 
the case of such objects as periodic Comet Encke, whil| 
is usually observable, during any specific apparitid 
only before, or only after, perihelion passage. To fl 
best of my knowledge, the U. S. Naval Observator 
Flagstaff Station, is the only observatory in the world | 
which astrometric observations of comets fainter tha 
mag. 17 are being made regularly at the present tin 


if 


i 
MEASUREMENT AND REDUCTION OF PLATES 


Since the field of the 40-inch reflector, or a, 
similar instrument, is fairly small, reference stars ¢ 
most invariably have to be taken from the Astrographi 
Catalogue. With the epoch of most plates of the Cats 
logue 40 to 60 years ago, proper motions have had 4 
appreciable effect upon the star positions. The call 
tion of the plate constants of the Catalogue often it 
volves meridian positions of reference stars from 2 
even earlier epoch, used without application of pro 
motions. The great improvement in the accuracy 
positions reduced from the portion of the Paris 
for which new plate constants have been calculated 
Heckmann, Dieckvoss, and Kox (1954) is particulai 
noteworthy. Even with reference stars taken from 
Astrographic Catalogue, there is an appreciable differe 
between the magnitudes of the moving objects bel 
investigated and those of the reference stars. Hen 
there is likelihood of systematic errors. 

Reduction of the measures may be carried out by 
of several methods: by dependences, by a modi 
plate constant method, or by accurate orientation 
the plate in measurement and reduction by means 
an assumed value of the plate scale. The dependen 
method, with at least six reference stars and a leas 
squares solution, is used for the reduction of most | 
the results from Flagstaff plates, the computations b 
ing executed with the IBM 650 of the U. S. Nav, 
Observatory, Washington, D. C. In emergencies, whe 
results are needed immediately, the plates are generall 
oriented accurately and the reductions carried out b 
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d at Flagstaff. In our experience these methods pro- 
: better checks in dealing with the Astrographic 
) logue. 


CALCULATION OF ORBITS 


inly very rarely has every precaution been taken 
Walculation of orbits to ensure that the calculations 
resent the actual motion of the object to the greatest 
wsible physical accuracy. Very few of the persons 
4g the computing have had observational experience, 
cially with large instruments. They are, therefore, 
7 always well equipped to evaluate the quality of the 
ervations on which the computations are based, or to 
i with the observations in the most effective manner. 
‘ey i is the attempt made to reduce the observations 
a consistent catalogue system, even if this can be 
ie. 
‘n general, a few positions carefully measured on 
perly exposed plates taken with an instrument of 
‘cient focal length and well spaced in time are worth 
re in practice in defining the orbit than great num- 
is of positions obtained only when the comet is 
ght. Computers would often have an easier time, 
w hout loss of precision and sometimes with gain, if an 
a reciable fraction of the total number of observations 
wte simply disregarded. This can be seen in those in- 
snees in which residuals of individual observations 
fm an improved orbit have been published in detail. 
Especially for short-period comets, the effort in 
i dictions has generally been to obtain sufficient ac- 
«racy for recovery. The efforts of the volunteer workers 
«the computing section of the British Astronomical 
‘sociation have been especially noteworthy and valu- 
4\e in making observations possible at repeated returns 
ta numerous members of the short-period Jupiter 
nily of comets. Necessarily, approximations have 
fen introduced in computing the perturbations, and 
“ne significant perturbations have been neglected, 
i to keep the volume of work within bounds. 
arting orbits have not always been the best available, 
‘d almost never the best possible. And, inevitably, 
stakes in the calculations have occurred from time to 
he. Because of misunderstanding, the idea has grown 
lat there are appreciable nongravitational forces opera- 
‘ve, and that precise predictions of future positions of 
mets are impossible. 
It should be noted that observations at one apparition 
‘not, however, determine the value of the mean mo- 
m sufficiently well to expect an exact prediction at 
je next succeeding perihelion passage. Furthermore, 
e probable error in the mean motion deduced from a 
ist-squares correction to the elements may be largely 
‘titious. This fact is well demonstrated by Merton’s 
1927) work on periodic Comet Schaumasse. A least- 
juares solution by G. Fayet and A. Schaumasse from 
}e observations in 1919 gave a value of the mean daily 
jotion, n=439"476+0"7. By linking the apparitions 
1911 and 1919, Merton found the value of the mean 
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motion to be ~=446"811. This value of m made the 
predicted time of perihelion passage in 1927 some 49 
days earlier than the 1919 observations alone indicated. 
In spite of decidedly unfavorable circumstances the 
comet was recovered very close to Merton’s prediction. 
Another example is the recovery of Comet Gale by 
Cunningham in 1938 some 5” 23™ west of what appeared 
to be the most likely place predicted on the basis of 
observations in 1927 (Davidson 1939). 

Few orbit methods give a good indication of the real 
physical determinacy. One notable exception is Cun- 
ningham’s method, in which a change of the coordinate 
system yields explicitly the deviation from great circle 
motion and the length of arc represented by the funda- 
mental observations. A similar insight into the physical 
accuracy may be obtained by computing variational 
orbits, on each of which the observations may be 
represented. This becomes practicable now with the 
application of large digital computers. 

The correction required to bring the predicted time 
of perihelion passage into accord with observation for 


TaBLeE I. 
Obs. AT 

Comet returns (days) Computer 
Reinmuth 2 3 0 E. Rabe 
Finlay 7 —1.2 M. P. Candy 
Comas Sola 5 (small) H.Q. Rasmusen, J. M. Vinter 

Hansen 

Harrington 2 —0.6 C. Dinwoodie, B. G. Marsden 
Brooks 2 10 (small) A.D. Dubiago 
Encke 46 0 S. G. Makover 
Schwassmann- 

Wachmann 2 6 —0.3 H. Q. Rasmusen 
Borrelly 7 +0.2 M. G. Sumner, M. P. Candy 
Wirtanen 3 0 P. Herget 
Forbes 4 +0.5 3B. G. Marsden 
Tempel 2 13 +0.4  B.G. Marsden 


periodic comets reobserved in 1960-61 is shown in 
column 3 of Table I. The total number of observed 
perihelion passages of each comet appears in column 2. 
It should be noted that if two perihelion passages have 
been observed and an effort made to obtain an accurate 
prediction, the divergence of calculations from observa- 
tions is very small. The predictions for comets Reinmuth 
2, Encke, and Wirtanen are examples. 


SECULAR ACCELERATIONS 


It is evident from the activity photographed in the 
head and tail of comets that some nongravitational 
forces do indeed exist. But the existence of measurable 
nongravitational accelerations of the center of mass 
should be regarded as an open question. The amount 
of any such nongravitational accelerations is vital to 
theories regarding the structure of comets and their 
rates of disruption. See, for example, the papers by 
Whipple (1950), and van Woerkom (1948). 
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It is a revealing experience to look critically into the 
published computational work, on the basis of which 
the existence of nongravitational secular accelerations 
has been asserted. Divergences from observations are 
quite small, and the quality of the physical parameters 
(such as planetary masses) and accuracy of computation 
that have entered into even the most precise orbits are 
disappointing in comparison with the size of the 
residuals. 

As an illustrative example, let us take periodic Comet 
Wolf 1, 1884 III. I choose this object only because the 
published details of the calculations are quite complete 
and readily available, and also because I am responsible 
for all of the accurate positions obtained of this comet 
at its most recent return in 1958. Kamienski, who has 
devoted more than 50 years to study of the motion of 
this object lists (Kamienski 1953, 1959) 72 papers he 
has published to report his findings. The quality of the 
computations, carried out by hand methods, is quite 
remarkable, and this object is nearly unique in that 
perturbations at the most critical times have been 
checked by independent calculations by completely 
different methods. The agreement of these duplicate 
calculations generally has been quite good, i.e., prac- 
tically to the limit of the accuracy carried. Nevertheless 
the small differences that do exist probably are signifi- 
cant in view of the size of the observational residuals, 
on the basis of which a secular acceleration is deduced. 
The system of planetary masses used involves differ- 
ences in the mass of Mercury of some 50%, and of Jupi- 
ter of 0.1% from the best modern values. These differ- 
ences cabice be considered negligible when the attempt 
is made to represent observations to 1” accuracy. 

The observed motion of P/Wolf 1 may be considered 
in two parts, separated by a close approach to Jupiter 
(Amin=0.12 a.u.) in 1922. Prior to the close approach 
the comet moved in an orbit of perihelion distance 
1.61 a.u., and afterward in one of perihelion distance 
2.45 a.u. Kamienski (1933) found it necessary to intro- 
duce an acceleration in the mean motion to represent 
observations in the smaller orbit: 


An=—0"000 00042 (t—%), 


where (i—¢o) is expressed in days from t=1884 
September 24.0 BMT. The apparitions of the comet in 
1925, 1933/34, and 1942 were linked with one set of 
elements, which, with the inclusion of perturbations 
due to six planets, represents also observations of the 
comet in 1950-51 and 1958 with maximum residuals of 
—0%3. Kamienski suggests that the residuals may be 
caused partly by uncertainty in the amount of per- 
turbations by Jupiter in the mean anomaly, one part 
in 8000 being sufficient to reduce the residual to zero, 
and partly by a continuing secular deceleration. In 
the larger orbit Kamienski finds (1961) a value for the 
deceleration term of 


An=—0"000 00012 (t—t») 
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-computations, and physical constants capable of i 


(with time counted in the same manner as tome | 
the deviation from theory in 1958 is attributed e 
to nongravitational forces. This acceleration, one-il h 
as large as that found in the smaller orbit, is basedy). 
an observational residual of (a—a,) coss=—O\y. 
=—4"65 on June 22, 1958. oa 
As the person responsible for the observations, I a | 
not feel that an acceleration deduced from such resi 
is convincing. Kamienski has implied the camel | 
clusion by suggesting that judgment should be withhi 
until observations can be obtained in 1967. Tat | 
portant point to be made here is that gravitai 
theory alone,-with borderline numerical accuracy 


provement, is adequate to represent observations © 
the entire time interval from 1925 to 1958 with 
residual larger than 4’’—a residual quite within t 
limits of possible observational error. There is, therefo} 
no evidence for appreciable nongravitational fori) 
operating on this comet during this interval. Da: | 

The observational residuals are somewhat a i 
the motion in the smaller orbit, prior to 1922, dt 
residuals are larger for Comet Encke, whic ha as 
perihelion distance of 0.3 a.u. But both computati 0 
and observations are more sensitive to effects that | i 
to the appearance of a run of residuals. There is m } 
likely to be a systematic displacement of the center 
light from the center of gravity in the brighter, asj 
metrical image of a comet near the sun. The time inte 
val of steps in calculating perturbations must 
shortened, leading to a rapid increase in the numbel 
steps, and hence! to an even more rapid buildup 
roundoff errors in numerical integration. The appro: 
mation of a barycenter is less valid. And the pertur } 
tions by planets whose masses are subject to consi¢ + 
able uncertainty (Jupiter, and, formerly, Mercut 
and Venus) may be appreciable. 

With the emphasis that has been placed recently © 
comets as probes in studying the interplaneti 
medium, and the inferences that have been drawn f 
garding comet models, it may not be out of order 
suggest that a great deal of work needs to be done bef 
much confidence can be placed on conclusions in 
realms. A definitive study of the motion of seve 
well-observed objects would be extremely useful. Obje 
suspected of a secular acceleration might supply es| 
cially interesting information. Such a study must st 
with a rediscussion of original observations, with effo: 
being made to emphasize observations carefully plann | 
and executed with moderate- to large-sized instrumen 
An attempt should be made to reduce positions to 4 
consistent catalogue system. Calculations should — 
carried out to full physical accuracy, and with sufficient 
guard figures to ensure that such factors as round¢ 
error in long-continued numerical integrations do 
destroy significant results. It would be extremely in 
esting to compute perturbations around a family 
variational orbits that do not do an injustice to observ 


), rather than for the one orbit that gives the best 
hematical fit to observations. The differences in 
orbits at aphelion, where the perturbations by 
ter are largest, may be important. A true estimate 
“els to be made of physical uncertainties in predic- 
ics, to ensure that an adequate search is made. Peri- 
comets undoubtedly have been lost in the past 
“41 are being lost at present) because of inadequate 
“@ches. It should be possible to verify results by 
‘@putations by different methods, so that mistakes 
omputation may be discovered and corrected, as 
as to compare the relative merits of different 
iniques. A great deal of information undoubtedly 
» be extracted from existing observations with vigor- 
and imaginative use of large-scale computing facili- 
now available. And in the near future precise orbit 
§i undoubtedly will be useful in planning the direct 
| of the properties of comets by means of 


ce probes. 
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lite 
'Wurppte: This is a very fine work. I should like to 
/n Miss Roemer in her plea for recalculation of some 
-yits. Encke’s comet especially needs it, for in its early 
“ys especially the perturbations could not have been 
iplied very accurately. Cunningham has data, but has 
»t reduced it as yet. If I had suspected that ere may 
“no secular changes in the period, I don’t know whether 
would ever have published my icy-conglomerate 
ddel! (Laughter.) 
ErcuHorn: If you say the accuracy limit is imposed 
‘the Astrographic Catalogue (can’t do better than 
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1” or 2’), do you attribute this to intrinsic inaccuracy 
in the positions in the AG, or to proper motions, or to 
the uncertainty in tying the astrographic measures 
into a certain fundamental system? 

RoEmeER: It is really a combination of the quality of 
the reference star positions (primarily observed for the 
AGK1i and used without proper motions) and the long 
time interval involved since. 

E1cuHorn: Are any systematic corrections applied? 

RoEMER: This is practically impossible to do, and it 
would probably cause damage. Let the person using the 
positions apply the corrections he considers desirable. 
Of course there is the possibility of using field plates, 
reduced using Yale stars or AGK2 stars, but this more 
than doubles the work. 

——: How large are the typical cometary nuclei that 
you measure? 

RoEMER: One practically never observes the nucleus 
directly. I would give very great weight to the visual 
observations of van den Bos, Baldet, and V. M. Slipher 
of the nuclei of comets close to the earth, the indicated 
size being from 4 to a few km for a typical object. 

HERZBERG: You mentioned an object with two nuclei. 
The motion with respect to each other has been 
measured ? 

Roemer: I think one can get the velocity of separa- 
tion, which gives some index of the energy of processes 
going on in the nuclei. This comet was at such a great 
distance from the sun that the disruption was not a tidal 
effect, but something that happened internally in the 
nucleus. 

HeERzBERG: Would you care to comment on the re- 
semblance between comets and asteroids? 

RoEMER: The Jupiter family comets cannot be dis- 
tinguished from asteroids on the basis of orbits alone. 
When I observed Comet Arend-Rigaux on its last ap- 
parition I found that it was completely stellar on all 
plates. The orbit is similar to that of a minor planet, 
but the object was designated as a comet because it on 
occasion showed some diffuseness. When Baade dis- 
covered Hidalgo he was undecided whether to call it a 
minor planet or a comet, but he decided on the former 
simply because more people were observing minor 
planets at the time and it would be better taken care of! 

(Laughter.) 
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cS ee recent appearances of several bright comets 
have resulted, happily, in an increased interest of 
astronomers with access to photoelectric equipment in 
comet research. The new observational data have, in 
turn, provided us with a better understanding of the 
physical conditions in comets. Most of the observations 
have been made in broad spectral regions defined, for 
example, by the standard UBV system. A number of 
astronomers have since acquired narrow-band filters in 
the event another bright comet puts in an appearance 
in the near future. The use of interference filters or 
spectrometers is to be very much preferred over broad- 
band photometry, of course, because most comets 
exhibit spectra composed of both molecular emissions 
and continuous radiation. When using broad-band 
photometry, one cannot always be sure whether changes 
in intensity (as a function of time, of plane of polariza- 
tion, etc.) are to be attributed to the emission bands or 
to the continuum. There are times, however, when 
virtually all the radiation is in either the bands or in 
the continuum. Such a circumstance may occur in a 
comet tail made of only gas, or of dust, or in the heads 
of some comets which have a very weak continuum. 
Excellent broad-band observations have been made, 
for instance, by Bappu and Sinvhal (1957), who have 
used interference filters, by Boiko and Kharitonov 
(1957), by Hoag (1958), by Johnson (1960), and by 
Thiessen (1958). Their observations have given us a 
fairly clear picture of the brightnesses, colors, and the 
percent polarization of the two bright comets of 1957. 
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It is particularly interesting to note that Thie) 
found a close correlation between the brightne 25 
Comet 1956 / and solar activity, in direct contradic i 
to the results of Beyer (1952), who has been carn 
out visual observation for many years at the si 
observatory. 

My own work on comets has been almost exclusiyjy 
in the field of narrow-band photoelectric photom 
The instrument used was a photoelectric spectrophoti|) 
eter (see Liller 1957) located at the focus of a 244) 
telescope at the University of Michigan Observatc} 
Comets studied so far include 1956 # (Arend-Rolar|} 
1957 d (Mrkos), Comet Encke, and 1959 & (Burnhai\? 

The remainder of this paper is divided into th\} 


of comet head spectra; part II describes briefly cor 
tail observations; and the last section discusses mi of 
chromatic isophotal contours made of Comet 1959 i 


I, HEAD SPECTRA 
A. Time Variations of 1956 / 


The discovery of Comet 1956 / well in advan 
perihelion passage made it possible to observe 
changes of the spectrum of this comet as it approa: 
the sun, passed perihelion, and then receded from t 
sun. Figure 1 shows a spectral scan made shortly af 
Comet Arend-Roland passed perihelion. As was ¢ 
pected, the intensities of the emission bands of C2 
CN relative to the continuum varied, in general, 


Fic. 1. Photoelectric scan of 1 
spectrum of the head of Comet 1956 
(Arend-Roland). A yellow filter wi 
inserted in the light path at way 
lengths greater than \5000 to eliminat 
effects of overlapping orders. Sca 
made at the focus of the rti 
Schmidt telescope, Portage Lal 
Michigan. F 


OBJECT: COMET 19574 
DATE: 8/15/57 
TELES 


SLOT WIDTH: 198 A 
SCAN SPEED: 540 A/min 


3. 2. Photoelectric scan of the red 
nfrared spectrum of the head of 
et 1957 d (Mrkos) made at the 
of the Curtis Schmidt telescope. 
lla, dotted line; comet, dashed 


| Ups 
iely as the heliocentric distance. In units of contin- 
i brightness, the intensity of C2 (A5165) increased 
toximately as the 1.31 power of the heliocentric 
ance. However, it was particularly interesting to 
that C2 did not reach maximum intensity, compared 
he continuum, until almost exactly one month after 
helion passage. 


B. Variations of Spectra with Comet 


‘he molecular emission-band spectra of the four 
iets which have been observed with the photoelectric 
etrophotometer are quite similar. In the photographic 
on of the spectrum, the bands of C2, CN, and C3 
‘most conspicuous, and appear with approximately 
same relative intensities. The sodium D lines were 
erved in Comets 1956 # and 1957 d for a brief period 
ime near perihelion. 
‘he greatest difference between the head spectra of 
se comets is the strength of the underlying contin- 
n. Comets 1956 # and 1957 d were similar in that 
y had relatively strong continua. For example, at a 
locentric distance of 0.68 a.u., the peak intensity of 
45165 (C2) band was just twice the intensity of 
continuum at the same. wavelength. For Comet 
cke this intensity ratio was more than 5 at r=0.68; 
jas over 16 for Comet 1959 & at a heliocentric distance 
).89 a.u. This last value is the only surprising one, 
se it has been fairly well established that comets 
ich have passed close to the sun many times are less 
isty”’ than are very long-period comets. According to 
segawa (1960), the orbit of Comet 1956 h is hyper- 
ic, while the period of revolution of Comet 1957 d is 
000 years, according to Schrutka (1961). Presum- 
y, both of these comets are comparatively ‘“‘young.” 
the other hand, Comet Encke is an ‘‘old” comet with 
eriod of 3.3 years, and it has probably made over 100 
olutions about the sun. The observations of Comet 
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1959 k would, therefore, lead one to believe that it is 
very “old.” It will be very interesting to learn the ele- 
ments of the orbit of this comet. 


C. Spectrum of the Head of Comet 1957 d 
in the Infrared 


Shortly after the discovery of Comet 1957 d, the 
writer obtained two scans of the spectrum of this comet 
out to a wavelength of approximately 10000 A. A 
reproduction of one of the original scans and a compari- 
son of it with the spectrum of Capella (Sp. G8III) is 
shown in Fig. 2. As can be seen, the comet spectrum 
is primarily of a continuous nature; the greatest differ- 
ence in brightness between the comet and the solarlike 
comparison star at approximately 1, is possibly due 
to the (0Q—0) band of red system of CN, according to 
Swings (1958). 


Il, TAIL SPECTRA 


A more detailed account of the observations of the 
spectra of the tails of Comets 1956 /# and 1957 d appears 
elsewhere (Liller 1960). In brief, the results were that 
the observed energy distributions in the continuous 
radiations from the tails can be explained by assuming 
the tails to be made of spherules of iron with diameters 
of approximately 0.6 » separated on the average by three 
or four meters. Reasonable estimates of the volumes of 
the tails of these comets, the rates of mass loss, and 
the relative abundance of the other elements leads to 
total masses for these comets of 10'* and 10” g, 
respectively. 


III. ISOPHOTAL CONTOURS OF COMET 1959 k 


Norman H. Schultz of the University of Michigan 
Observatory and the writer were able to obtain several 
monochromatic intensity profiles of Comet 1959 k by 
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sweeping the telescope east and west along one declina- 
tion circle and allowing the comet to move slowly 
northward or southwards. The isophotal contours show 
clearly that the images in the light of \3884 (CV) and 
45165 (C2) are very nearly equal in size if one considers 
contours of intensity measured in terms of the intensity 
of the brightest pomt of that head isophote. Contours 
measured in the feeble light of the continuum of this 
comet showed that this image was from 1.5 to 2 times 
smaller than the images in the light of the emission 
bands at a heliocentric distance of approximately 1 a.u. 

A single intensity profile of considerable interest was 
made on April 28, 1960, when Comet 1959 & was 1.0 a.u. 
from the sun and only 0.2 a.u. from the earth. The diam- 
eter of the \5165 (C2) image appears to have changed 
by nearly a factor of 2 over a time interval of 70 min. 
According to Barber (1960), at about the same time the 
comet was fluctuating by more than a half magnitude in 
brightness, and on the following night, a very bright 
aurora was observed in Michigan. Obviously, additional 
observations of a similar nature are urgently needed to 
make clear the association between solar activity and 
comet appearance and brightness. 
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_DISCUSSION 


Oprx: Concerning the interpretation of color | 
polarization, in addition to the two parameters ¢} 
sidered here (the refractive index and the part) 


ably change the color, and the composition of 
particles which need not be pure substances at all. T. 
there are a few dielectric nonabsorbing particles, wit 


for absorbing particles there are so many different ¢ 
binations that it may be difficult to unscramble the. j 
and tell what particular kind of particles caused 
One can maké a hypothesis for the most probable lee 
position of such dust, but the distribution of the ré 4 
must be taken into account as well as the polarizat 
and color curves. i 
Litter: I might add that a frequency distribution) 
the squared shape was taken in the calculations of 
scatter curves of the individual dielectric and iron p| 
ticles; the variation is by a factor of 2. 
Donn: Another parameter perhaps most difficul 
all to take into account in relation to the interpreta 
of these curves is the shape of the particles. All one 
do theoretically is to consider the cases of sphere 
spheriods. Actually these kinds of particles are not 
kind that it is easy to see forming in a comet. One sh« 
know to what extent one can relate the characterist 
of a randomly oriented group of irregular particles W 
spherical particles. ft 
Herzperc: In the slide for Encke’s comet the 4 
group was practically absent. Is this effect real?. — 
Liter: This is one band that appears somewk 
different from one comet to another. In Comet 195% 
on one scan it was very strong. But whether this 
indeed a night-to-night or a photometer variation shou 
be determined. — 
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‘OR nearly a century a comet was generally visual- 
| ized as a “sand bank.’ Solid particles in a loose 
regate of varying sizes were assumed to carry with 
‘m a certain amount of gas, to provide the activity 
ithe cometary head and tail, under solar activation. 
snerally the mass was considered to possess gravita- 
fie coherence, that is, to be stable or nearly stable 
der the influence of solar tidal disruption. Slow, and 
vasionally rapid, disintegrations of comets were ex- 
ined as arising from solar heating, collisions, or tide- 
‘sing effects, the loose swarm of particles being con- 
ered, in fact, as only in quasi-stable equilibrium for 
jited periods of time. 
Even the requirement of gravitational coherence 
s more recently abandoned by Lyttleton (1953), who 
acluded that “each particle of the comet must de- 
‘ibe a separate, independent orbit around the sun.” 
te sand-bank theory of comets was further strength- 
ed by an ingenious concept due to Levin (1943), 
lo applied the laws of deabsorption of gases to solids 
‘cometary nuclei. By the effects of solar heating he 
counted rather well, quantitatively, for the variation 
gaseous output and brightness of comets as a function 
‘solar distance. 
In 1950 the author demonstrated some consequences 
an icy-conglomerate comet model in which the 
icleus of the comet is postulated as a discrete mass of 
lids which may have been formed from typical inter- 
ellar material at temperatures of the order of 10- 
°K. In the vicinity of the sun, because of solar 
ating, the materials that sublimate profusely in a 
icuum at temperatures of the order of 150-500°K 
try away from the comet with their vapors the less 
latile constituents in aggregates variously sized from 
tremely fine dust to pieces perhaps meters in dimen- 
m. The icy-conglomerate model simulates almost all 
the observed characteristics of comets. No critical 
mparison among various comet models has been made, 
that the icy model is here examined critically. The 
ajor portion of the present discussion is devoted to the 
servational and theoretical evidence relevant to the 
tual existence of a discrete cometary nucleus, and to 
> nature. 


GAS AND SOLIDS OF A BRIGHT COMET 


As Lyttleton has pointed out, direct observation does 
it prove the existence of a discrete cometary nucleus. 
so, the orbital calculations made from cometary obser- 
tions show that these observations do not always 
yply exactly to the center of mass of a body moving 
ecisely under the law of gravitation. 
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INTRODUCTION The fact that comets are observed to lose gas and dust 


to the interplanetary medium is no demonstration that 
a comet has a discrete nucleus other than a diffuse 
volume of solid particles with possibly absorbed gases. 
The hypothesis of a diffuse nucleus, however, sets defi- 
nite limits on the amount of gas that might be held by 
solid particles with the observed surface area. Some of 
the arguments here follow the principles set by 
Vorontsov-Velyaminov (1946) and lead to similar 
conclusions. 

Wurm (1943) calculated that Halley’s comet must 
have contained some 1.510” molecules of C2 near 
perihelion with lifetimes of 10° sec. For the total passage 
this adds up to some 2X 10*4 molecules or 4X10" g. 

Recent calculations by Wurm for the CO* molecule in 
Comet Arend-Roland, 1956 / lead to an estimate of 
some 8X10** CO* molecules in the tail at one time, a 
mass of 4X10” g. The corresponding values for the tail 
of Halley’s comet would undoubtedly come out larger 
in view of its greater dimensions. Although we do not 
know very accurately the lifetime of the COT molecule 
in space, the material in the tail would surely be re- 
placed at least 10 times during a solar approach, giving 
the total loss of CO* during one passage as a minimum 
of 10%* molecules or 510% g. 

The apparent surface area of the solid material at 
the head of Halley’s comet at great solar distances had 
a minimum value of some 10" cm’, that is, an apparent 
radius of 20 km, derived from the assumption of lunar 
albedo (Wurm 1943; Vorontsov-Velyaminov 1946). 

Thus the gaseous layer of C2 lost by the comet during 
its perihelion passage would amount to roughly 210” 
mol/cm? of superficial surface area of the solid aggre- 
gate. If a monomolecular layer requires some 10” 
mol/cm? the effective thickness of the absorbed C, 
layer was originally 210° molecules, or roughly 10~ 
cm. For the later estimate of CO* by Wurm, the cor- 
responding depth of the COt layer comes out 3X10’ 
molecules or roughly 1 cm. 

There is no reason to believe that the several ob- 
served molecules and atoms in the head of a comet are 
the parent molecules originally present in the comet 
nor that CO* is fundamentally the most abundant 
molecule. Hence there is every reason to believe that a 
solid surface layer, centimeters in depth, is a low esti- 
mate of the gas loss in one perihelion passage by Halley’s 
or a comparable large comet. 

That cometary solids could accumulate 10” or even 
10° molecules of gas per square centimeter of surface 
area during a single revolution about the sun, to main- 
tain quasi-equilibrium in gaseous input and output, is 
quite unthinkable. If we assume a solar corpuscular 
radiation of 1000 atoms per cubic centimeter, moving 
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1000 km/sec and impinging upon Halley’s comet con- 
tinuously at this rate for one revolution, we find 10” 
atoms/cm? rev. Such a rate is inadmissibly high, even 
at the earth’s distance. When we allow for many revolu- 
tions (more than 32 for Halley’s comet) the observed 
loss of CO* in all is the order of a meter thickness to be 
applied to the present apparent surface area. 

We must conclude that the solids in a bright, new 
comet must carry a large reservoir of gas, the order of 
10% to 10% mol/cm? and probably very much more. 
The corresponding total original gaseous mass of the 
comet must be at least 10 g, probably a gross 
underestimate. 

Many of the meteoritic solids ejected from comets 
exist for a considerable period of time in orbits about 
the sun; probably they are the major constituent of 
the zodiacal cloud responsible for the zodiacal light and 
the Fraunhofer corona. The Harvard Meteor Program 
(Jacchia and Whipple 1960) has shown clearly that 
visual meteors arise almost solely from comets. Van de 
Hulst (1947) and Allen (1947) estimate the density of 
the particulate zodiacal cloud in the neighborhood of 
the earth at 10-8 to 10-7! g/cm’. 

The writer (1955) has shown that the zodiacal cloud 
can be maintained in quasi-equilibrium by the continu- 
ous ejection of approximately 1 ton/sec of meteoritic 
material, if this material is not destroyed until it spirals 
in close to the sun under the action of the Poynting- 
Robertson effect. The rate of mass replacement required 
is independent of the particle dimensions (excepting 
rather extreme limits), shapes and densities. In fact, 
however, a major fraction of the mass will be eliminated 
by other processes, perturbations by the planets, inter- 
stellar wind, particle collisions, and sputtering by the 
solar corpuscular radiation. Thus the actual injection 
rate necessary to maintain the zodiacal cloud is at least 
10 g/sec, and possibly as high as 10° g/sec, probably 
from cometary debris. 

To date we have no good method of comparjng quanti- 
tatively the amounts of meteoritic material produced 
by different comets. Let us assume, fairly reasonably, 
that the meteoritic material lost by Halley’s comet in 
one perihelion passage is adequate to supply the zodiacal 
cloud for one year. At a rate of 100 tons/sec, the total 
mass lost in one perihelion passage would amount to 
3X10 g. The “new” Halley’s comet, when it made its 
first perihelion passage, may have contained a mass 
more than 100 times as great, and possibly today con- 
tains 10 times this amount, or the order of 10" to 10" g 

Our knowledge of the masses of meteoroids is still 
rather uncertain, but for a meteoric body of moderately 
faint visual magnitude and a velocity of 60 km/sec, such 
as an Orionid meteor which may arise from Halley’s 
comet, the mass may be the order of 0.01 to 0.4 g. If 
one meteor of 3.5 mag. is to strike the atmosphere and 
be visible to a single observer per 10-minute interval, 
the corresponding space density of meteoritic material 
will then be in the range 4X10-°6 to 14X10 g/cm! 
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TABLE I. Mass* estimates for a “new” Halley’s comet. | 


Source of estimate Mass (g) 

C2 in the tail (Wurm) 4x10 

CO* in the tail (Wurm) 5X 1015 

Dust in the zodiacal cloud 10'§ to 10 

Meteor streams 10'6 to 1018 H 
Dust in the tail (Liller) 5X10" to 7X ia 
Total mass (minimal) 10" to 10% 8. ie 
Radius (density =1.3 g/cm*) 3.1to14.6km | : 


* 100 times the mass loss in the recent perihelion passage. i 
for a one-magnitude range in meteoric size. We do |} 
know how to correct from a range of one magnitude| 
find the total space density. Let us adopt a facton 
about 10 times; the total space density arising from | 
cometary stream will then be some 10-** to 10-% g/c} 

If, as a very rough estimate, the cross section o} 
large stream is 1/5 a.u. in each dimension near } 
earth’s orbit, the area is 10° cm? and the total voluj 
10" cm’. i 

Thus, for Halley’s comet the mass of meteoritic nr} 
terial in the stream, if it does indeed produce the Orioy 
and the 7 Aquarid meteor showers, may be the order} 
10'6 to 108 g in total. This mass would then represé} 
a minimum initial mass of the “new” Halley’s com} 
Similarly, for Tuttle’s comet of 1862, Voronts¢| 
Velyaminov calculated a mass of 51016 g from { 
Perseids. - 

Liller (1960) has made a very significant determir 
tion of the rate at which fine dust is lost by Come 
Arend-Roland (1956 #) and Mrkos (1957 d). From t} 
colors, polarizations, brightnesses, and configuratil 
changes of the particulate tails, he deduces, respectivel 
5X10" and 7X10" g as the total initial masses | 
fine dust for these comets, in the general “class” | 
Halley’s comet. These estimates can be compared wil 
those based on gases, the zodiacal cloud, and meteor 
streams, as shown in Table I. The agreement is’ fair 
satisfactory, considering the gross uncertainty in tI 
basic data and assumptions. Furthermore, these est 
mates all concern specifically different cometary mat 
rial. We conclude that a typical, minimum value fi 
the total mass of a new large comet is the order of 10 
g and may easily be 10" g. The corresponding shel 
radii, for density 1.3 g/cm’, are respectively 3.1 an 
14.6 km. These values are in good agreement with t 
present estimate of 15 km for the maximum radius 
Halley’s comet, based on observed magnitude at ual 
solar distances and on the assumption of a very lor 
albedo, and with the broad expecting of the ic 
comet model. 


DIFFICULTIES WITH THE “SAND-BANK” COMET MODE 


Gas Content 


No diffuse sand-bank type or comet model yet pre 
posed can provide the calculated large ratios of gas t 


i) Leoritic material. Suppose that the amount of meteor- 
material is grossly underestimated, or that the 
eous losses are grossly overestimated. If the meteor- 
loss in one period from a large comet, such as 
lley’s, is assumed to be the order of a thousand times 
ater than 10’ g/rev the comet’s minimum radius 
y)m becomes 15 km, comparable to the value from 
| ghtness-albedo determinations at great solar dis- 
_ tices. Since the comet must then be effectively a single 
dy, we have already abandoned the diffuse sand-bank 
‘Tpothesis in making the assumption of such a large 
Ss. 
| The gas-to-dust ratio in observed large comets ef- 
‘tively eliminates the diffuse sand-bank hypothesis 
the cometary nucleus, unless the gaseous loss has been 
erestimated by a factor of at least 100. In this case, 
e hypothesis can be sustained only on the basis that 
e meteoritic material is highly absorptive of the gases 
volved, the order of 0.001 to 0.1% of the meteoritic 
pass. 


Meteoritic Masses in Comets 


For the minimum dimensions of individual particles 
i the diffuse sand-bank hypothesis, let us take the 
_esent mass of Halley’s comet as one-tenth its initial 
‘ass, or 10"? g, and the density of particles as 3.5 g/cm’. 
47 e then find that the average radius of the particles for 
comet with an apparent effective area corresponding 
) a 20-km radius sphere, as observed, is 20 m. Such a 
onclusion as to the size of particles in the sand-bank 
presents no major difficulty, except in the Lyttle- 
on version, which ‘‘is wholly constituted of small dust 
jarticles.” 
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Motions in Comets 
| 


Plummer (1906) has shown that on the sand-bank 
jaodel the secular reduction in period of Encke’s comet 
‘annot be explained by the Poynting-Robertson effect 
jecause the relative changes in the elements e and a 
ire not consistent with expectation. Furthermore, cer- 
jain comets show increases in their periods instead of 
he shortening observed for Comet Encke, e.g., Comet 
Wolf I, Comet D’Arrest, Comet Pons Winnecke, and 
dossibly even Halley’s comet. No process that solely 
‘educes angular momentum, such as the action of a 
‘esisting medium or the Poynting-Robertson type 
sffects can be invoked to explain both decreases and in- 
creases in cometary periods. We can conclude, therefore, 
that the observed secular changes in the periods of 
comets have not yet been adequately explained on the 
sand-bank model. The writer doubts that an explanation 
can be found. 

The existence of negative period changes in comets 
is also attested by the number of comets and meteor 
streams with orbits lying entirely within Jupiter’s 
orbit. Perturbations by the planets could not so reduce 
the aphelion distances for so many objects in the short 
periods of time available. 
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COMET SYMPOSIUM 
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The icy model allows for period changes in rotating 
comets by jet action of the sublimating gas normal to 
the radius vector. Delays in heat transfer, expected 
for nucleus surfaces containing meteoritic material, 
lead to the observed period changes if reasonable physi- 
cal assumptions and masses are accepted. The required 
radii of the nuclei are the order of 1 to 10 km, which also 
allow adequate lifetimes for larger comets. 

Hamid and Whipple (1953) have reanalyzed the defi- 
nitive orbits for 64 long-period comets to derive a 
weighted mean reduction in effective solar gravity of 
5k? /k? = —0.53(+0.10 m.e.) X10-*. If the effective solar 
gravity is not to be reduced by a factor of more than 
—0.53X10~> by solar radiation pressure, the radius of 
the particles must be of the order of centimeters or 
greater for densities applying to ordinary stony or 
metallic particles. 

For the sun-grazing comet, 1882 II, the derived value 
of 5k?/k? is —0.7(+2.1 m.e.)X 10-5. If we adopt three 
times the mean error as an acceptable limit in 5k?/k?, 
the particle dimensions must have been the order of 
a centimeter or the effect of solar gravity would have 
been greater than that observed. 

We must conclude, therefore, from the motions of 
comets that (a) the sand-bank hypothesis has not and 
probably cannot account for the observed secular 
changes in cometary periods, and (b) that the effective 
action of solar gravity indicates that the particles in 
comets probably average at least a centimeter in dimen- 
sion. The first of these conclusions violates the diffuse 
sand-bank model. Both conclusions violate the Lyttle- 
ton sand bank as he has proposed it. 


Sun-Grazing Comets 


Russell (1926) has shown that for the great sun- 
grazing comet, 1882 II, with perihelion passage less 
than one solar radius of the sun’s surface, solar radiation 
would have vaporized individual particles up to 30 cm 
in diameter. Lyttleton (1953) admits that such a comet 
constructed on his model would have vaporized com- 
pletely. He adopts a temperature of 2000°K which, for 
a mean molecular weight of 50, leads to a gaseous ther- 
mal velocity of about 1 km/sec. He suggests that be- 
cause this velocity is so small in comparison to the linear 
velocity of the comet, some 500 km/sec, “any tendency 
of the gaseous material to expand because of its high 
temperature would lead only to a small!divergence away 
from the standard orbit, measured by an angle of the 
order of 2X10~ radian.” He then presents calculations 
by Hoyle to show that recondensation would occur after 
the solar passage at a linear rate of 10-* cm/hr. 

This argument leads to two major difficulties: (a) 
expansion of the gas along the orbital path, and (b) 
the effect of light pressure on the small particles while 
they are in the process of growth. The change in velocity 
AV required to increase the semimajor axis of the comet 
from the initial value of a=83.4 a.u. to parabolic, if 
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applied along the orbit of the comet is 111 cm/sec. Al- 
most half the material would attain an orbital velocity 
in the positive sense exceeding the parabolic and all 
but about 2% remaining would likewise receive an equal 
velocity change in the negative sense along the orbit, 
for which the reduced semimajor axis would be less 
than one-half the original. Thus, on Lyttleton’s 
hypothesis not more than 1 or 2% of the comet could be 
available for condensation in approximately the original 
orbit. 

Tf we accept the condensation process assumed by 
Lyttleton, we see that the particles grew for 50 days 
with diameters less than 10 cm. It is quite obvious 
that these particles would have been moved into hyper- 
bolic orbits while moving under a reduced gravitational 
attraction caused by light pressure. 

There can be no alternative to the argument that the 
processes suggested by Lyttleton would actually trans- 
form his entire sun-grazing comet to the state normally 
associated with the extended coma and tail of a comet; 
here the material has been removed completely from 
the “sphere” of cometary attraction or motion. 

The great comet of 1882 II is one of a group of six 
comets with extremely small perihelion passages and 
comparable periods, strong evidence that a parent 
comet survived such a close perihelion passage at least 
once in the past, and had split up nto a number of com- 
ponents that also endured. No “sand-bank” comet could 
survive such an experience unless the particles compos- 
ing it were originally at least meters in diameter. 


Brightness Variations of Comets 


In general the great bursts in brightness are extremely 
difficult to explain in terms of a “sand-bank” model 
and appear to require some further, but not necessarily 
difficult, development of the icy comet model. The regu- 
lar brightness variations of “old” and “new” appear to 
contradict expectations for the “sand-bank” model, 
while other observed properties of the two classes of 
comets are more or less possible with both models. 

The several limitations or difficulties outlined above 
are definite contradictions to the Lyttleton version of 
the “sand-bank” model. If one assumes the icy comet 
model, none of these difficulties appear, except possibly 
the large cometary bursts in brightness at great solar 
distances. 

The compact “sand-ba: model encounters most 
of the difficulties of the diffuse model. There is real 
difficulty in obtaining the required gas-to-dust ratio 
unless the particles are assumed to be extraordinarily 
small, or endowed with great gas-containing capacity, 
or unless we have greatly overestimated the gas-to-dust 
ratio. Cometary period changes would remain unex- 
plained as would the observed brightness variations 
generally. 

In the author’s opinion, following Schatzman’s (1953) 
calculations, a compact “sand-bank” model would 
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eventually become a single structure, the adjacent) 
ticles tending to adhere. Although the structure m 
never become strong, the nucleus would lose = ‘ 
tremely slowly after a few solar passages and woul | w 
be recognized as a comet. Such nuclei with extrenjy® 
small perihelion distances might act as ober 
case of the asteroid Hidalgo we do observe a mya) con 
ary type of orbit without cometary activity. 

Certain characteristics of a cometary nucleus 
developed below, on the assumption that the nu ¢ 
is of an icy character with embedded solids of metec 
character. 


—~ 


A COMETARY NUCLEUS REQUIRED TO WITHSTAND 
TIDAL DISRUPTION | ee 


What strength of materials in the icy nucleus 
required to prevent tidal distortion and breakage fo! 
comet such as 1882 II, with an extremely small pe 
helion distance? Calculations show that the nuclal 
such comets may possess maximum dimensions of 
of kilometers if their internal tensile strength is 
of dirty (HO) ice near the freezing point, the orde| 3 
10° to 107 dynes/cm?. Large deviations from spheri 
could reduce the maximum allowed dimension be 
bending forces might be introduced. 


THE DISAPPEARANCE AND SPLITTING OF COMETS) 


Mokhnach (1956) has criticized the icy comet me 
on the grounds that in disruption a great deal of crus! 
ing of the ices in the nucleus should result in an ou} 
pouring of material and an unobserved enormous bur 
in luminosity. In the case of a breaking nucleus ‘ 
must postulate the detailed nature of the breakup b 
fore he can draw conclusions as to whether or not a co: 
tary burst should occur. If the nucleus breaks beca 
of rapid rotation induced, for example, by the 
wheel-rocket effect, then one should agree with Mo 
nach that a large area of new surface should be expo 
producing a burst to accompany the splitting of t 
nucleus. 

On the other hand, if one assumes that an irreg 
shaped cometary nucleus sublimates away (like ar 
irregular iceberg melting) and splits under a relatively 
slow rate of rotation, I can see no need for a bursifl 
and brightening to accompany the splitting. : 

In the case of the sun-grazing Comet 1882 IJ a large 
quantity of material certainly constituted the luminous 
“thread” and temporary “beads” after perihelio 
passage. 


ON 


COMETARY BURSTS 


No adequate theory of cometary bursts has yet fe el 
presented. The general concept of collisions between 
asteroidal-type bodies and cometary nuclei is a rathet 
ad hoc type of explanation, and highly improbable. 
The minimum amount of material required for the: 
Schwassman-Wachmann outbursts is not enormous, the: 
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‘er of 10* tons of fine dust. Whitney (1955) has shown 
\t the amount of solar energy that the comet re- 


duce the effect. He suggests some stormg mecha- 
4)m, perhaps a “boiler” type of outburst. 
daser (1955, 1956) has suggested that free radicals 
Sh as OH, CN, etc., imprisoned in crystals such as 
0, NH; and CH, or in solid hydrates, might provide 
je energy for cometary bursts. Donn and Urey (1956) 
tiilarly have suggested that free radicals (originally 
luded by the writer in the icy comet model) might 
triggered by corpuscular radiation, by solar far 
)raviolet or by some other form of high-voltage 
ergy. More recently they have withdrawn this ex- 
‘mation and suggest the possibility that relatively 
jible molecules such as acetylene (CH:CH) might 
‘mbine with other molecules or atoms of the nucleus to 
) ovide considerable explosive energy. 
|A possible explanation for cometary outbursts may 
llow from the icy model without the assumption of 
ze energy in the ices, if the original deposition of 
‘atter in the nucleus is taken to be highly irregular. 
Also, nonradial forces from sublimating material may 
isily arise in a cometary nucleus, producing great 
‘langes in angular rotation, the ‘“pin-wheel” effect. 
configurations can easily be visualized in which the 
‘mtrifugal force established by such a pin-wheel type 
tforce could detach a large or small volume by a mas- 
)ve rupture. 
Thus it is not clear whether the basic principles of 
je icy model need the addition of novel hypotheses to 
ccount for cometary outbursts. 


SOME PROPERTIES OF COMETARY MATERIALS 


ie 
|Delsemme and Swings (1952) have made a very 
mportant contribution to an understanding of the icy 
fomet nucleus by their concept that solid hydrates 
uch as CH,+6H,O may occur in some quantity. The 
roblem of the enormous differences in vapor pressure 
yetween CH, and H.0 is thus resolved in normal comet 
ictivity. Also, the restriction of an extraordinarily low 
emperature required in comet formation to solidify 
SH, is relieved. 

The general and fundamental problems of the precise 
chemical and physical structure of an icy comet nucleus 
ire too complex and too little understood for proper 
reatment at present. In previous papers (1950, 1951, 
1955) the author has discussed some physical and chemi- 
sal aspects of the icy comet model and will not repeat 
these discussions and calculations here. 

One point is of considerable interest ; the photographic 
meteor studies at Harvard have thrown considerable 
ighton the nature of cometary solids. McCrosky (1955), 
ior example, finds that many meteoroids are so weak 
structurally that they crush at a dynamic pressure of 
).02 atm. Also it is clear that photographic meteoroids 
show an enormous range in character. Jacchia’s (1956, 
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1958) work distinguishes between fragmentation or 
fragility of meteoroids, on the one hand, and a function 
of density and luminous efficiency on the other. Each 
meteor stream, representing an individual comet, shows 
individual characteristics which differ from stream to 
stream. Statistically averaged, the characteristics of 
stream meteors are like those of the sporadic meteors, 
practically all of cometary origin (Jacchia and Whipple 
1960). 

The Taurid meteoroids, originating in Encke’s comet, 
are relatively rigid structures, showing practically no 
fragmentation but possessing average mass-luminosity 
properties. Also, there are many Taurid fireballs, show- 
ing that fairly large structures have survived. Encke’s 
comet must have once been a very large comet, for it is 
known to have contributed huge streams of meteors 
of which only a few are now observable (Whipple 
1940). Hence, we are observing material from the central 
part of a very “old” comet, in the Oort sense. 

The Giacobini-Zinner comet, on the other hand, ap- 
pears to be small and probably “new.” Its meteoroids 
are featherlike, both extremely fragile and probably 
of surprisingly low mean density (~0.01 g/cm‘). 

The Geminid meteoroids, of unknown cometary 
parentage, are of average strength but appear to be 
extraordinarily dense (or of very high luminous effi- 
ciency). Their continued existence with a perihelion 
distance of only 0.14 a.u. must be related to their high 
density. Possibly the parent comet was quite large; 
else how could it have had time to attain such a rare 
orbit with aphelion distance only 2.6 a.u.? Since large 
perturbations by the terrestrial planets are unlikely, 
the aphelion must have been reduced by systematic 
period shortening, a slow process. 

The meteoric data suggest that the centers of “‘old” 
comets, in the Oort sense, spew out large and rigid or 
relatively dense meteoroids, while small comets or the 
surfaces of ‘new” comets emit feathery material, weak 
and of low density. The spectra of comets complement 
the picture. ‘“‘Old” comets show primarily the spectra 
of gases and “‘new” comets the spectra of scattered sun- 
light (Oort and Schmidt 1950). Thus the meteoritic 
material near the centers of “‘old” comets must generally 
cohere, so that fine dust is rare, while the unconsolidated 
dust in small comets and near the surface of “new” 
comets predominates over the gases in producing 
light. 

The conspicuous occurrence of C3 molecules in comets 
with gaseous spectra (Rosen and Swings 1953; Swings, 
Fehrenbach and Woszczyk 1957) adds an interesting 
piece to the jigsaw puzzle of cometary nuclei. Can car- 
bon be the first material to condense in a cooling hot 
cloud of gas, and therefore be more prevalent in centers 
of the large comets? Or is the carbon in small comets in 
the form of dust, so that it does not radiate as a mole- 
cule? Some important chemical process must be indi- 
cated, but what? 

An enormous, if not definitive, insight concerning the 
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nucleus and the evolutionary problems of comets could 
be gained by a space probe made to land on a cometary 
nucleus. Cores of the nucleus should be stratified like 
geological sedimentary strata and should give the oldest 
and least disturbed material record of ancient processes. 

The writer wishes especially to thank Karl Wurm 
and William Liller for their kindness in providing 
manuscripts of their important articles on comets prior 
to publication. 
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DISCUSSION 


——: I am surprised that Whipple has been able to 
’ obtain an upper limit to the size of a comet passing 
close to the sun as a result of perturbations, because 
this is a case of Roche’s limit, which depends on the 
ratio of the density relative to the size of the body and 
to the integrated density in a sphere surrounding the 
perturbing body and passing through the center of the 
comet. That is, one would expect to derive, not an upper 
limit to the radius, but an upper or fewer limit to the 
density. 

WHIPPLE: The icy-conglomerate idea was partly 
to give enough internal strength to support the body. 
If forced to depend on gravity alone this kind of nucleus 
might not survive. Whether or not the density were 
below the Roche limit, cracks could cause destruction. 


PR ED WE EP Bs 


“when close to the sun it is questionable whether thi} 


Oprk: When comets have grazed the sun, disrup|)y 
by tidal force is very simply illustrated. Distrit\) 
the mass of the sun over a sphere — to the distz 1) 
of a comet; this density must be 4 or less than thal 
the comet nucleus for the clei to hold togethegl 


1) 
|! 
aly 
\ 


to hold itself together by gravitation under this col 
tion, hence to have high density. i 
Donn: I think Whipple has left out one great | | i 
vantage of his icy-conglomerate model, in that one 
make calculations and predictions of whee one wall f 
expect to happen and do laboratory work on the str) E 
ture of the nucleus. A basic difficulty concerning 
sand-bank theory is that if small grains once evapor||if 
could ever recondense, because the density is so I} ! 
that some of the recombination mechanisms are vé\y_ 
likely not active. It may be interesting to calculi) 
what would happen in the center of a comet like Enck|# 
starting off with a density reasonable for an icy nl | 
cleus. When you build up to about 10-km radius tH 
internal pressure gets high enough that you stél i 
crushing the low-density snow (in this connection thej§ 
are estimates on the packing of snow in glacieri\} 
Beyond 10-km radius the density would increase toil, 
value of a few tenths, and then there would be changif 
due to heating effects also. The preliminary calculatios} 
are reasonable; one can go farther with such calculation|} 
and also carry out experiments. I hope to do this. | 
: Is there a perfectly simple mechanism fd) 
changing the periods because the individual planeijj 
are in different positions at each perihelion passagi |) 
so one could see how one could get changes in angula|} 
momentum and energy, and perhaps an estimate of how} 
many would be thrown out of the solar system, an|j 
how many would be captured by Jupiter? | 
WuirPte: That, of course, is the basis of thinking bj 
Russell, and Newton also, and Oort takes this ‘int 
account. But does the period change after one haj 
taken all the known perturbations into account, ai} 
Miss Roemer queried? This needs further investigation} 
——: Suppose comets originally have gravel-bank 
structure. Would not the particles in that gravel bank 
go “through” one another? Would they not eventually 
kill their kinetic energy and cluster in a compact body 
with the pieces all touching each other? 
WHIPPLE: Schatzmann looked into that in 1953. If 
the dispersion is too great the time is too long. But it 
depends on the starting assumptions whether it will 
happen: I rather think it will. 


_ tal change in close encounters with the planets 
{pik 1951) are further- developed, auxiliary tables 
culated, and the theory applied to comets and aster- 
lal populations. | 

‘The probabilities depend on the relative velocity 
mponents of the particle (comet, asteroid, meteorite) 
len crossing (not necessarily intersecting) the orbit 
‘the planet. Crossing is defined as equality of the 
tliocentric distances. Secular variation of node and 
tihelion lead repeatedly to intersection. In a Jacobian 
ime rotating with the planet (in circular orbit), the 
locity components, in units of the planet’s orbital 
locity, are then 


sf Ges (2-Al—¢)—-Ao} , (1) 
eadial), 
i U,=+[A(1—e) }} cosi—1 (2) 
“ ransversal in plane of planet), 
U,=+[A(1—e?) }} sinz (3) 
“transversal and normal to the plane), 
U={3—2(4 (1—&) }} cosi— A}? (4) 


i) otal relative velocity vector). 

"Here A is the semimajor axis of particle, in units of 
Hie planet’s heliocentric distance, e the eccentricity, 
id 7 the average inclination of particle’s orbit to that 
' the planet. ‘ 

The mathematical expectation per orbital revolution 
‘| the particle, to be called “probability,” of an event 
ith cross section 7s? is 

; p=sU/(r|U,|sini), (5) 
| . . 
here s is the target radius. 
| For physical collisions 

s=S=R[1+2 u/(UR)}, (6) 
‘here u is the planet’s mass in solar units, and R the 
‘lanet’s radius. 

| For 90° average accumulated angular deflection in 
‘epeated encounters 


s=o={B In{ (D*+B)/(R?+B) }}}, (7): 
vhere 
B=16 p?/(mU*) (8) 
nd 
D= (3 u)}. (9) 


ro” is the cross section for achieving isotropic distribu- 
ion of the direction of the U vector. 


*Supported by a grant from the National Aeronautics and 
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Modifications for an eccentric orbit of the planet are 
also used. . 

The angular deflections define the rate of orbital 
change and elimination by ejection to infinity or to a 
crossing with a larger planet. 

Table I gives some results of the calculations. Ob- 
jects crossing Jupiter’s orbit are chiefly eliminated by 
orbital change and ejection, and have a relatively short 
lifetime.. Those confined to the region of the terrestrial 
planets are long-lived, especially those crossing the orbit 
of Mars alone, and are chiefly eliminated by physical 
collision. The Mars asteroids have lifetimes greater 
than the age of the solar system and represent about 
50% of their original population. The other groups con- 
sist of transient objects, their present number depending 
upon the balance between elimination and supply from 
other sources. 

As distinct from physical disintegration, the survival 
time scale of not too small solid remnants of a comet 
depends solely on the orbital elements and is found to 
be very much longer than estimated time scales of ap- 
parent disintegration. It is suggested that residual 
nuclei of some dead comets are not dissolved into me- 
teors but survive as asteroids until they are removed 
by collisions and orbital change. 

Comets could not have condensed during the lifetime 
of the solar system in the region of 50 000-150 000 a.u.; 
Oort’s hypothesis of their origin as asteroids ejected 
from the inner portions of the solar system can be ac- 
cepted. The genetic link between the two classes of ob- 
jects implies that two main ingredients of their material 
may occur in both, although in different proportions: 
I, Whipple’s icy conglomerate; II, compact meteoritic 
material. The second is especially capable of survival. 

Three contemporary populations of stray bodies are 
compared and their balance of exchange and decay esti- 
mated: (1) group M, the asteroids crossing the orbit 
of Mars; (2) group A, the Apollo-type asteroids cross- 
ing the orbit of the earth and other terrestrial planets; 
(3) group C, the comets. Group A, of an average lifetime 
of 108 years, cannot have well existed from the beginning 
and must be currently replenished; cumulative pertur- 
bations by Mars on group M, lifetime 610° years, 
can supply only a fraction of the loss, the rest being 
apparently furnished by extinct comet. nuclei. A two- 
stage supply, first through capture by Jupiter, then 
through perturbations by the terrestrial planets in 
close encounters, appears to be sufficient to feed into 
the space inside Jupiter’s orbit a number of objects, 
equal to the number of same size lost by the Apollo 
group. The main supply comes from comets whose 
aphelia are close to Jupiter’s orbit. 

The statistical estimates agree as to the order of mag- 
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Taste I. Dynamical elimination probabilities and life expectations of crossing objects. J; relative probabilities of elimina | 
J, by collision; J.., by ejection; Jz, by deflection to Jupiter crossing; Jz, by deflection to earth crossing. } 


ERNSPe 


OPIK 


Object Comet 1939 IV Comet Giacobini-Zinner Comet 1866 I 
Vaisala (Giacobinids) Tempel 
(Leonids) 
Crossings Jupiter only Jupiter-earth earth-Uranus 
Lifetime, 10° yr 0.34 0.50 17.0.5 
E 0.026 0.028 0.062 
Ja 0.974 0.972 0.938 
Object Comet 1862 IIT Comet 1942 IT Comet Encke 
Tuttle Vaisala (Taurids) 
(Perseids) 
Crossings earth-Neptune Jupiter-Neptune —. Mercury-Mars 
Lifetime, 10* yr 162. oes 265. 
Je 0.057 0.036 0.975 
Je 0.943 0.964 ee 
Jz ao see 0.025 | 
Object or group Apollo Hermes Icarus Geographos Apollo group el 
(1951 RA) average average) 
(8 objects) (34 obje 
Crossings Venus-Mars Venus-Mars Mercury-Mars earth-Mars Mercury-Mars Mars only! 
Lifetime, 108 yr 64. 39. 165. 152: 102. 6450. 
Je 0.884 0.926 0.976 0.781 0.871 0. P 
a 0.116 0.074 0.024 0.219 0.129 see Ff 
B =o sure See — ee 


nitude if asteroid diameters (d, km) are calculated with 
lunar albedo, 


logd=3.63++logrA—0.2m, (10) 


where m is the apparent magnitude reduced to full 
phase, r the heliocentric, and A the geocentric distance 
in a.u., and those of comet nuclei from a formula based 
on two observed cases, 


logd=2.18—0.2me, (11) 


where mo is the integrated magnitude reduced to r=1, 
A=1 and a 3-inch telescope (Bobrovnikoff-Schmidt 
reduction). As an independent check, the same formulas 
yield the correct predicted number of meteor craters 
in the lunar Mare Imbrium (Opik 1960). * 
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DISCUSSION 


Waxrete: I have a question for Dr. Opik about the 
Apollo group of asteroids. Of those that have been ob- 


served, I believe that five out of seven or four out\} 
seven are highly selected objects; they’re the onesje 
which the motion is such that when they come near 
earth their angular velocity across the sky is rat 
small—in meteor terms their elongation is near 90°.(f 
is quite possible that the probability of encounter w 
those asteroids would be higher: In other words, havi : 
this motion the perturbations by the earth will act me i 
on them than on the others. I wonder whether calcu 
tions might be carried out for the unknown remaind 
(at least the two or three others that don’t come in 2 
this category); perhaps most of the ones not observ. 
are in the second category, and do not have this” pre 
erty. It might be that the probability will fit in ie 
with as observations on the Mars asteroids. — oY 


Apollo 64 million ue (there are three or four of the ; 

and the others fall nearer Comet Encke with times |) 
about 150 to 200 million years. In the probability . 
discovery I allowed for the relative disproportion ¢ 
the photographic plate, so that the photograph} 
impression is not inversely proportional to the squall 


earth. 
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feNG the last few years, the team at the 
4 §’ Astrophysical Institute of the University of Liége 
l\@ succeeded in obtaining a number of high-dispersion 
‘tra of some relatively bright comets, which show 
varticular the structure of the CH, NH, CN, and 
“bands with considerably more detail than on pre- 
5 isly reported spectrograms. Owing to the weakness 
/¥he underlying continuum, the spectrum of the comet 
‘nham (1959 k) was particularly appropriate for a 
_Ncise determination of the position of as many as 55 

ssion features in the 4050 ae aed system of C3. 


la in a forthcoming paper giving details of their ob- 
‘vations, Swings, Fehrenbach, Dossin, and Haser 
61) stress the need of new laboratory investigations 
-Jcerning the Cs bands. 


ily excited in different kinds of laboratory sources 
{ that the main 4050 A band, which undoubtedly 


i f some molecular constants have been obtained. But 
}: rotational analysis of the whole system is far from 


alysis have been so far completely unsuccessful; in 
/t, in the whole region between 3800-4200 A, con- 
# ning in laboratory emission spectra many thousands 
) lines, regular rotational structure can be observed 
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only in the (0,0) band itself and even the position of 
other band heads have been differently reported by 
various authors and do not always coincide with re- 
ported cometary features. 

A few years ago, G. Herzberg observed that the 
strongest cometary band at 4050 A appeared in absorp- 
tion after flash photolysis of diazomethane. The best 
conditions for the development of this absorption have 
been quite recently carefully studied by Herzberg and 
Shoosmith, and excellent spectrograms with very high 
dispersion (about 0.3 A/mm) were obtained a few weeks 
ago. All the cometary features are clearly recognizable 
in the absorption spectrum and the intensity distribu- 
tion is very similar in both kinds of spectra. Regular 
rotational structure is seen not only in the main band 
but also in several other bands; the analysis is now 
under way. The complete vibrational and rotational 
analysis of the whole system will take a long time. In 
fact, there are definite indications of complications due 
to vibronic interaction and to Fermi resonance. The 
only (preliminary) result obtained up to now concerns 
the vi’ frequency, which can be evaluated to about 
890 cm. But we hope that the available laboratory 
spectra will be adequate for a complete analysis and 
might give, together with new high-dispersion cometary 
spectra, a useful contribution to cometary physics. 
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Star counts on six Schmidt plates, covering the Sculptor dwarf galaxy and the surrounding area, provide 
information on the distribution of the brightest stars in the system. Equidensity contours drawn from the 
counts are smooth ellipses with ellipticity 0.35. Density profiles along the major axis do not agree with 
Hubble’s or de Vaucouleurs’ interpolation formulas for luminosity profiles of giant ellipticals. The size of the 
Sculptor system as observed agrees well with the computed tidally limited size caused by the attraction 


of our Galaxy on its outer stars. 


I. INTRODUCTION 


HE distribution of luminosity across giant ellip- 
tical galaxies has been given a great deal of 
attention (most recently by Mrs. Liller 1960). Similarly 
the distribution of light and stars in globular clusters is 
being studied in detail. (King, in preparation). The 
present paper is one of a series which examines the 
distribution of luminosity with radius for objects inter- 
mediate in their characteristics between the giant 
ellipticals and the globular clusters, namely, the dwarf 
elliptical galaxies. 

The first paper (Hodge 1961) examined the structure 
of the Fornax dwarf galaxy by means of star counts, 
the surface brightness of the object being far too low 
for the standard photoelectric or. photographic methods 
of measurement. The present paper gives the results of 
study of a similarly faint and extended object, the 
Sculptor dwarf elliptical. 

Since Shapley’s discovery of the Sculptor dwarf 
(Shapley 1938, 1939), relatively little has been pub- 
lished about its characteristics. Baade and Hubble 
(1939) used the 100-inch telescope to corifirm Shapley’s 
conclusion that the Sculptor system resembled an 
enlarged globular cluster by their discovery of RR 
Lyrae variables in it. Later Thackeray (1950) obtained 
an important and extensive series of plates from which 
he was able to determine periods and light curve 
characteristics for a number of the variables. 

In his announcement of the discovery of the Sculptor 
system, Shapley (1938) tabulated the results of star 
counts on the discovery plate. The writer has extended 
Shapley’s studies by means of star counts on six new 
plates. These have provided fairly complete data on the 
stellar distribution, which may now be compared with 
those for other dwarf ellipticals and for giant galaxies. 


II. OBSERVATIONAL MATERIAL 


The six plates used for the present study were 
obtained with the ADH Baker-Schmidt telescope at the 
Boyden Observatory, South Africa. The Schmidt cor- 
recting plate was diaphragmed to 24 inches to avoid 
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vignetting. Three of the plates were exposed on 103a 
ale Sbebind a~GG-11 filter and three on 103aj 
emulsion behind’-a GG-13 filter; the color system 
equivalent to V and B of the UBV system. 
The exposures represent a range in magnitude lin m 
exact magnitudes are not available, but it can be 
mated that the deepest plate, ADH 5059, reach 
approximately B=19.1. This estimate is based on ti) 
fact that at the limit of the plate lie several of | 
variables identified by Baade and Hubble (1939) ; thi 
estimate of the mean magnitudes of these objects | 
maximum is 19.12 (pg). 
A 2.5-mm grid was used for counting all stars in 4 


system. The grid was oriented VS and EW, 
nearly parallel with the axes of the object 
below). Each plate was counted at one sitting so 
uniformity would be maintained. The plates 
counted to the limiting magnitude and each w 
counted once. From experience with similar star coun 
made on the Fornax dwarf elliptical, it was decid 
that the reproducibility of counts on a particular pla 
is greater than from plate to plate, so that it is bet 
to. count many plates each once than one plate mat 
times. 
In the richest portion of the galaxy the star densi | 
on the deepest plate, ADH 5059, reaches a maxinn 
value of eight stars per square minute of arc. The imag | 
diameter on this plate is less than 0.05 mm=0.06 6 
thus the probability of a star being obscured by anoth¢ 
on this plate is very small, of the order of 0.02. We a 
conclude that saturation of images, as found on one 
the plates used for the Fornax galaxy (Hodge 156 
does not occur in the present material. | 
i 
TaBLe I. Observational data. | 
See SS SSS SSS ee 


Stars counted -Foregroun 
Plate - Exposure Members of (stars | 
No. Color (minutes) Total Sculptor per oO’) 
5109 V 10 2281 254 0.28. 
5088 V 20 2878 565 0.2193 
5111 v 30 2737 509 0.31 9 
5058 B 60 3350 882 234.4 
5061 B 60 3676 1104 0.35 
5059 B 60 4630 1529 0.42 
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able I lists the plates used and the total numbers 
stars counted on each. The number given there for 
: total number of the member stars of Sculptor was 
tained by subtracting the expected numbers of back- 
jund stars, as determined from surrounding more 
| tant areas, from the number counted within the area 
|the galaxy itself. As pointed out below, the Sculptor 
jtem has a definite limiting radius where the star 
‘asity reaches zero; for this reason the number of 
¢imted members for each plate could be obtained 
tambiguously. 


DISTRIBUTION OF STARS 
Lines of Equidensity 


Contour diagrams illustrating the counts for each 
Jite in the form of equidensity lines are given in Fig. 1. 
Jich was constructed directly from the recorded counts 
‘enveloping within a particular curve the squares with 
tparticular count, in steps of five stars per counting 
suare. The contours plotted in Fig. 1 include the 
dreground stars in their counts. The outer equidensity 
le is in each case that for 12 stars per counting square 
5 stars per square minute of arc), and the contour 
iterval is five stars per counting square (0.6 stars per 
{uare minute of arc). 

)The diagrams in Fig. 1 were used to determine the 
jneral properties of the shape of the system: orient- 
ition, ellipticity, and symmetry. The actual profiles of 
jellar density and subsequent analysis are based 
rectly on the counts. 


Tsophotes 


| 
{[t would be of considerable interest to have an 


sophotal diagram of the Sculptor system so as to be 
le to compare its properties with those of galaxies 
ir which luminosity density is directly observable. It 
possible to artificially construct such a diagram if the 
sumption is made that the luminosity function of the 
tulptor system is similar to that of M3. We know from 
aapley’s (1938) work that this is probably true, at 
ast for the top 1.5 mag. (pg). Similarly, Baade and Miss 


on of the Draco dwarf system also resembles that of 
[3. The fact that Baum (1959) finds that dwarf ellip- 
cals (unlike giants) have integrated colors similar to 
1ose of globular clusters is further evidence that the 
art of the luminosity function responsible for most of 
1e light is probably similar to the normal globular 
uster luminosity function. 

To convert star densities into isophotes, it is necessary 
) know what fraction of the total luminosity of the 
7stem’is contributed by the unresolved stars in the 
ystem. For this purpose, we can use Sandage’s (1954) 
dmputation of the fractional light for stars at M,, and 
ighter in M3. The depth of the counts in Sculptor on 
late ADH 5088, the deepest visual plate, is approxi- 


’ 


SCULPTOR DWARF GALAXY 


wope (1961) find that the top of the luminosity func- 


Fic. 1. Equiden- 
sity contours drawn Si! ’ 
from the star counts. 
The outer contour 
line represents a pro- 
jected density of 1.5 
stars per square 
minute of arc. The 
contour interval is 
0.6 star. 


mately 1.6 mag., as judged from photoelectric sequences 
on similarly exposed ADH plates. Thus the counts 
reach M,,~—1.4, and the counted stars account for 
0.28 of the total light of the system, if our assumptions 
are correct. 

Figure 2 shows the artificially constructed isophotal 
map of the Sculptor system, based on the star counts 
on plate ADH 5088 (smoothed) and Sandage’s lumi- 
nosity function. The map indicates that in visual light 
the luminosity density at the center of the system 
reaches 23.9 mag. per square second, while at a distance 
of 30’ from the center along the major axis this value 
has dropped to 27.7 mag. per square second. 


Ellipticity 
The ellipticity of the Sculptor system, defined by 
e=1—b/a, (1) 


has been derived by measuring it for all curves in the 
raw isodensity diagrams. The mean value for all 39 


Fic. 2. Artificially constructed isophotal map of Sculptor, based 
on assumptions of a globular-cluster-like luminosity function, 
invariant across the galaxy. Contours are labeled in units of 
magnitudes per square second. 
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Taste II. Profiles along the major axis.* 
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a Average number of stars per counting square (8 D’). 


curves is 0.3540.05. A plot of ¢« against r shows that 
there is no detectable variation of ellipticity with 
distance from the center, as had been found for the 
Fornax dwarf (Hodge 1961) and for giant ellipticals 
(Liller 1960). 


Orientation 


From measures of the orientation of all the curves 
in the isodensity diagrams, the orientation of the 
Sculptor dwarf is found to be 98° from the north, 
measured eastwards. This value has a probable error of 
about 2°. 


Symmetry 


Unlike the Fornax dwarf system, the Sculptor 
galaxy seems to be perfectly symmetrical. Neither from 
the isodensity diagrams nor from the detailed profiles 


a 


Fic. 3. Major axis profiles for the six counted plates. The projected 
density is in units of stars per square minute of arc. 
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is there any suggestion that the object is not a perfect}: 
symmetrical ellipsoid. | 


Density Profiles 


i 
4 
Table II lists the star densities, determined for di), 
ferent distances along the major axes of the system), 
for each of the six plates. The values given in Table I}, 
are averages of the number of stars per counting squar}) 
for the two or three squares intersected by the maja|) 
axis and/or a line 3’ long perpendicular to the axis ani}, 
centered on it at at distance r from the center of th 
object. The values in Table IT include foreground stars) 

The east and west halves of the major axis ar! 
averaged together in Table III, which gives profiles i} 
terms of stars per square minute of arc, with the fore 
ground stars (Table I) subtracted out. Figure 3 shows ¢| 
plot of these points, with smooth curves drawn through} 
them by eye. 1 

The six curves can be combined to give a meat 
curve only if the luminosity function and the color: 


Fic. 4. The difference between the normalized density for 
each plate and the mean density at that point. 
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Fic. 5. The outer mean contour points, showing the adopted 
| curve drawn through them. 


e center of the system. To test this, the profiles given 
|Table III have each been multiplied by the appro- 
jiate factor to bring the smoothed curves in Fig. 3 
gether at r= 6’. The mean of these normalized profiles 
is been taken, giving the three deepest plates double 
eight. Finally, the normalized points NV for each plate 
.we been subtracted at each given distance r from the 
ean profile M. The difference M—VN is plotted in Fig. 
as a function of r. The figure shows that there is no 
dpreciable slope for the lines of any of the plates. Thus 
e conclude that the tops of the luminosity function 
ad color-magnitude diagram are sufficiently invariant 
ith distance that we may have confidence in a mean 
lajor axis profile derived from all six plates. 

A mean profile, computed as described above, is 
lven in the last column of Table III. The tabulated 
alues are the computed ones out to a distance of 
=15’, beyond which they are smoothed values read 
om a curve drawn through the computed points (Fig. 
). This profile, derived from all six counts, is considered 
le most representative of the galaxy, and is therefore 
sed in all of the following analysis. : 


| << diagram are invariant with distance from 


Radius 


In the case of the Fornax system, it was found that 
1e galaxy had an actual radius, determined by the fact 
iat the stellar density reached a definite zero value. 
his radius was read from a curve showing the density 
lotted against the reciprocal of the radius, as done for 
1e globular cluster M15 by King (1961). Such a plot 
f the outer parts of the Sculptor galaxy is shown in 
ig. 6. The radius of the system derived from this plot 
r=46/+3' for the major axis. This means that if the 
istance modulus is m—M=19.6 (Baade and Hubble 
939), then the linear radius of the Sculptor dwarf 
long the major axis is r=1.1 kpc. 


Shapley’s Star Counts 


In his initial paper on the Sculptor system, Shapley 
1938) gave a table of counts made on the discovery 


DWARF GALAXY 


387 
25 


-20 


Fic. 6. The mean density distribution plotted in 
units for the radius determination. 


plate, A18005, for an area of 4° centered on the galaxy. 
The plate used was of very unusually high sensitivity 
so that the limiting magnitude of the counts was 
approximately 19.5 pg. In order to compare Shapley’s 
data with that of the present work, a major-axis profile 
was computed from his counts (Table IV). First the 
background count was established by averaging the 
counts from 70 squares in the four corners of his grid, 
where the squares were all at a distance appreciably 
greater than the radius derived above (46’ for the 
major axis and 30’ for the minor). The background was 
found to be 0.58 stars per square minute of arc. Table IV 
gives the star density along the E-W axis, which is 
very nearly the major axis, obtained by averaging the 
two sides and subtracting the background. The table 
indicates that either foreground fluctuations or counting 
irregularities happen to depress the foreground count 
along the major axis in the very outer two steps of the 
grid, so as to make the major axis look larger than the 
value obtained when a reliable foreground count is used 
to determine the zero point. This would explain Shap- 


TaBLeE III. Major-axis profiles and mean profile. 


Zee D109) 5088) Sk 5058 5061 5059 Mean 
Stars per square minute Profile 
0’ 1.06 2.47 2.80 4.12 4.67 (ieehes ioe: 
3 0.83 BBS 1.98 3.91 4.25 5.80 5.94 
6 0.78 1.97 1.34 2.69 3.46 4.49 4.70 
9 0.65 0.94 1.07 1.53 Me), 3.14 3.12 
12 0.30 0.72 0.78 1.02 1.48 2.41 2.09 
15 0.35 0.48 0.55 0.84 0.98 1.77 1.47 
18 0.25 0.35 0.79 0.58 0.56 0.94 1.13 
21 0.10 OTs 0.98 0.29 0.19 0.44 0.83 
24 0.14 0.07 —0.01 0.16 0.26 0.19 0.58 
27 O19N 10,295) 05140512) 0:26 #0228; 710.38 
30) mont haa. 1 OFT pee, OL 2d Onl 8 O13 0226 
33 —0.01 MOS Ost ae Os Om 0 T1ON Oto 10.16 
36 —0.03 0.05 —0.05 0.02 0.10 0.06 0.09 
39 —0.06 0.14 0.02 0.05 —0.02 0.09 0.04 
42 —0.10 0.22 —0.01 —0.05 0.08 0.04 0.02 
45 0.11 0.27 —0.01 -—0.15 0.02 —0.02 
480 OS tie 0515) 0712 0.18 
51 0.09 0.11 
54 —0.14 —0.05 
57 —0.06 —0.05 
60 =0.03 
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ley’s statement that the diameter of the Sculptor system 
might be as large as 2°. 
Figure 7 shows the outer portion of the curve derived 
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Fic. 7. The outer density distribution from Shapley’s counts. 
The background was established from larger areas at greater 
distances than shown here. 


from Shapley’s counts. The curve drawn smoothly 
through these points has been used to find the radius 
of the galaxy to compare it with that derived above 
from the present work. Figure 8 shows the curve used, 
which gives a value for the radius (measured 8° from 
the major axis) of r=40’. This value is very reasonable 
when compared with 46’ obtained above for the major 


TasLe IV. Major-axis profile after Shapley. 


Tr 


p 
Minutes of arc Stars per D1’ 


0’ 12.70 
5 9.76 
10 6.28 
15 4.38 
20 3212 
25 1.84 
30 0.84 
35 0.16 
40 0.12 
45 0.00 
50 —0.16 
55 —0.16 


axis, indicating satisfactory agreement between the 
two independent studies. 


Transparency 


The images of galaxies lying in the direction of the 
Sculptor system have been counted on two plates. The 
projected density of these objects is tabulated as a 
function of distance from the center of the Sculptor 
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galaxy in Table V. The data show no apprec| ' 
diminishing of the background galaxy count be 
the Sculptor system, especially when it is taken | 


10 44 
20 41 
30 42 
40 50 
50 55 
60 50 
70 49 


account that a distinct cluster of galaxies lies dire i 
south of the system at a distance of about 50’. 


IV. DISCUSSION tT 
. . } 
In general shape, the major-axis profile of the Sculpip 


Al8005 


Fic. 8. The curve used to determine the radius 
from Shapley’s counts. 


i 
(Hodge 1961), with the exception that the densi! 
gradient of the outer portions of Sculptor is somewhi 
steeper. Figure 9 compares the observed density di 
tribution with the formula 


N=N)/ (r/a+1)’, ( 


where NV is the projected density, No is the centr 
density, r is the distance along the major axis, and 
is a scale parameter, equal to rat N={NVo. This formul 
is Hubble’s luminosity law for giant ellipticals, writte 
in terms of star density on the assumption that th 
luminosity function is constant throughout the galaxy 
The observed curve does not fit the computed one, an 
adjusting the values of a or No does not improve the fi 
As in the case of Fornax, the curve fits the predictio 
moderately well out to a distance of 1¢ (13), after whic 
it falls off much more steeply than the prediction. 
Similarly the Sculptor dwarf does not conform t 


_|Vaucouleurs’ law (de Vaucouleurs 1948). For his 
te ; oat Mee ae ee, 

-jaula a galaxy’s luminosity distribution forms a 
ight line when plotted in the units of Fig. 10. The 
ve there cannot be considered a straight line. 


Gravitationally Limited Radius 


: 

|| ae 
‘he similarity of the Sculptor with the Fornax system 
vests that its radius, too, may be determined by the 
vitational attraction on its outer stars by the 
axy. The expression for the gravitational limited 
ius is (approximately) 


r= R,(M/2M,)* (3) 


n Hoerner 1957, Burbidge and Sandage 1958), where 
the tidal radius, R, is the distance to the Galaxy, 
is the mass of Sculptor, and M, is the mass of the 
axy. 

tis necessary, before computing the gravitationally 
ited radius, to have an estimate ofthe mass of the 

tem. In the case of Sculptor the mass cannot be 

asured directly. Instead we are left with two possible 

ans of obtaining a rough guess at its mass. 

first, we may compare its luminosity function with 
t of a globular cluster of known mass. Such a com- 


computed 


is 


a. 
observed \ 


x 


' 

| 

i 

} 

‘| 
i log (r+a) 


“4 9. The density distribution plotted in units in which Hubble’s 
law (dashed lines) gives a straight line. 


the two kinds of systems. Shapley’s luminosity 
jaction (1938) indicates that Sculptor is 15.2 times as 
h as the globular M3, if we compare the two in the 
‘erval between the brightest stars and 0.5 mag. 
low the horizontal branch. This is necessarily uncer- 
in because of the lack of knowledge of the similarity 
dissimilarity of the two systems. However, the com- 
rison ‘can give an estimate of the mass of Sculptor, 
ing Sandage’s mass for M3 (1954). The result is 
=2X10° Mo. 

Another estimate can be made from Sculptor’s 
ysolute magnitude, assuming a mass-to-light ratio for 
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jrison assumes equality of the luminosity functions ~ 
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the system. Shapley (1939) gives its apparent magni- 
tude as approximately 9.0, photographic. If its mass-to- 
light ratio is the same as for globular clusters, of the 
order of 2 (Sandage 1954), then the mass of the Sculptor 
dwarf is M=4X10° Mo. A larger mass-to-light ratio 
leads to a proportionally larger estimate of the mass, 
but enlarges the estimated limiting radius only by the 
cube root of this increase. Using Brandt’s (1960) value 
for the mass of the Galaxy, Baade and Hubble’s (1939) 
value for the distance of Sculptor, and 310° Mo as 
the mass of Sculptor, the computed tidal limit is 
v=1.7 kpc. This is slightly larger than the observed 
radius of 1.1 kpc, and is comparable to it considering 
the uncertainties involved in the mass estimates and 
considering the fact that the tidal radius applies to the 
situation of nearest approach of the system to our 
Galaxy. Uncertainties in the distance are large, but 
not very important in this connéction, as they tend to 
cancel each other out in the comparison of the two 
radii. 
V. SUMMARY 


The Sculptor dwarf galaxy is very similar in nearly 
every respect to the previously studied dwarf in Fornax. 
Its stars are smoothly distributed in a symmetrical 
elliptical configuration of ellipticity 0.35. The distri- 
bution of star density with radius does not agree in 
shape with that of giant elliptical galaxies. The radius 
of the Sculptor dwarf is 46’, approximately 1.1 kpc; 


Fic. 10. The density distribution plotted in units in which 
de Vaucouleurs’ law gives a straight line. 


this agrees well, within the uncertainties, with the 
computed tidal radius imposed on the system by the 
mass of the Galaxy. 
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Multicolor Photoelectric Photometry of Bright Galaxies.* I 
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Photoelectric observations are presented for 57 galaxies and 52 stars as obtained with a blue-sensitive 
photocell. The response bands have effective reciprocal wavelengths of 2.66, 2.39, 2.07, and 1.68 (micron 
units). In addition, 8 galaxies and 19 stars have been observed in the red and infrared. The effective 
reciprocal wavelengths for the red response bands are 1.55, 1.26, 1.14, and 1.00. From one to five aperture 
sizes have been used in the photometry of each galaxy, with an average of three per galaxy. All the obser- 
vations have been reduced to an absolute energy color system by means of direct comparisons between the 
response of the photocell and filter combinations and the response of a thermocouple. A subsequent paper 


will contain the bulk of the discussion of the data. 


INTRODUCTION 


UMEROUS investigations involving extragalactic 

photometry have been carried out; among the 
most extensive investigations are those of Stebbins 
and Whitford (1952), Pettit (1954), Holmberg (1958), 
and de Vaucouleurs (1961). Some six-color observations 
were obtained by Stebbins and Whitford (1948); 
however, most of the available photometry on extra- 
galactic systems has been limited to two colors and, 
in the case of photoelectric work, to only one or at 
best relatively few aperture sizes for each galaxy. This 
present paper extends photoelectric observations, 
similar to the Stebbins and Whitford six-color data, to 
a larger sample of galaxies with up to five aperture 
sizes used in the measurement on each system. Some 
measurements have been made in eight colors using 
two photocells; however, the vast majority of the 
observations were taken in four colors using a blue- 
sensitive photocell. The results of the work are pre- 
sented in two papers; this paper includes the obser- 
vations and describes the photometric system and 
procedure, while a later paper will contain a discussion 
of the observations. 


INSTRUMENTATION AND TECHNIQUES 


Photoelectric Equipment 


The blue sensitive photometer employed utilized an 
unrefrigerated type 6094 EMI photocell. The 11-stage 
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end-on EMI photocell was normally operated with| { 
total applied voltage of 1300 v, and the output fi), 
through a dc amplifier and recorded on a Brow 
recorder. Red and infrared observations were obtaini) 
with a type 16 PM-1 Farnsworth photomultiplier tul}) 
The 16-stage end-on photocell was operated und|) 
dry-ice refrigeration with a total applied voltage |} 
2200 v from an electronic power supply. The first of} 
servations were obtained with the Farnsworth phot . 
in conjunction with a chopper stabilized amplifi¢ 
Later observations were taken with an electromet 
tube amplifier. As with the EMI photocell the amplifi 
output was recorded on a Brown recorder. | 

The sensitivity curves of the EMI and Farnswor') 
phototubes were determined as a function of waveleng#} 
for use in absolute energy calibration of the coli} 
system. The calibration was carried out by comparit} 
the response of the phototubes with that of a thermi) 
couple. The source of illumination was a tungsten lam} 
viewed through a Hilger-Muller double-pass quari} 
monochromator. Both the thermocouple and photoce} 
entrance windows were placed in the same position, ¢| 
close as possible to the exit slit of the monochromato} 
Because of the great sensitivity of the photocell} 
compared to the thermocouple, it was necessary t 
introduce light attenuation, in the form of diaphragm) 
in the light beam before it entered the monochromatc} 
when the photocells were used. The thermocoupl} 
output was read with a sensitive galvanometer, an 
the photocell output was fed through an amplifier an’ 
recorded on either a Brown or Varian recorder. 

The calibration was carried out in three sessions 


}: Farnsworth photocell under dry-ice refrigeration 
‘tween A5500 and A16 000, the EMI photocell between 
»500 and 8000, and the EMI photocell between 
? 500 and 4900. The photocells were operated at 500 v 
‘xen from an electronic supply. Throughout the cali- 
J ation the entrance and exit slits of the monochromator 
‘re set equal and the center slit at about twice the 
+ dth of the other two. For the Farnsworth and longer 
jivelength EMI calibration the exit slitwidth corre- 
jonded to 240 A in the emergent spectrum at \5850. 
/ wr the ultraviolet EMI calibration the exit slitwidth 
i rresponded to 260 A in the emergent spectrum at \4045. 
“i.e measuring procedure consisted of scanning the 
tire range of interest with the thermocouple, then 
_peating the scan with the photocell, and then repeat- 
'g the thermocouple scan as a check. The two thermo- 
‘uple scans agreed satisfactorily in all cases. The 
-avelength scale of the monochromator was calibrated 
', visual observations and photoelectric scanning. 
ossible errors in the wavelength calibration range 
‘om 5 A at 44500 to as much as 50 A at 410000. 

In view of the technique employed, the only impor- 
‘mt differences in the light seen by the thermocouple 
‘ad photocell arise from the light attenuation dia- 
-aragming, and slight differences in the viewing at the 
“sit slit, particularly the viewing of light scattered from 
e jaws of the exit slit. Both of these effects are prob- 
ly small. Visual inspection of the emergent light 
‘iggests that a small amount of light is scattered into 
ie beam from adjacent spectral regions, probably by 
“vattering at the jaws of the central slit. The effect 
‘ppears to be small and in any case the wavelength 
“ifference of the light from the main bundle is small. 
‘hotoelectric scanning in the ultraviolet beyond the 
\hototube cutoff and in the far infrared gives zero 
esponse, implying that scattering over large wavelength 
|ifferences is virtually zero. Table I gives the measured 
‘hototube response relative to the thermocouple as a 


i 

| 

| & 

TABLE I, Normalized response of photocells. 
| 


1/r Wg 1/r ie 1/r P 
EMI Photocell 
moo .626 2.45 1.000 2.00 ~ .707 1.55 .00790 
fees .699 2.40 .999 1.95  .632 1.50 .00234 
Pee sOseueg oo 6,904 “91.901 ..550 1.45 .00059 
eee coteeeso. 980) 1.85. 2447 1.40 .000157 
1.70 (SOU mee eos 962 1,80 £328 1.35 .000033 
feet. O1Qenz ON 933° 1.75.- .222 1.30 .0000068 
moo 659460 2:15 .895 1.70 .144 1.25 .0000012 
os 1,974 2.510.839 1.65 .0673 
50.993. 2.05 776 ~=1.60 0224 

Farnsworth Photocell = 

MeO) 329%) 991-40.- .587 1.10- 1.000 0.80 .00275 
Mooroloneooe! 675 )06=— 1.05. .953 0.75 .00013 
OO 7.550). tc00 2761 1.00  .767 0.70 .000015 
foo .o97F- 1:25 .840 0.95 .468 
1.50 .448 1.20 .909 0.90  .159 
45512 93515 2966 0:85 .026 
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TaBLeE II. Filter system. 


Band Filters 


Corning 9863--approximately Corning 7380 

Schott BG-12+-Schott GG-13 

Corning 5030 (1.6 mm)+Corning 3387 

Corning 3480 

Federal Engineering special interference filter 

F we Engineering special interference filter+Corning 


ei Engineering special interference filter+Schott 
-10 
Heimann No. 205 
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normalized to read 1.000 at maximum sensitivity. The 
values are not corrected for the effect of finite slitwidth. 


Filter System 


Corning and Schott glass filters or filter combinations 
were used to define four response bands in conjunction 
with the EMI photocell. Special interference filters 
made by Federal Engineering (Madison, Wisconsin) 
were used for three of the four red and infrared pass 
bands. For the far infrared a Heimann filter was used. 
Clear glass was cemented to three of the four glass 
filters, and the far-infrared filter, to bring them all to 
the thickness of the thickest glass filter. In the case of 
the ultraviolet filter the added glass cuts out light 
below A3400 normally transmitted by the filter. The 
interference filters were not adjusted in thickness but 
are nearly the same as the others. Table II lists the 
filter components and thickness when different from 
stock values. The filters and response bands they define 
are designated by serial numbers in order of wavelength. 
This numeral notation is used throughout the paper. 

Filter transmissions were determined as a function 
of wavelength with a Hilger-Muller double-pass quartz 
monochromator by viewing the emergent beam alter- 
nately through a filter and directly. The measured filter 
transmission F, combined with the photocell response P 
from Table I, and the reflectivity of aluminum MM, 
according to Allen (1953), define the response function 
G. Using the reciprocal wavelength units, 


G(x)=P(«)F(x)M?(~), wx=1/d. (1) 


The response functions include two aluminum reflec- 
tions, since all observations were made at the Newtonian 
focus of the Mount Wilson reflectors. Table III gives 
the eight G(«) functions; Fig. 1 illustrates the response 
functions. 


Observing Program and Technique 


All observations were made at the Newtonian focus 
of the 60-inch and 100-inch telescopes at Mount Wilson 
during 1957 and early 1958. All red observations were 
made with the 100-inch telescope, while the EMI 
photocell was used at both reflectors. Five circular 
diaphragms were used in the program, from 0.8 to 7.0 
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TaBLE III, G(«) functions. 

1/r Gy 1/r Ge 1/r G3 1/r Gs 
2EO 0a 007, Dajan we OOD 230") 007, 1.80 .002 
2289" 2050 2-40) 056 2.25 .200 densi O41 
2.80 .166 24655 eelo9 2,20) 2356 Lef0r) 7.082 
Dio «292 2.60 .240 2.457 41S 1.65 043 
DOr aor0) 2255 .338 2.10 .418 1.60 .014 
2.65 .388 2.50) ~ 7o9Z 2.05 ~ .384 1555 2005 
2.60 .346 2.45 .419 Deo eaZo 12507 + 500% 
PAROS) CLPAY 2.40 .428 1.95 200 
2,00" 1093 Pp hile) 1.90 .174 
2.45 .025 2.30 ~ .384 1.85 . 102 
2.40 .006 B25 .336 1.80 .047 
Desor) 002 Pee AA if LoS .017 

Plbaa Role” ee fest 110 2005 
2105 92083 
2.05  .019 
2.00  .003 

ae ce UX? Gs iy Gs Ue een 
1.700 .007 1.450 .004 1.350" 3001 1.150 .004 
1.675 .064 1.425 .028 15325." 3003 ia Uy¥L 
1.650 .145 1.400 .093 1.300 .010 1.100 .103 
150255) AGL 1.375 .198 1,275) 2029 A075: eeoo 
1.600 .189 1-350) 2307 1.250 .080 1.050 .369 
LSi5 eke 1.325) | 3360 L225) 84 1.025 .446 
1.550 .149 1.300 .371 £200) .312 1.000 .441 
1.525 .148 1.275 .404 LAS) ere AO 0.975 -.384 
1.500 .199 12250; 2455 L50> 2385 0.950 .298 
1.475 .208 1.225 .486 PAST Fst hs: °0.925 .198 
1.450 .087 1.200 .453 1.100 .250 0.900 .107 
1.425 .024 7S: 312 1.075 .242 O78i5, .052 
1.400 .009 ISO eno, 1.050. .274 0.850 .018 
1.375 .004 SASy (05:5) 150255270 0.825 .008 
1.350 .002 1.100 .018 1.000 .043 0.800 .002 
1625) 002 1.075 .007 0.975 .010 
1.300 .001 1.050 .003 0.950 .002 
1.275 .001 1.025 .002 
de250) 001 1.000 .001 
1.200 .001 0.975 .001 
1.150 .002 
1.100 .006 
1.050 .007 
1.000 .004 
0.950 .002 
0.900 .001 


mm in diameter in steps of about three in area. The 
diameters, measured in a measuring engine, ‘are given 
in Table IV along with their actual angular size in the 
focal planes of the 60-inch and 100-inch telescopes. 
Galaxies were selected for observation according to 
the following restrictions: (1) With the exception of 
M31, M32, M33, and NGC 205, all galaxies lie between 
+30° and +90° galactic latitude. This restriction was 
put on the selection to reduce effects of galactic ab- 
sorption. Proper calibration of the galactic absorption 
effect would have required observations of a much larger 
sample of galaxies. (2) The majority of the galaxies 
observed are bright noncluster field galaxies with small 
redshifts. A few representative galaxies from the Virgo 
and Coma clusters were also included. (3) All the Hubble 
extragalactic types are represented, including examples 


over the full range of inclination from edge-on to, 


face-on. (4) Whenever possible, galaxies were taken 
from the list of systems for which Morgan and Mayall 
(1957) have given spectral types. To complete repre- 
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Fic. 1. Response functions for the photometric system. 
sentation of structural types, additional galaxies a a 
selected from the Humason, Mayall, and Sandaj), 
(1956) redshift catalogue. (5) With one accidental el) 
ception, all galaxies are free from interference by | ‘ 
field stars. i 


observed to determine stellar colon on the Shotonl i 
system. The majority of the stars observed range fro 
spectral type F to M and include representatives of f 
main sequence, normal and high-velocity giants, 
subdwarfs. A few blue stars were measured. 

Stellar settings and settings on galaxies with wel 


regions of galaxies were measdred by offsetting from tk | 
nucleus or-a nearby reference star. Nearly all galaxi) 
measurements were made directly above sky by of 
setting the photometer between the galaxy and | 
preselected clear sky region. A typical observing ruil} 
on each galaxy with each diaphragm size employe} 
consisted of five deflections, one in each of the fou} 
colors in order of increasing wavelength followed by | 
repetition of the first color as a check. For faint galaxie 
multiple deflections were taken to increase accuracy} 
Throughout each measurement run an offset guide | 
was used to keep the object centered while bein} 
measured and to permit recentering after offset to sky| 

Two stars were used as primary reference and ex 
tinction stars when using the EMI photocell; thes) 
were the very blue star Feige 95 (Feige 1958) and ; 
nearby reddish star. Eight other stars were used a 
secondary standards. One star was used for a standart 
in the red work. Extinction measurements generally 
were taken at the beginning and end of each night an¢ 


Taste IV. Measuring diaphragms. 


Diameter in minutes of art 


Diameter in mm 60-inch 100-inch ; 
0.809 0.365 0.216 | 
1.328 0.600 0.354 © 
2.264 1.023 0.604 
4.032 1.822 1.075 
7.028 3.177 1.874 | 


ntervals of one to three hours throughout the night. 
ntifications, magnitudes, colors, and observational 

bable errors for the reference stars are given in 
ie VII. 

)bservations were reduced to outside the atmosphere 
4 fairly standard techniques. For both the four blue- 
+ and four red-cell colors, reduction was carried out 
‘terms of one magnitude and three-color indices. Ex- 
§ ction coefficients were allowed to vary from night 
‘night, and all the available observations were con- 
+ ered together to derive the best set of extinction coef- 
fents. A careful investigation of the dependence of 
‘ inction on color showed significant effects only for 
#: color index 2-3 and the magnitude 3 


Calibration of the Color System 


Once atmospheric extinction has been removed from 
} otoelectric measurements, instrumental magnitudes 
o be represented by 


m= —2.5 log [ £(0)G(2)ae-+ const. (2) 


| ; 
Is the absolute energy curve of the object observed, 
id G has been defined in Eq. (1). If we expand E(x) 
Eq. (2) in a Taylor series about an arbitrary x0, we 
in reduce the result to the form 


=—25 log| E(x) fear 
dE 
0 fatino 
ae 


Tye 


pal) 


fo- x0)’G(x)dx+.. +c. (3) 


low define the set of moments y: 


W1= X03 nf Gear f wear, (4) 


4 f Glidden i. PG aids N=2,3,4,.... (6) 


nserting the moments in Eq. (3) the first derivative 
2rm vanishes, so we have 


Ck 


1 
1= —2.5 lo {BOn)+5—— vt..| 
2 ; 2 dx? 71 3 


—2.5 log 4p G(x)dx+C. (6) 


‘o the extent that terms in y2 and higher derivatives 
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TaBLE V. Calibration constants of the photometric system. 


Band A v1 12 - Color H 

1 2.683 2.664 0.006 1-3 —0.343 

2 3.371 2.392 0.023 2-3 0.345 

3 3.026 2.068 0.016 ; 

4 0.141 1.682 0.009 3-4 2.885 
: 3-6 4.830 

5 1.696 1.554 0.005 5-6 —0.922 

6 2.618 1.262 0.005 

7 2.280 1.136 0.003 6-7 0.338 

8 1.001 0.002 6-8 0.348 


2.270 


may be neglected we have 


m= —2.5 logE(y1)—2.5 log f GGa)dx+C. (7) 


Thus, given the integrals of the G functions, we can 
Bctercaine the shape of the energy curve. It is of course 
obvious that observations with response functions G 
cannot show detail in the energy distribution over 


- intervals of « significantly less than the bandwidth of 


the G functions. As the monochromatic energy distri- 
bution £ is smoothed over larger and larger intervals 
Ax, until Ax approaches the bandwidth of the G 
functions, values of E(y1) computed from Eq. (7) will 
be better and better approximations to the energies 
read from the smoothed energy distribution at the 71 
points. This is the case since smoothing the energy dis- 
tribution reduces the terms neglected in Eq. (7). 
Except possibly in the ultraviolet and violet, most 
observed energy distributions smoothed over Aw inter- 
vals comparable to the G function bandwidths have 
fairly small second and higher derivatives. Thus, rela- 
tive energies, computed from Eq. (7), generally give a 
good determination of the shape of the true smoothed 
energy distribution. Errors in the determination of the 
G functions generally set the limit of accuracy attain- 
able. Further discussion of the validity of the effective 
wavelength approximation is given later in this section. 
Transposing Eq. (7) and writing 


mp=—2.5logE(yi1), A=2.5 log f Gar, (8) 


_ we have 


ME=MorstA+C. (9) 
If we now define 
then color indices may be placed on an absolute energy 
basis by taking 
Cr (a,b) =Covs (a,b) + A (a,b). (11) 
Values of A, H, y1, and y2 for the photometric system 


are given in Table V. The absolute values of A in 
Table V have no particular significance since no attempt 


’ was made to determine C in Eq. (9). The A values do 


show, however, the relative “‘constant energy” areas 
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Taste VI. Computed and observed colors for the sun and 10 Lacertae. 


Color BT) Obs. A 
1-3 —0.255 -—0.254 —0.001 
23 —0.156 —0.156 0.000 
34 —0.261 —0.23 —0.058 
3-6 —0.@4 a ~as 
34 —0.318 —0.314 —0.004 
6-7 —0.174 —0.233 0.059 
6-8 —0.388  —0.589 0.201 


A=Compated minus observed 


of the bands withm the systems for each photocell. 
Since the constant C is different for the EMI and 
Farmsworth photocells, Eq: (10) does not apply for 
the determination oi H(3,6) in Table V. H(3,6) was 
found by using the Minnaert (1953) solar energy curve 
to calculate the 3-6 color index expected for 1 the sun 
according to the general color-index equation 


(A—B)=—2.5 el [EGC / f Ga(2)az] 
+25 oe f E(x)Ga(x)dx if f Gale)ds|. 


(12) 


Comparing the calculated solar color with the observed 
instrumental 3-6 color index of a G2V star, interpolated 
from the available observations, yielded the value of H 
given in Table V. 

In Table I11 it will be noted that the response function 
for band 5 shows a secondary maximum to the infrared 
side of the primary maximum. Corrections to remove 
the effect of the secondary maximum have been applied 
to all the 3-6 colors. This has been done by calculating 
the correction for various energy distributions and 
expressing the correction as a ‘function of the observed 
3—i color index. The correction was related to 5~7 since 
the secondary maximum in response band 5 corresponds 
closest to band 7. The corrections in no cases exceeded 
0.12 mag. and are probably quite well determined. In 
calculating the data for band 5 in Table V the response 
function was cut off at r= 1.225, the correction having 
accounted for the remainder of the function. 

After application of the special corrections to the 5-6 
color indices, all color Imdex measurements on the 
natural imstrumental system were reduced to an 
absolute energy system by use of Eq. (11) and Table V. 
It is im this form that all color indices are given in 
Table VIII to XII Band 3 instrumental magnitudes 
have been reduced to approximate visual magnitudes 
by addition of a constant determined by comparing 
stellar instrumental band 3 magnitudes with published 
visual photoelectric magnitudes. Because of the limited 
number and type of stellar observations obtained, no 
attempt has been made to relate the magnitudes or 
colors more closely to any conventional system. De 
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-discussed linear relationships of the color system to { 


-Lacertae. Using the G(x) functions absolute energ 


un a 

Comp. Obs. G2 A AS ; 
1.14 1.148 0.006 -0.00 . 
0.506 0.492 0.014 0.01 4 
0.302 0.387 —0.085 —0.07 #4 
0.41 0.41 0.00 0.00 : 
0.021 0.038 —0.017 —0.01 : 
—0.027  —0.082 0.055 0.06 : r 
—0.013 —0.204 0.191 0.20 i 
i 

3 

3 


Vaucouleurs (1961) using the observations contail 


in this paper, as originally tabulated in thesis form, } 


Johnson-Morgan U, B, V system. In using such f 
tionships, however, one should bear in mind that 
systems differ widely in response functions tra asfi 
mations are generally nonlinear, especially in the shon# 
wavelength regions. Furthermore, the transformatieg 
will differ according to the type of object considere 
so transformations based upon stars will not neces 
apply strictly to galaxies. If it is desired to t 
the observations in this paper to any other sys 
while retaining the maximum precision, it Is rec 
mended that direct integrations be performed over 1 
desired response bands. Equation (12) can be used ff 
this purpose with Ga and Gz referring to the d esire 
system and E(x) taken from a smooth plot oi each & 
dividual object, or group of very similar objects, tak ; 
from the data in this paper. 

As a final item relating to the color system calibrati : 
we consider deviations of the calibration from a tru 
absolute energy system. Two energy distributiol 
available for comparison with the observations af 
those of the sun and a blackbody approximation to i 


color indices were computed, for the Minnaert (19 
solar energy distribution and a 20000° blackbo 
according to Eq. (12). Table VI shows the compari 
between computed and observed color indices. 

color indices of 10 Lacertae are direct observations w 
the color indices compared with the solar values 2 
interpolated among available observations to represent! 
a G2V star. ' 

The agreement of the residuals from 10 ieee an 
the sun, according to Table VI, is quite good. Relativel 
small calibration errors are indicated for all color indice 
except 6-8, with only 34, 6-7, and 6-8 showing si; 
nificant residuals. Since 10 Lacertae has no observe 
3-6, and the solar color index has been matched 
calibrate 3-6 initially, this material gives no check 0} 
the calibration of 3-6. 

Additional comparisons of the color system 
outside information are afforded by comparison with 
absolute energy calibrated flux scans according to Code 
(1958 personal communication). Magnitudes were read 
off the scans at the effective 1/A points, and cc 
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TaBLE IX. Stellar colors determined with the blue photocell. 


Star® Spectrum 1S) es, 3-4 3 

u Cas G5Vp 1.16 0.56 0.45 5.18 
53 Cas B8Ib 0.51 0.20 0.27 5.48 
HD 81192 G7III+-hv 1.69 OG 0.62 6.66 
6 UMa FoI 0.91 0.36 0.30 eis) 
a Leo _Maill 3.26 1.35 1.24 pedi2 
¢ Leo FOIL ORS ORs 0.18 S200 
+51°1696 GOsd 0.86 0.43 0.40 9.87 
HD 103095 G8sd 1.26 0.59 0.50 6.52 
HD 108910 K3-+III hv 2.93 12 0.98 tees 
31 Com GOIII 1.19 0.54 *0.44 4.99 
B Com GOV 1.06 0.45 0.36 4.29 
HD 122563 GOIIL 1.34 0.64 0.64 6.40 
HD 124752 KOV 1.63 0.71 0.56 8.61 
a Boo K2I1Ip 2.42 1.05 0.82 0.25 
6 Boo F7V 0.94 0.40 0.31 4.02 
p Boo K31 PAS 1.11 0.85 3.89 
a Boo F2V 0.73 0.28 0.23 4.35 
HD 128902 K2I1I1 hv 2.96 1.25 1.01 Onto 
+26°2606 F4sd 0.64 0.30 0.31 9.68 
HD 134063 GSIlI-++hy 1.65 0.78 0.60 7.95 
HD 134439 KOsd 1.26 0.60 0.52 9.13 
5 Ser A F8III-IV hv 1.03 0.44 0.32 5.00 
HD 140283 F5sd 0.69 0.32 O53) eeu 
d Ser GOV 1.05 0.43 0.46 4.50 
6 Oph MII Se do 1.08 3.20 
y Her A9OTIT 0.89 0.25 0.13 3.60 
n Drab G8III 1.74 0.75 0.58 2.90 
HD 148349 M2-+II1: hv 3.42 1.43 Al. 25 5.69 
n Her G4III 1.70 0.75 0.57 3.63 
48 Her KiII-III hy 2.38 1.04 0.81 6.84 

y Dra KSIII 3.20 ibe Sih 1.03 2.68 
HD 168322 G8IIT hv 1.82 0.79 0.63 6.29 
110 Her F6V 0.89 0.35 0.28 4.17 
x Cyg KOUI 1.87 0.79 0.60 3,98 
o Dra KOV 1.64 0.67 0.50 4.78 
16 Cyg A G2.5V 1,22 0.52 0.40 6.00 
16 Cyg B G4V 124 102545 0542 6.26 
HD 186776 M3III hv 3.09 1.28 1.04 6.81 
55 Cyg B3la 0.41 0.22 0.30 4.78 
HD 199191 G8III hv 1.74 0.78 0.61 Gaol 
61 Cyg A K5V 2.30 0.96 0.87 5.49 
61 Cyg B K7V 2.48 1.04 1.03 6.44 
HD 201891 F9sd 0.81 0.36 0.33 135 
10 Lac O9V —0.25 —0.16 —0.20 4.59 


8@ UMa, HD 103095, 31 Com, HD 124752, BD +26°2606, HD 134063 
and HD 134439 were measured twice. HD 128902 and 10 Lac were measured 
three times. p Boo and 110 Her were measured five times. All other stars 
were measured only once. 

b Measures include light of a faint companion. 


agreement was found to be excellent in bands 1 and 3 
while in band 2 the integrated magnitude was a few 
hundredths of a magnitude brighter than the smoothed 
scan predicted. For the Minnaert (1953) solar energy 
curve, comparison of integrated mean magnitudes in all 
response bands with magnitudes read directly at the y1 
points showed excellent agreement for band 1 and all 
the red and infrared points. The integrated magnitudes 


for bands 2, 3, and 4 fell slightly below local fluctuations 


in the solar energy distribution. However, the smooth 
curve drawn through all the integrated magnitudes at 


WET AS NI GS AUN GD pia eg bi: 


The general conclusion would seem to be that the coll} 
.and magnitudes given in this paper are quite rep} 


the y1 points averages out local fluctuations very w 


sentative of the actual energies, locally smoothed, \| 
the 71 points. It also seems probable that the 
system, including the mean corrections in Table \j 
does not deviate very appreciably from an absoli|) 

energy system. Hl | 


OBSERVATIONS 


Stellar Observations 


observations on each object, for all ‘tae used for "| ] 
ence and extinction. Stars are listed in order of rig} 
ascension. oe hago common designations are ine|) 


the table. 
Table [X contains colors and band 3 Byer | fi 
all the stars observed with the EMI photocell —_ 


Keller en These stars tend to show cect f 
characteristics. HD 122563, although not listed bi 
Keenan and Keller, is an extreme example, color}, 
metrically, of this class. Giants not marked “hv” 


“pb” are more or less normal population I objects. Tabll 


3-6 5-6 6-7 


; t 
Star Spectrum 6-8 | 
HD 81192. G7III-+hv 0.92 0.24 —0.00 —0. 
a Leo M2III 2.17. 0.87, Ose 
¢ Leo - “FOIL” 0.05 —0.08 -—0.12 —0. 
HD 103095 G8sd 0.63 0.18 —0.03 —0.1 
HD 108910 K3-+III hv 1.67 0.48 0.07) O18 
31 Com GOIII 0.45 0.05 —0.06 —0.16} 
6 Com GOV 0.35 0.02 —0.09 —0.19 
HD 122563  GOlIII 0.97 0.24 0.00 —0.03 
6 Boo F7V 0.15 —0.02 —0.10 —0.25 
p Boo K3III 1.23 0.28) OL020 cous 
« Boo F2V —0.01 —0.05 —0.11 —0.30! 
HD 128902 = K2III hv 1.54 0.46 ‘0.08 0115} 
+26°2606 F4sd 0.15 —0.00 —0.07 —0.21. 
HD 134063 GSIII-+hv 0.77 0.20 —0.01 —0.04; 
HD 134439  KOsd 0.62 0.13 —0.05 —0.08} 
d Ser GOV 0.55 0.01 —0.10 —0.26 
5 Oph MIIII 2.05  .0.71 7 0.09; nes 
y Her AOIIT —0.13 —0.12 —0.14 —0.32 
10 Lac O9V —0.31 —0.23 =0.59) 


a HD 81192 was measured 20 times and served as standard reference 
star for the group. 8 Com was measured twice and ¢ Leo, HD 108910, 31 
Com, and o Boo were measured twice in 6-7 and 6-8. All other values are 
single determinations, 


Object 


TABLE XI. Extragalactic colors determined with the blue photocell. 


Type 


(Crd IN OC WIR TEL (OAT COIN 1B, AD TRL 
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Field 


NGC 205 
NGC 221 


NGC 224 


NGC 598 
NGC 2903 


Haro: 22 
NGC 3031 


NGC 3115 


NGC 3351 


NGC 3368 


NGC 3510 


NGC 3810 


NGC 3998 


NGC 4216 


NGC 4245 


NGC 4274 


NGC 4278 


NGC 4283 


NGC 4374 


Ep 
Ep 


Sb 


Sc 
Sc 


Sb 


E7 


SBb 


Sa 


SBc 


Sc 


SO 


Sb 


SBa 


Sa 


Ei 


EO 
SO 


ion 
oO 


1.02 
0.36 


1.02 
36 


0 

3.18 
1.82 
1.02 
0.60 
0.36 


1.02 


Or 
Ww 
nN 


[oe} 
is) 


ROR GOFF 
w& 
rox 


Pw Corr 
nN 
nN 


CRF w& 
ioe) 
iS) 


i 
i=) 
iS) 


ec 
8 


1-3 2-3 3-4 3 
1.14 .50 54 11.34 
1.07 -46 .40 13.04 
tol! 74 . 66 9.28 
1.61 14 .69 10.12 
1.80 .80 518) 6.60 
1.83 .83 .70 7.31 
1.85 .83 73 8.20 
1.84 .83 AS) 9.10 
1.87 85 74 9.93 
0.70 a) -45 Te 52 
1.09 sey? 54 9.65 
1ysalal OU! sel 10.30 
1.06 you .60 10.99 
1.01 .48 .56 11.60 
0.88 -42 aoe) 12.16 
0.55 -Sye .10 14.84 
1.76 .81 74 9.49 
1.87 .85 shel 10.62 
1.69 81 64 9.80 
1.68 .82 - 66 10.14 
1.79 .84 .67 10.95 
We 2 7, .60 .63 10233 
133: . 66 . 66 10.86 
1.30 . 66 .64 11.34 
isis) .61 sO LS ee 
1.44 . 68 . 66 9.88 
1.50 AYP) .68 10.25 
1.59 .76 .69 LOBEL 
1.64 78 Bho 11.34 
0.68 34 eS 14.04 
0.67 BOS) .36 14.65 
0.57 BOO) alle 14,94 
VOn87 -42 -49 10.86 
0.98 -48 £93 11.80 
122 .59 . 64 13.18 
1.65 .80 .65 10.98 
1.67 .80 .66 11.26 
1.69 .80 .68 11.58 
17S .84 .67 11.90 
1.78 82 14 11.26 
1.92 .88 .80 12.06 
1.34 60 .59 12.62 
1.47: Ue tee 1176 
1255 AN .65 12501 
1255 74 .66 1257 
1.49 aril .68_ 13.53 
1.56 ofl S7fi\ 10.85 
1.64 silif Bi) 11.72 
1.62 .78 .64 10.61 
1.60 .76 .64 10.81 
iene .78 . 66 1 GPa tab 
1.67 81 65 11.45 
1.70 .80 -62 P2255 
1.72 .81 .69 10.20 
1.76 .83 Wil 11.01 


Remarks 


M32 


M31 


M33, a 


M81 


M95, Leo group 


M96, Leo group 


Haro 26 
8"0 W, 32°03NV 


Leo group 


UMa cloud 


Virgo cluster 


542 E, 136” N 


13"6 S,b 
Centered 
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TABLE XI (continued) 


Object Type Field 1-3 2-3 3-4 3 Remarks 
NGC 4406 E4 1/82 1.68 ave 64 10.30 M86, Vir cl. 
0.60 1.74 eu ae 66 11.30 : 
NGC 4449 I 1.82 0.51 24 ars) 10.28 c 
0.36 0.46 20 -23 12.33 d 
Anon EO 0.60 1.57 84 78 13.58 e 
NGC 4486 EO 3.18 1.68 78 .69 9.62 M87, Vir cl. 
1.82 eval 81 68 9.97 
1.02 1.72 .82 68 10.45 
0.60 1.69 .80 wi 11.02 
0.36 1.71 .82 i172 11.66 
NGC 4501 Sc 3.18 137 64 65 10.09 M88 Vir cl. 
1.82 1.44 69 -69 10.58 
1.02 1.59 75 74 11.25 
0.60 173 81 78 11.93 
0.36 1.82 84 .79 IZ ot 
NGC 4548 SBb 3.18 1.39 68 Bik! 10.74 Virgo cluster 
1.82 nip | -69 61 11.23 
1.02 1.61 ays 3 11.80 
0.36 1.67 Bie .70 12.89 
NGC 4565 Sb 1.02 NSCS) 83 Atif 11.38 
0.36 1.95 92 88 12.78 
0.60 1.39 -67 66 14.31 6972 W, 82°7 N 
1.02 1.09 54 58 13.82 174” W, 193” N 
NGC 4579 Sb 3.18 1.44 .69 62 10.10 Virgo cluster 
1.82 1.48 .70 58 10.55 
1.02 1.58 8: -68 11.04 
0.60 1.60 78 - 68 Ha begoyl 
0.36 1.61 .719 71 11.97 
NGC 4627 0.60 1.05 45 42 14.08 f 
NGC 4631 Sc 1.02 1.03 -46 54 1252 Offset, g 
1.02 0.47 «20 .36 12.24 Offset, g¢ 
1.02 0.74 38 -4l 11.97 Offset, g 
NGC 4643 SBO 1.87 165 76 69 11.04 
0.60 172. 81 68 11.68 
0.22 1.74 83 .70 12.72 ‘ 
NGC 4698 Sa 3.18 1.54 .70 66 10.90 Virgo cluster 
1.82 1.56 74° 66 11.20 
1.02 1.64 78 62 11.60 
0.60° + 1.64 76 66 12.04 
0.36 1.64 76 ahd 12.51 
NGC 4699 Sb 3.18 1.41 67 66 9.88 
1.82 1.48 72 64 10.19 
1.02 Beattie 55, 75 66 10.38 
0.60 1.60 76 66 10.78 
0.36 1.65 76 67 11.24 
NGC 4713 Sc 1.82 0.60 30 34 11.86 
1.02 0.62 aod mo 12.38 
NGC 4762 Sa 0.60 1.61 75 -58 11.82 Virgo cluster 
1.02 Loy 63 -66 14.16 65"1 E, 108” V 
NGC 4865 E6 0.60 1.64 84 64 14.05 Coma cluster 
0.22 1.64 77 68 14.53 
NGC 4874 SO 1.08 1.70 84 68 13.00 Coma cluster 
0.60 1.69 85 68 13.49 
0.35 eS .89 69 13.98 
NGC 4881 El 0.60 153 .86 60 14.17 Coma cluster 
0.22 ies 86 65 14.84 
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TABLE XI (continued) 


SL) Nao aad 


Object Type Field 1-3 2-3 3-4 3 Remarks 
NGC 4889 E4 0/60 1.76 .89 .66 12.94 Coma cluster 
0.22 1.81 .90 .68 13.90 
NGC 4907 S33) 1.08 1.43 .70 64 13.80 Coma cluster? 
0.60 1.45 Bie .61 14.16 
} 0.35 1.65 719 sift: 14.77 
0.22 iit 85 .69 15.25 
Ic 4051 Ei 1.08 1.68 .84 -63 13.84 Coma cluster 
0.35 1.70 83 .70 14.48 
NGC 4911 Sb 0.60 1.24 58 .60 Sig Coma cluster? 
0.22 iL Stl 80 5S 15.04 
NGC 4921 Sa 487 1.48 .70 .62 12.86 Coma cluster? 
1.08 157 78 .67 13.42 
| 0.60 1.64 82 .68 14.03 
| 0.35 1.70 84 .69 14.49 
NGC 5005 Sb 1.82 1.41 68 .62 10.44 
| NGC 5055 Sb 1.02 1.21 .58 .56 12), i04 54”2 N, M63 
| 1.02 1.14 soy .62 12.69 5472S: 
NGC 5194 Sc 1.82 1.09 le .56 10.10 - M51 
; 1.02 1.13 55 sey 10.61 
0.60 1.26 58 .62 11.18 
0.36 1.36 .62 61 11.85 
1.02 0.99 .46 POD 12.11 461 E, 57"0}S 
| NGC 5195 Ep 1.02 1.69 .80 .80 Abit M51 companion 
| NGC 5248 Sc 3.18 1.02 47 56 10.62 
| 1.82 1.09 52 .56 11.08 
1.02 1.14 57 .60 2 
0.60 1.14 5 64 12.28 
0.60 0.84 .37 46 13.99 57"0 E, 4878 S 
NGC 5363 I 1.02 tends 83 66 11.10 h 
NGC 5457 Sc 1.82 1.14 54 oBve 10.90 81 W, M101 
0.60 et? AGW .56 IDES i Centered 
NGC 5560 (SBc) 1.02 1.14 56 61 13.10 
NGC 5566 SBa 1.82 1.56 74 64 11.22 
0.60 1.66 80 .67 11.90 
0.36 pay 0) 79 tli 1222 
NGC 5576 E4 1.02 1.56 76 63 11.60 
NGC 5633 Sb 1S 02; < 0.94 43 .50 12.67 
; 0.36 1.05 53 OU 13.62 
NGC 5850 SBb 1.08 1.69 81 .66 12.42 
0.35 i733 82 afl 13.32 
NGC 5907 Sb 0.60 1.60 .76 .73 12.91 
1.02 1eSi .67 sol 12.83 NOLO IIe ty ens 
1.02 1.16 .54 .56 13.29 46"1 EF, 111” $ 
1.02 0.86 38 47 13.76 84/1 E, 198” S$ 
NGC 6239 (Sc) 1.82 0.71 36 sei 12252 
: 1.02 0.72 39 rO2 WHS 
0.36 0.64 39 aah) 13.90 
NGC 6503 | Sc 1.82 1.06 Sil 56 10.86 
0.60 1523 .59 62 1227 


® Centered on brightest patch, slightly SE of photographic center. 
b M84, Virgo cluster. : 
© Centered on the over-all structure, 
b Centered on stellarlike nucleus slightly E of the general center of the galaxy. 
© Virgo cluster dwarf elliptical galaxy, 7’ N, 4!4 W of NGC 4486. 
f Small galaxy just N of NGC 4631, This galaxy resembles an elliptical galaxy except that it has no significant nucleus. 
® Since this galaxy shows no distinct nucleus, offsets have been made from a star near the edge of the galaxy and directly WN of the center. Taken in 
1e order listed in the table the offsets are 54”2 E; 103” W, 48"8 S; 78"6 E, 27%1 S. Although no real nucleus is seen, the measure 54” E of the offset 
ar may include light from the outer part of a central bulge. Both other regions measured include bright visual patches, probably H11 regions. 
h A fairly bright star near the nucleus is included in the measures. 
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TABLE XII. Extragalactic colors determined with the red photocell. 
Object Type Field 5-6 6-7 6-8 3-6 6 Remarks 

NGC 224 Sb 1/08 +53 .09 val 1.36 6.76 M31 
0.60 .54 11 24 1.34 Gels) 
0.35 54 12: B25 fie32 8.67 

NGC 2903 Sc 1.08 .36 .10 .19 0.98 9.95 
0.60 .38 Bea .19 0.99 10.61 
(OST 9) 39 .09 oa, 1.00 11220 

NGC 3115 E7 1.08 -40 .10 a20 sil 8.79 
0.60 -40 .08 .20 oo OMlZ. 
0.35 42 .08 19 i io 9.46 

NGC 3351 SBb 1.08 40 .08 ppl SII) 10.17 M95, Leo group 
0.60 41 .10 Pll) 1.18 10.59 
0.35 35 .07 .18 cae 11.07 

NGC 3810 Sc 1.08 27 03 ike) 0.89 10.84 Leo group 
0.60 «29 .08 «20 tee 11.43 
0.35 sel .10 12 1.14 12.09 

NGC 4486 EO 1.08 -40 10 19 ie: 9.28 M87, Vir cl. 
0.60 -40 10 vip) ds 9.89 
0.35 .38 11 20 1.10 10.60 

NGC 4631 sie 1.08 58) 16 30 0.80 NAA a 

NGC 5248 Se 1.08 530 10 20 1.01 10.66 | 
0.60 33 09 aly 1.10 fay ‘y 
0.35 532 11 .20 He 11.70 4 


®78"7 E, 27"2 S from reference star. See Table XI. 


Extragalactic Observations 


Table XI gives the color indices and band 3 mag- 
nitudes for all galaxies observed with the EMI photo- 
cell. The NGC number or other identification is given 
for each system along with the Hubble structural type 
taken from the Humason, Mayall, and Sandage (1956) 
redshift catalogue, wherever possible. One type, given 
in parentheses, was determined by the author. The 
galaxy Haro 22 was taken from a list of blue galaxies 
given by Haro (1956). The aperture employed for each 
measurement is given in minutes of arc, and unless 
otherwise indicated the diaphragm was centered on 
the nucleus. The last column of Table XI contains 
remarks, including references to notes following the 
table. 

Table XII contains all extragalactic observations 
obtained in the red. The material is arranged as in 
Table XI with the addition of the color index 3-6 which 
ties together the blue and red observations. In order to 
formulate 3-6 it was necessary to interpolate among 
band 3 magnitudes to derive the value of 3 corre- 
sponding to the field size used to measure 6. This 
problem arose because of the different scales of the 
60-inch and 100-inch telescopes, and the fact that all 
galaxies measured in the red at the 100-inch had been 
measured in the blue at the 60-inch. However, for the 
diaphragms employed in measurement, a given field 
size at the 60-inch telescope corresponds closely to the 
field size of the next larger diaphragm at the 100-inch 
telescope. Thus the interpolation corrections were small 


i 
and it is unlikely that serious errors are present in ‘ 
due to the interpolation. Stellar colors in Table X arej 
of course independent of diaphragm size. ue 


J 
4 

Observational Errors : | 
: 


Observational errors may be divided into two cate-| 
gories, internal random errors of measurement on each 
individual galaxy or star, and tie-in systematic errors 
between different galaxies or galaxies and stars of} 
between different nights. Some information is available! 
for estimating the random errors of measurement on; 
individual galaxies by comparing the first and last. 
measurement made with each diaphragm size. The fir, 
and last measures were almost always taken with the 
ultraviolet filter when using the EMI photocell; thu 
at least two determinations of 1-3 are available for each 
diaphragm size on nearly every galaxy. An investigatio 
of the 1-3 residuals about the final mean value shows, 
as expected, a very good correlation with the band 3 
magnitude, due allowance being made for the different 
light-gathering powers of the 60-inch and 100-inch | 
telescopes. Brighter than 3=12.5 at the 60-inch or 13.7 
at the 100-inch, 1-3 measurement errors are consistent 
less than 0.02 mag. Extra deflections were taken for the 
faintest systems, so it is doubtful if random errors as 
great as 0.05 mag. exist in 1-3 at any magnitude. 3-4 
random errors should be about the same as in 1-3 since 
deflections in 1 and 4 were quite similar in most cases. 
Because of the increased sensitivity of bands 2 and 3, 
random errors in 2-3 probably do not exceed 0.02 or 


3 mag. at any magnitude. Relatively little informa- 
n is available for estimating errors in the red obser- 
tions; however, what repeat observations are avail- 
le rarely show deviations of individual measures from 
2 mean larger than 0.02 mag. for any of the red 
otocell color indices. 

Tie-in systematic errors between galaxies and stars 
id between different nights are much more difficult to 

‘termine than random errors. Except for a few cases 
| the red, no repeat measurements were made on 
' laxies using the same field diaphragm. Reference stars 
vere referred to frequently but provide no check on 
‘sidual errors after reduction since the standards 
ovided the basis for extinction determinations. For 
ve or six galaxies some observations using different 
aphragm sizes were obtained on different nights. On 
te basis of these observations, assuming smooth color 
vadients with radius, night-to-night systematic errors 
irely exceed 0.03 to 0.04 mag. in all colors. It is 
oubtful if any exceed 0.1 mag. 

Systematic errors for magnitudes in band 3 or band 
are larger than systematic errors in the color indices. 
‘his is the case primarily because magnitude extinc- 
‘on is much more variable than color extinction and 
arder to determine accurately. Differential magnitudes 
etween various diaphragm sizes are probably as reli- 
ble as the colors, but the actual values of individual 
nagnitudes are definitely less accurate. In view of the 
ywer accuracy of magnitude determinations, the color 
adex 3-6 relating the measures with the two photocells 
s the least reliable color. Correlations of 3-6 with 
ther colors imply that few if any discrepancies as 
arge as 0.1 are present in the 3-6 colors. 

Two factors preclude any precise comparison with 
nost other extragalactic observations. The uncertainties 
£ color system transformations has already been men- 
ioned and will enter into any comparison of dissimilar 
ystems. A second consideration arises from aperture 
ffects. Since color varies with radius in most types of 
alaxies, and different observers have used widely 
aried apertures, a direct comparison between observers 
s difficult. Errors arising from the above sources are 
nore or less comparable with normal photoelectric 
recision. De Vaucouleurs (1961) has shown that the 
tellar observations contained in this paper are quite 
onsistent with standard U, B, V data. Furthermore, 
here are no obvious inconsistencies in the extragalactic 
bservations in comparison with de Vaucouleurs. A 
nore detailed comparison would not seem to be war- 
anted at this time. 


DISCUSSION 


Comparison of Stellar and Extragalactic 
Color Indices 


As an introduction to the analysis of the extragalactic 
olors given in this paper we consider first a comparison 
yetween stellar and extragalactic color indices. Vari- 
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TABLE XIII. Comparison of high- and low-velocity star color 


indices. 

Color G8IIIThv  GS8III G&8sd K2.5+]iI. hv K3III 
1-3 1.78 1.74 1.26 2.94 2.67 
2-3 0.78 0.75 0.59 1.23 aaa 
3-4 0.62 0.58 0.50 1.00 0.85 
3-6 1.60 123 
5-6 0.47 0.28 
6-7 0.08 0.02 
6-8 0.10 0.04 


G8III hv mean of HD 168322 and HD 199191 
G8III Eta Draconis 

G8sd HD 103095 

K2.5+III hv mean of HD 108910 and HD 128902 
K3III Rho Bootis 


ations of color index with such parameters as radius, 
inclination, velocity, and galactic latitude are neglected 
here and are considered in detail in the second paper in 
this series. Figures 2 through 6 present five color-color 
diagrams of the stellar and extragalactic observations. 
In the first three figures, 1-3 is plotted as ordinate 
against the progressively longer wavelength color 
indices 2-3, 3-4, and 3-6. In the last two figures 5-6 
is used as ordinate with 6-7 or 6-8 as abscissa. In each 
figure stellar observations are shown to the left. Lines 
have been drawn to represent particular types of stars, 
viz., main-sequence stars, normal giants, high-velocity 
giants, and subdwarfs, whenever separation into groups 
is suggested. The red observations are insufficient to 
make significant distinctions. The lines representing 
the different types of stars are repeated in the right 
portion of each figure, where the extragalactic obser- 
vations are shown. Different symbols have been used 
for the various basic extragalactic structural types. 

In Figs. 2 and 3 in particular, the separation of stars 
into distinct groups is quite clear. Particularly inter- 
esting is the separation of the high-velocity giants from 
the normal giants in the same manner as the subdwarfs 
deviate from the main sequence. The high-velocity-giant - 
point deviating most extremely in Figs. 3 and 4 is 
HD 122563, which shows extreme population IT spectral 
characteristics. Another comparison between the high- 
velocity giants and the normal giants is shown in 
Table XIII, where color indicates at equal spectral 
type are compared. The Table XIII examples, as well 
as others that may be taken from the stellar color index 
measurements, indicate that the high-velocity giants 
are consistently redder than the normal giants at the 
same spectral type. If equality of spectral type is to 
mean equality in effective temperature, then similarity 
of the far-red energy distributions, where line absorption 
has little effect on the energy distributions, should 
correspond to equal spectral class. Metal-poor stars, 
such as are generally considered to occur among the 
high-velocity stars, should show, if anything, less line 
blanketing than the low-velocity stars, and hence if 
they deviate at all in color index from low-velocity 
stars they should show bluer color indices at short 
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Fics. 2-4. Color-color relationships for stars (left) and galaxies (right). The lines are discussed in the text. 
Coma cluster galaxies are not plotted, and in the most crowded portions of the figures some points have been 
omitted for clarity. The most extreme red points are nuclei of spirals seen edgewise, hence heavily reddened. 


wavelengths. This effect is well known in the subdwarfs, 
and since the same apparently holds true for the high- 
velocity giants it would appear that the high-velocity 
giants have been classified too early in spectral type. 
In each of the color-color relationships the extra- 
galactic colors are seen to lie in a fairly distinct sequence 
displaced from, but nearly parallel to, the stellar 
pattern. The displacement is readily explained as a 
natural consequence of the composite nature of the 
galaxies. Consider, for example, that a certain type of 
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GALAXIES 


galaxy has an energy distribution resembling a G4II 
star over the spectral region from band 1 toJband 3 
i.e., 1-3 for the galaxy is equal to 1-3 for a G4ITI star. 
Now as we examine the galaxy’s energy distribution 
further to the red, energy contributions from stars of 
later and later spectral types become increasingly im 
portant. Thus the galaxy will yield a color index 3-4, 
which resembles a star of spectral type later than 
G4III, say K1III. Passing still further to the red, we 
would find that 3-6 for the galaxy resembles 3-6 of a 
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Fics. 5-6. Color-color relationships for stars (left) and galaxies (right) in the infrared. 


\ 
<3III star. Thus when we plot one color index against 
nother for the galaxy, the longer wavelength color 
adex is found to be too red to correspond to a star 
vhich matches the galaxy in the shorter wavelength 
olor index. This is just what is observed and illustrated 
n Figs. 2 through 6. Note that as progressively longer 
olor base lines are employed the difference in the 
pectral type of the stars matching the two color 
ndices increases, hence the deviation of the galaxy 
olors from the colors of any one type of star grows. 
The reddest color one would expect to observe in a 
falaxy, neglecting reddening by interstellar or inter- 
ralactic matter, should be representative of the latest 
pectral type giants. Figure 6 shows that the 6-8 
olors begin to resemble M-giant star colors, clearly 
ndicating the considerable importance of M giants, or 
ven cooler stars, as contributors to the integrated 
mergent energy from extragalactic systems. For the 
jluest galaxies it is interesting to compare the measure- 
nents, particularly in Fig. 3, with the blue galaxy 
measurements of Hiltner and Iriarte (1958). These 
ralaxies appear to form the blue end of the extragalactic 


sequence, and presumably represent galaxies containing 
a large number of young blue stars, but still enough 
late-type stars to produce the composite effect. There 
is no indication that the Hiltner and Iriarte observations 
do not fit smoothly on the sequence indicated by the 
data in this paper. Further consideration of the syn- 
thesis of extragalactic colors is given in the second paper 
of this series. 

In addition to demonstrating the composite structure 
of galaxies, Figs. 2 through 6 show the well-known effect 
that structural type is correlated with color. The Sc 
and SBc spirals are clearly bluer than the Sa or Sb 
systems, and the elliptical galaxies form a fairly 
homogeneous group among the redder systems. 
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The linear equations are presented which describe the observed color indices of an eclipsing binary during 


the eclipse as functions of the difference of limb darkening coefficients, difference of intrinsic colors, and 
corrections to the rectification constants. The ratio of radii of the components of the eclipsing binary, the 
ratio of their luminosities in one color, and the coefficients of limb darkening in this color are the only 
elements of the system upon which the coefficients of these equations depend. A least-squares solution of 
the equations describing the color indices of YZ Cassiopeiae observed by Kron gives for the star eclipsed 
in primary minimum a difference of limb darkening coefficients 4° —15™=0,35+-0.02 (m.e.) which is 
about twice as large as the value 0.16--0.08 resulting from independent solutions in each color. 


I, INTRODUCTION 


HE problem of determining the elements of an 
eclipsing binary from two-color photometric ob- 
‘-rvations can be, in general, reduced to a least-squares 
lution of the system of equations with corrections to 
ie elements on one side and the observed brightnesses 
their differences on the other side. Part of the equa- 
fons describes the changing brightness in one color 
nly. The remaining equations can be written in one of 
wo forms; they may contain either the brightnesses in 
nother color or the color indices. 
The second form of these equations seems better 
cause the color index is as a rule measured more 
ccurately than the brightness. This is especially true 
or the differential photometers containing several 
aultiplier tubes which are becoming more and more 
ommon. The color indices measured by such photom- 
ters are to a large extent free of the disturbing influence 
if atmospheric scintillation. 

The equations describing the color indices are almost 
ndependent of equations describing the brightnesses 
n one color. Each of these systems of equations contains 
. Separate set of unknowns: In each of these systems the 
erm containing the corrections to unknowns deter- 
nined by the other system are negligibly small. Hence 
ach group of equations can be solved separately. 

The most interesting quantity obtained by solution 
f equations describing the color indices of an eclipsing 
inary during the eclipse is the difference of limb 
larkening coefficients in two colors. This difference can 
ye determined with much higher accuracy than can 
ach of the coefficients separately ; the reason is that in 
he equations describing the brightnesses of eclipsing 
inary in one color the coefficients of the correction to 
he radius of eclipsed star are almost exactly propor- 
ional to those of the correction to the coefficient of limb 
larkening (cf. Irwin 1947). Hence the accuracy of a 
imultaneous determination of limb darkening and 
adius of eclipsed star is very low. 

Methods of determining the differential limb darken- 


* Present address: Warsaw University Observatory, Warsaw, 
oland. 


ing from two-color photometric observations have been 
proposed by Irwin (1947), by Kopal (1946), and by 
Russell and Merrill (1952). However, these methods do 
not pretend to be of high accuracy and do not permit 
us to estimate the error of the result. In an attempt at 
a somewhat more accurate approach we are considering 
in this paper all the first-order effects upon which the 
color index of eclipsing binary may depend during the 
eclipse. 


II. CHANGES IN COLOR CAUSED BY DIFFERENCES IN 
INTRINSIC COLORS AND IN LIMB DARKENING 


A. Fractional Loss of Light 


The surface brightness J of a given point on the 
surface of spherical star is with fairly high accuracy 
given by the equation (Kopal 1949) 


J (y)=J(0)(1—u—v+u cosy+0 cos’y), (1) 


where vy is the angle between the direction to observer 
and the normal to the surface of a star at a given point 
while ~ and v are the limb darkening coefficients 
fulfilling the inequality 0<«w<1. 

The fractional loss of light f of the eclipsed star is 
defined as the ratio of integrals: 


f= | fom comda|| fren cosnio] (2) 


where s is the area of the eclipsed part of the stellar 
surface, S is the total area of that hemisphere of the 
eclipsed star which is directed towards us, and do is 
the surface element. Evaluating the integrals in (2) we 
find (cf. Kopal 1949) 


Ff L(l—u—o)f°+3uf?+zef?]; (3) 


\6—2U—S0 
here f/ is the value which f assumes for ~=v=0, f? that 
assumed for w= 1 and v=0, and f® the value assumed for 
u=0 and v=1. These values are proportional to the 
associated alpha functions a,’ tabulated by Kopal 
(1947) : 

fl=ah; fP=za%; f?=2a2?. (4) 
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We shall denote by 6 the apparent distance of the 
projected centers of the two component stars and by 
Ya, 7, the radii of these stars; throughout this paper 
the subscript a will always refer to the smaller star and 
b to the larger star. We denote further, in accordance 
with usual notations, 


k=1ra/T», (5) 
p= (6—15)/Ta, (6) 


and we define the normalized fractional loss of light 
a(k,p) by the equation 


We shall now consider two cases: occultation (the 
eclipse of the smaller star) and transit (the eclipse of 
the larger star). 

Occultation. In this case we have fa(k, —1)=1 and 
hence aa(k,p)=fa(k,p). Introducing the modified limb 
darkening coefficients 


4g 3 : 
VG TES i non ee). aa (8) 
6—2u,—344 6—2uUg—32z. 


and making use of (3) we write the normalized frac- 
tional loss of light for the occultation in the form 


q=al+ (aq? —a”)atg+ (aa2—a”) ya, (9) 
where a”, a, and a® are defined by formula (7), in 
which f is replaced by fY, f?, or f®, respectively. Writing 
(9) for two wavelengths \; and \» and subtracting we 
find 


da! — ag = (ag? — a”) (44! — xq) 


+i(@a?— a) (yaht—y,"*)s (10) 
Transit. Equation (3) gives for this case 
fol (k, —D=P, NE Degen) 
fo? (k, —1)=3R'®(k), (12) 
fo®(k, —1) = (4—3k)R3, (13) 


where the function ®(%) is tabulated by Tsesevich 
(1940, Table 5). Introducing modified limb darkening 
coefficients 
UyP (R) 
ee ? (14) 
1-—%4— WtuUpP(k) +4 (4— 3k) Rvp 
3(4—3k) kop 
ve ? (15) 
1—m— VotUpP (Rk) +2 (4— 3k) kv, 
we have, by (3), (7), and (11)-(13), 


ap=aU+ (ap? —a”) 25+ (as°—a") yp, 


(16) 


and finally, by subtracting Eqs. (16) written for wave- 
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lengths A; and Az, we have 


ap — a= (ay? — a”) (ary — arp”) 


B. Relation between a*!—a and Color Index of a 
Eclipsing Binary 

Denoting by / the rectified fractional light of tl 
eclipsing system expressed in terms of the combine 
luminosities of its constituent stars and by Lg, Lo, thi 
fractional luminosities of the smaller and the large 
component, we have for an occultation the fraction 


loss of light a 
fa= (1—1)/La, (18 


f=0-Di cr 


and for a transit 


and / for internal tangency in the case of total or annula 
eclipses (in this case pp>=—1 and ap=1) and for 1 
deepest eclipse in the case of partial eclipses. Taki Di 
now the ratio of the fractional loss of light for a give 
p to that for po we get for an arbitrary type of eclips 
the normalized fractional loss of light 


a=[(1—1)/(1—A) Ja, (26 


and hence 
ay»! ay? ‘ 
o—q’= (1- ae oe DP), (21 
1—\™ : 


q—q2 


Equation (21), after adding and ub the tert 
1 =) 


ay P2(1—d™), is 
<(1—P»). = 


AL 
We shall consider the two types of eclipses i 
Occultation. Since a,=fa, we have, by (18) and (20) 


aoa/(1—Ag) = 1/La, (23 


where aa is the value of a for occultation. Let u 
introduce the parameter Z) depending only upon th 
difference in intrinsic colors of stars constituting th 
eclipsing binary system and defined by the equation 


Ly?/Ly™ 
Zo=1———_——_. 
LL? / L,™ 


ei ee i! 
l= AN 


(24 


Tf the color indices C of the constituent stars are define 
by the equation 


C=—2.5 log(L2/L™), (25 
their difference is related to Zp by the equation 
Ca—Cy=2.5 log(1—Z). (26 


fultiplying (22) by L.™ we have, by (23)-(24), 
P= (ag —a.)L.4—(1—P*)(1—L.™)Zo. (27) 


tilarly, multiplying (10) by Z.™, we get, by using 

), the expression for the difference in light of an 

psing binary in two wavelengths 

-[/M= (ag? —a¥) (aq! — aq?) Lg 

ais (@2—a") Ca— yo) La 

—(1—P*)(1—L.)Z. (28) 

e system of equations (28) written for different 

wses of the eclipse can be used for determining the 

intities x.1—x,., ye!— 2, and Zo by the method 

least squares. : 

ransit. Introducing the function Y(k, —1) defined 

the equation 


b, —1)=fs(k, —1)/# 
1) =folk, -1)/ @ 


= ile aa errr a 


have, by (19)-(20), 


c2oe/ (1—)s) = [RY (k, — 1) 2s 4, (30) 
l hence, by (22), 
~)i= f2y™ (k, ‘as 1) (a3— a2) Ly™ 
—{1-[YN(k, —1) LM/ 
LY*(k, —1)L>"]}(1—-2*). (31) 


ce ys<<1, and since |[1/36(%)]—1|<0.17 for 
0.21, the following approximate expression is valid 
en k is not small: 


(k, me 1 
1), Fan 


= i| (ap— xp”) 


: 1 a a 32 
ae Jo: oa? yay) 2) 


Jefining the parameter Z, depending on the differ- 
e in intrinsic colors, 
LL? / L,™ Zo 
L/L Zo—1 


| multiplying (17) by k?Y*(k, —1)(1—L.™) we have, 
(31)-(32), 


= DM = a (4) — ay?) + de (yo™!— yo) 


Z.= 


, (33) 


| —(1-2)L§Z,, (34) 
sre the coefficients are 
BY (k, = 1) (1 Le) (?—2”) 
+(1—L.MZ,)[1—1/3#(k) ](1—P*), (35) 


-RYM[(k, —1)(1—La™) (a»°—a”) 
+(1-LMZ,)[1—1/(4—3k)k](1—™). (36) 
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The approximation is fairly good for k not smaller than 
0.20. The terms containing Z; in (35)—(36) are in most 
cases negligibly small and successive approximations 
are not necessary for determining Z;. 

Simultaneous solution of Eqs. (28) and (34), ie, a 
simultaneous solution for both primary and secondary 
eclipses, is possible when Z and Z; are replaced by one 
unknown, which can always be done if only we know 
the approximate value of, say, Zo. Denoting this 
approximate value of Z» by Zo and denoting 


Z.=%0/(Zo—1), 


the expansion of Z; in the Taylor series with only the 
first two terms retained is 


ZixZ1—[(Zo—Zo)/ (Zo— 1)? ]. (37) 


After substituting (37) into (34) we transfer the constant 
term (1—)2)L,"Z, to the left side of (37). Similarly we 
substitute Zop=Zo+(Zo—Zo) into (28) and transfer 
the constant term (1—/2)(1—L.™)Z» to the left side 
of this equation. The unknown Z —Zp is now deter- 
mined by the simultaneous solution of Eqs. (28) and 
(34). Tf (1—Ap) Lo > (1—Aa”) (1—L.™), it is better 
to express Zp in terms of Z; and not the reverse; the 
difference between the two procedures, however, is very 
small if the approximate value Zp is sufficiently close 
to the true value. 

The simultaneous solution of Eqs. (28) and (34) for 
primary and secondary eclipses can be made not only 
in the case of a circular orbit but also in the case of 
elliptical orbit of small eccentricity. The geometrical 
elements po, Ya, 7, and 7 (inclination of the orbit) 
obtained in the case of an elliptical orbit from each 
eclipse separately on the assumption of a circular orbit 
are somewhat different, but the ratio k is the same for 
both solutions. And & is the only geometrical element 
upon which Eqs. (28) and (34) depend. 


C. Computation of a?—aY and a2@—aY 


To make a solution of Eqs. (28) and (34) we must 
know the elements k and LZ, from a solution based on 
the observations in one color, say \i. We must also 
assume (or determine from the observations in one 
color) the coefficients of limb darkening %™, uy, vo, 
and 2‘! in this color. Knowing these elements we 
compute for every normal point of the light curve the 
normalized fractional losses of light for wavelength 1, 
using for partial phases of the occultation the equation 


agt= (1—)/L.™; (38) 
for totality we have a,'=1 and for transit 
ay t= (1—)/PY™(k, —1)(41—L™). (39) 


The values a,™ and a, found from (38) and (39) are 
entered into Tsesevich’s (1940) table of the function 
p(k,a) for the assumed value of u™, The values #, 
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k 
Ne 1.00 0.95 0.90 0.85 0.80 0.75 0.70 0.65 0.60 0.55 0.50 0.45 0.40 0.35 0.30 0.25 0.2 
1,00 0 0 0 0 0 0 0 0 0 0 0 (0) 0 (0) 0 0 ( 
0.95 | —34 — 32. — 30 =—28 —25 =23, —20 —18 =16 —13 °—11 — 9 = (8 = 6) 4a 
0.90| — 81 —78 =—7 —69 —63 —57 —51 —46 —40 —34 —29 — 24 = 20 —15 — 1200 eee 
0.85 | —129 -126  —121 —113 105 96 86 77 68 58 —50 —41 — 34 —26 —20 —14 — | 
0.80 | —172 -—173 -—167 —158 -—147 —135 —122 -—109 -—96 — 84 —71 '— 60 — 48 \— 33=- — 29 — 21) ee 
0.75 | —208 -—215  —211 —201 189 174 159 142 -126 -109 —93 —78 —64 —50 —38 —27 =f 
0.70 237 251 250 242 228. —212 —-194 —175 —155 -135 -116 -—97 —79 ~62 —47 —34 —@% 
0.65 | —258 —282 —285 —279 -—265 —248 -—228 -—206 -—183 -160 -—137. -115 —94 —74 —56 —40 —2 
0.60 | —271 -—307 —316 -—312 —299 -—281 —259 -—235 —211 —184 -—158 -—133 ' -109 — 86 —65 —47 = 
0.55 | —275 —325 —341 —341 -—330 -—312 -—289  -—263 -—236 -—207 —178 -—150 —123 —97 —74 —53 —@ 
0.50) —271 —336 -—361 -—366 -—357 -—340  -—316 -—289 -—260 -—229 -197 -—166 -—136 -—108 — 82 —S9 —@ 
0.45 | —260 -—342 -376 -—386 -—381 -—364 -—341 -—313 -—282\ —249 -—215 —182 -149 -118 —90 —65 —@ 
0.40 240 341 386 402 —401 —386 -—364 —335 -—302 -—268 -—232 -196 -—161 -—128 —98 —70 — 
0.35 | —213 —334 -—391 —415 417 405 383 354 —321 —285 —247 —209 -172 -137 -—104 — 75 | —j5 
0.30) —179 —322 -—391 —423 -—430 -—421 -400 -372 -338 -—300 +261 —221 -—182 -—145 —111 — 80 —\5 
0.25 | —139 —304 —387 —427 -—439 -—433 -—414 -—386 —352 -—314 -—273 -—232 -192 -—153 —116 — 84 —'§ 
0.20| — 93 —282 —378 —427 -—446 -—443 -—426 -399  -—365 —326--—784 —242 -—200 —159 -—122 — 88 —5 
0.15| — 43 -—255 -—366 —424 —448  —450 -—435 -—409° -—375 -—336 —294 -—250 ~-—207 -—165 -—126 —91 —6 
0.10 | + 12 225 350 418 448 —454 —442 -—417 -—384 -—344 -—302 -—257 -—213 -170 -130 —94 —6 
0.05| + 70 +191 —331 408 445 455 446 423 -390 -—351 —308  -—263 -—218 -174 -—133 — 96 —6 
0.00 | +131 —-155 —309 396 440 454 448 426 —395 -—356  —313 -—268 —222 -178  -—136 — 98 —6 
—0.05 | +194 -116 -—285 -—381 -—432 -—451 -—447 -—428 -397 -359 -316  -—270  -—225 -—180 -138 —99 —6 
—0.10 | +258 — 76 259 365 422 445 444 427 —398 —360 —318 —272 -—226  -—181  -—139 -—100 — 6 
—0.15 | +322 -— 35 —231 -—346 -—409 -437 -—439 -—424 -—396 -360 —318 —272 —226 -—182 -—139  -—100. — 6 
—0.20| +385 + 6 —203 —326 395 427 432 419 -—393 -—357 -—316 —271 -—226 -—181 -—138 —-100 — 6 
—0.25 | +446 +46 —174 305 380 415 424 412 -—387 —353 —312 -—269 ~224 -179 -137 —99 —6 
—0.30 | +504 + 86 -—146 —284  -—363 —402 -—413 -—403 -—380 -—347 -—307 -—264 ~+220 -177 -—135 —98 —6 
—0.35 | +558 +123 117 262 345 387 400 393 -—371 -—339 -—301  -—259 -—216 —173 -—132 -—96 —6 
—0.40 | +606 +157 — 91 -—240 -326 -—370 -—386 -—380 -—360 -—329  -—292 -—252 -—210 -—168 —129 —93 —6 
—0.45 | +648 +187 — 66 —218 307 353 370 365 —346 —318  —282 -—243 -—202 —-162 -—124 —90 —5 
—0.50 | +680 +212 — 44 -197 -—287 -—334 ~—352 -—349  -—331  -—304 -—270 -—233 -194 -156 -119 — 86 —§5 
—0.55 | +704 +231 — 24 -177 -267 -—314 -—333 —330 —314 -—288 —256 -—221 -—184 -—147 -113 — 81 '—§ 
—0.60| +715 +244 -— 9 -159 —247 -—293 -—311 —310 -—295 -—271 -—240 —-—207 -1i172 -—138 -106 —76 -—5 
—0.65 | +712 +248 + 2 -143 -—227 -271 -—288 —287 —273 —251 -—222 -191 -—159 -127 -—97 —70 —4 
—0.70 | +694 +242 4+ 8 128 206 248 263 262 249 228 202 -174 —-144 -115 -— 88 —63 —4 
—0.75 | +657 +227 + 9 -115 -186  —222 —236 -—234 -—222 ~-—203 -—179 -—154 -—128 -102 —78 —S6 —3 
—0.80 | +597 +199 + 4 -104 -—164 -194 -—205 -—203 -191 -—174 -—153 -131 -—109 -— 8 —66 —47 —8 
—0.85 | +511 +4158 -— 5 —93 —141 -163 -—170 -—166 -—156 -—141 —124 -106 — 87 —69 —53 — 38 —2 
—0.90 | +393 +102 -— 19 — 80. —112 -—126 -—128 -—124 -115 -104 -—91 —77 —63 —50 —38 —27 —1 
—0.95 | +232 +37 — 28 —58 —72 —77 76 72 66 59 51 43 35 28 —21 —-15 —14 
—1.00 (0) 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
q 
0.10 Oo +11 +28 +46 +63 +78 +92 4103 +112 +4118 4121 4121 +118 +4112 4103 +092 +47 
0.20 0 21 52 84 115 144 168 189 204 214 219 218 212 201 184 161 re 
0.30 0 29 72 117 160 200 233 261 281 295 301 299 290 273 248 217 17 
0.40 0 36 89 145 198 244 288 321 346 362 368 365 353 332 301 262 21 
0.50 0 42 104 168 230 285 333 371 399 417 424 420 405 380 344 298 24 
0.60 0 46 115 187 255 316 369 411 442 461 468 463 446 417 377 327 26 
0.70 (0) 50 124 201 274 . 340 396 441 474 494 502 496 478 446 403 348 28 
0.80 0 52 130 211 288 357 416 463 497 518 526 519 500 467 421 363 29 
0.90 0 54 134 217 296 368 428 476 511 532 540 533 513 478 431 372 3¢ 
1.00 O° +54 +4135 +219 +299 +371 +432 +480 +516 +537 +545 +538 +517 4482 +435 +375 +36 


found for every normal point from these tables, are 
entered into the tables of a@?—aY’ and a&—a", which 
should be tabulated as functions of k and p. The tables 
of w?—a” can be easily obtained by taking the differ- 
ences of functions a?(k,p) and a¥(k,p) tabulated by 
Tsesevich (1939, Table III and 1940, Table D). The func- 
tion a? —a/ is tabulated for occultation by Irwin (1947, 
Table 4) and for transit in Table I of the present paper 
(multiplied by k?). Values of a?—a¥ for k=0.96, 0.97, 
0.98, and 0.99 can be obtained from the tables of aU 
published for these arguments by Merrill (1950) and 
tables of a? which can be obtained from the Princeton 
University Observatory on request. 

For annular phases of a transit eclipse, (a@?—a¥)k? 
is tabulated in Table I as a function of arguments k 
and q. The values of g can be found from the table of 
q(k,X) given by Tsesevich (1940, Tables XVI and 


XVIa), where X is computed from the equation 
X= (@—1)/A*(R); (AC 


the values A*(k) are also tabulated by Tsesevich (194( 
Tables XVIII-XXI), who denotes by « our u. Th 
values a?—aY for annular phases can be found als 
directly from the formula 


as? —a¥ =[A (k)/A*(k) \(ao—1), (41 


where A(k) is tabulated by Tsesevich (1940, Tabl 
XVII). 

The terms a®—aY are not yet tabulated. They ca 
be simply expressed in terms of the “boundary inte 
grals” I1,0° and J_1,.° tabulated by Kopal (1947). Fe 
occultation we have 


kp+1 
ae al (4: 
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where the boundary integrals are taken from # tables 
for the argument a defined by the formula 


cosa= 3 (k?p?-+ 2kp—k?+-2)/(kp+1). 


For the partial phases of transit eclipse we have 


(e+ p)k 1 
a [at 2(-+0 )n.2'| 
4—3k k 


1 
eos ————— a’ 
aan ee, 


(43) 


a,2@—al= 


where the boundary integrals are taken from # tables 
for a defined by 


cosa= 3 (kp’-+k+2p)/(kp+1). (45) 


| Finally, for annular phases of transit eclipse we have 


where the boundary integral is taken from a tables for 


a defined by 


cosa= { (1—k)@+-4[ (1+-3k)q—1]}/ 


{—k)¢P—4[(1—3k)g—1]}. (47) 


WI. INFLUENCE OF THE ERRORS OF RECTIFICATION 
A. Analysis of the Observations outside the Eclipses 


The harmonic analysis of the observations outside of 
eclipse is done by solving by the method of least squares 


the two independent systems of linear equations: 


— —1=A™— B™ cosp-+C™ sin’, (48) 
Iy—])' = (A*— A®) — (B*— B®) cos 
+(C’—C™) sin’?@, (49) 


where @ is the phase counted from primary minimum 
‘and Jo is nonrectified fractional light of the eclipsing 
binary expressed in terms of the mean of lights at 
‘phases @=90° and 270° taken as unity. 


The rectified fractional lights /“ and the differences 


Pp, which have to be entered into Eqs. (28) and 
(34) are computed, following Kopal (1959), from the 


equations 
IM=])1(1—A™—C™ sin?9)+ B™ cos 


+ (11) (!—b™ cosp—y™ cos’0), (50) 
Ps P= (Iy—1))[(1— A+) 
+6™ cos—(C™+-y™) sin’9] 
—hM[(A*—A™)+ (C¥—C) sin] 
+ (B’— B™) cosd-+ (1-19) 
XL(a—a+-7— -y™) — (b—5™) cosé 
+ (y2—y™) sin’@], (51) 
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where the constants A, B, C, and their differences are 
the values obtained from the least-squares solution of 
Eqs. (48)—(49), while the coefficients a, b, and y cannot 
be determined from the observations and must be 
estimated theoretically; the methods of such estima- 
tions are described by Kopal (1959). In general, 
lo|+|b|+ly|«]4|+|B]+/C]. 


B. Correction of the Errors of Rectification 


Among the terms in Eq. (51) the terms J)\2(A**— A™) 
and (B’—B™) cos? can have absolute values much 
larger than any of the other terms except the first 
because the coefficients J)’* and cos@ can be near unity 
during the eclipse; the other terms in (51) contain 
factors such as sin*@ or J)\2@—J), which are in most 
cases much smaller than unity during the eclipse. Since 
A™—A™ and B’—B™ obtained from the solution of 
Eqs. (49) are not exact values we shall denote the true 
values of these coefficients by A'”*—A™+6,4 and 
B’— BM-+63. The inexact value ?2—/™ obtained from 
(51) should be replaced in (28) and (34) by the true 
value D2—]1—])6 4-6 cos@. The corrections 64 and 
dp should be treated as additional unknowns in these 
equations which now take the final form 


Occultation : 
Pe— [M+ (1—-1") (1-La™)Z 
co (ag? —a NES (agt— Xa?) 
+ (cea) Ea (yoyo) (1=P') 


X (1—La®) (Zo— Zo) +10b.4—5z cosd; (52) 
Transit: 
—M+ (1—P2) LZ, 
= a1 (2% — ary) + a2 (yo — yo™) : 
+[(1—-P)Lad/ (Zo—1)"](Zo— Zo) 
+1 )64—6, cos#. (53) 


In most cases the quadratic terms y—y can be 
neglected and Eqs. (52)-(53) can be solved by the 
method of least squares for the five unknowns: vq\'!—1,™, 
ayia, Zo—Zo, 54, and dp. 

It should be emphasized that the purpose of intro- 
ducing the additional unknowns 6,4 and 6, is not so 
much to increase the accuracy of determining the differ- 
ential limb darkening as to increase the reliability of 
the mean errors obtained from the least-squares solu- 
tion. The mean errors obtained by neglecting these two 
unknowns would be systematically in error. 

It should be noted that the term 


aed ed NC peta Be nay oles Ga 
in (51) is varying proportionally to the terms containing 
Zo in Eqs. (52)—(53). This means that the error of the 
term a@—a'+y"—-y™ contributes to the error of Zo 


and does not affect the other unknowns. 
If the solution is not made simultaneously for both 
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eclipses, the unknown 6g should not be determined. 
Considering only one eclipse, we have 


5p cosd~ (1—)2) sgn (cos#)ép+/” sgn (cos#)dz, (54) 


where sgn(cos@) equals +1 for primary eclipse and —1 
for secondary. The two terms on the right side of (54) 
are varying proportionally to the terms in (52)—(53) 
containing Zp and 64, respectively, which means that 
dz does not need to be determined separately. Making a 
solution of Eqs. (52) without determining 6, we obtain, 
instead of Z) and 654, the quantities Zo+[sen(cos@) 
X (6g +65) —Sary]/(1—La™) and 64—sgn(cos6)dz, where 
6, and da,, are corrections to the terms 6%2—6)™ and 
a—qM+7'’—-~yM, respectively. 

The contributions to )’—/™ resulting from the effects 
of spherical harmonic distortion of order higher than 
the second (cf. Kopal 1959) consist of terms varying 
proportionally to cos6, those varying proportionally to 
Jo2, and terms independent of phase. The terms of the 
two first types contribute to 64 and 6g, respectively. 
The terms independent of phase are eliminated together 
with 6g only in the case when the solution is not made 
simultaneously for both eclipses. These terms are’ in 
most cases negligibly small. 

The rectification eliminates only the part of the 
photometric terms resulting from the aspherical shape 
of stars. The remaining terms which can not be elimi- 
nated are called perturbations within minima. An 
estimate of the influence of these terms (as given by 
Kopal 1959) on ?#—/‘ was made for moderately 
distorted stars for which r2+7,=0.50. It was found that 
for such stars the errors of *—/™ due to perturbations 
within minima should not exceed 0.003. This means 
they are negligible for the present accuracy of photo- 
metric observations. The methods described in this 
paper should not, however, be applied to eclipsing 
binaries for which ra+r, exceeds, say, 0.55. This 
maximum value should be diminished if the difference 
in intrinsic colors of stars constituting a binary is very 
large since in this case the influence of pertutbations 
becomes also larger. 


C. Weights of Observational Equations and Additional 
Equations Describing the Observations outside of 
Eclipses 


Let us suppose that, as is usually the case, the dif- 
ferences /‘“—/]\ must be expressed in stellar magnitudes 
to make their mean errors independent of the brightness 
of the eclipsing binary. If the considered normal point 
is formed of observations of unit weight, the equation 
(52) or (53) corresponding to this normal point must 
be multiplied by n}/(M+P2), The equations multiplied 
by such factors will have equal weight. Let the color 
index, expressed in magnitudes and corresponding to an 
observation of unit weight made during the eclipse, be 
characterized by the mean error e. It can be shown that 
the mean error of the left side of Eq. (52) or (53) 
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Fic. 1. Coefficients ie —aU)k2 and (a&@—a¥) F for primary clip 
of YZ Cassiopeiae (transit). 


multiplied by the afore-mentioned weighting facto: 
equals (In10)«/5 or 0.4605¢. The mean error of this lef 
side is usually an order of magnitude larger than th 
product of an unknown on the right side and the mea1 
error of its coefficient. The errors of these coefficient: 
result from the uncertainty in determining a fron 
observed brightness by means of Eq. (38) or (39) 
If only the amplitude of light variations is much large: 
than the amplitude of color variations, the errors of the 
coefficients of the unknowns in Eqs. (52)-(53) haw 
negligible influence on the results of solution. Thi 
condition is essential for validity of the method pre 
sented in this paper. 

In order that the unknowns 64 and éz may be reliably 
determined, Eqs. (52)—(53), following the procedur 
first introduced by Piotrowski (1948), must be supple 
mented by two additional equations containing infor 
mation about the accuracy with which AY—A™ anc 
B’— B™ are determined from the harmonic analysis o 
the observations outside of eclipse. These will be equa 
tions of the form 64=0 and 6g=0 multiplied by sucl 
factors that the weight of these equations is the sam 
as that of all the other equations of the form of (52)—(53 
multiplied by the weighting factors. Let the mean error 
of A’—A™ and B’— B® as obtained from the least 
squares solution of Eqs. (49) be esa and exp, respec 
tively. Now the two additional equations are 


0.4605 (€/ena)d4=0, 
0.4605 (€/ern)dn=0; 


(55 
(56 


we can easily see that if the mean errors of 54 and 6, 
in Eqs. (55)-(56) are esa and eaz, the error of the lef 
side of each of these equations equals 0.4605«. 


IV. EXAMPLE, DIFFERENTIAL LIMB DARKENING 
OF YZ CASSIOPEIAE 


The eclipsing binary YZ Cassiopeiae was selected a 
an example because the coefficients of limb darkenin; 
have been determined for this system with an accuracy} 
much greater than for any other eclipsing binary. Th 
binary was observed by Kron (1939, 1942) in blue anc 
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TaBLE II. Results for YZ Cassiopeiae. 


Obtained 


4644500 — 47,8700 74,4500 _ 76700 4500 6700 
= 7, 4500__ 4, 6700 
Ua Ua 


Zo B4500_— B6700 


Assumed 
Solution k Toa tq ®00 WB hay 
I 0.540 0.49 
0.45 0.063 
II 0.532 0.60 
Kron (Kopal, Shapley) 0.542 0.49 0.45 0.063 
+.012 +.06 +.16 +-.002 


0.330 0.368 —0.05 0.438 +0.004 

+.013 +.014 Safi +.012 +.001 

0.325 0.342 —0.09 0.435 +0.004 

+.011 =e 012 +.24 + .010 +.001 
0.16 —0.03 0.439 +0.003 
= .08 +.24 


red colors (effective wavelengths 4500 and 6700 A) and 
independent solutions were made by this author for the 
observations in each color. The geometrical elements 
obtained from these solutions are slightly different. 
According to Kopal and Shapley (1956) the quadratic 
terms in the law of limb darkening (our coefficients 2) 
must be taken into account if the observations in both 
colors have to be explained by the same set of geo- 
metrical elements. 

This opinion is not supported by the present analysis, 
which shows that the observations in both colors can 
be well explained by a single set of geometrical elements 
without introducing the quadratic terms in the law of 
limb darkening. The coefficients a?—a¥ of the unknown 
y"—y in (10) and (17) are, in the case of YZ Cassio- 
peiae, almost exactly proportional to the coefficients 
a?—a” of the unknown x—x™. This is shown on 
Figs. 1 and 2, where both these coefficients are plotted 
as a function of . 

The proportionality of the coefficients of y'—y* to 
those of «1— x excludes the possibility of independent 
determination of these two unknowns. For that reason 
the solution of Eqs. (52)—(53) for Kron’s observations 
of YZ Cassiopeiae was made on the assumption that 
yM=y=(. For estimating the influence of the assumed 
values of & and of the blue coefficient of limb darkening 
on the final results two least-squares solutions were 
made, each for a different set of initial elements. The 
first set of assumed values of k, La**™, and u*™™ is almost 
identical with values given for YZ Cassiopeiae in the 


0.00 


oa | -05 ° 
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Fic. 2. Coefficients a@?—aY and a?—aY for secondary eclipse of 
YZ Cassiopeiae (occultation). 


catalogue by Kopal and Shapley (1956) and with values 
obtained from the blue observations by Kron (1942). 
The other set is based on the value of the radius of the 
smaller star obtained by Kron from the observations in 
the red light only. This second set of elements does not 
represent the blue observations as well as the first set; 
the mean O-C is about 30% larger. 

The results of both solutions are given with their 
mean errors in Table II. We can see that a rather large 
change of the initially assumed elements has very little 
influence on the values of «4°°—46 and Zp» obtained 
from the solution. The difference u**— 1° is somewhat 
more strongly affected because it depends explicitly on 
the assumed value of u*°. The theoretical changes of 
color computed from the results of both solutions are 
compared in Fig. 3 with values of /5°—/]* resulting 
from Kron’s observations. The curves computed from 
both solutions almost coincide. This means that the 
mean errors given in the upper part of Table II should 
be increased only slightly to take into account the 
uncertainty of the assumed values of k, u*°, and 1,45”. 

The most probable value of the difference of limb 
darkening coefficients for star eclipsed in primary 


minimum 
144500 146700 — ().35+0.02 (m.e.) (57) 


is about two times larger than the value resulting from 


00 


-02 


Secondary 


‘00 05 10 5 


d 
Phase ~ 
Fic. 3. Color indices of YZ Cassiopeiae according to obser- 
vations made by Kron. The solid line represents solution I; the 
dashed line, solution II. Open circles are the observed values for 
primary eclipse, crosses for secondary eclipse. 
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the independent solutions made by Kron (1942) for 
each color separately. The value (57) is larger than the 
theoretical values of this difference for a star of spectral 
type A3 which are 0.29 according to Miinch and 
Chandrasekhar (1949), 0.23 according to Hosokawa 
(1955), and 0.27 according to van’t Veer (1960). 

The solutions for investigating the differential limb 
darkening of several other well observed eclipsing 
binaries will be made by the author with the aid of an 
electronic computer. The results will be published as 
the second part of this paper. 
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Photoelectric observations of 171 yellow and 169 blue magnitudes of V Tucanae are presented. The 
available evidence indicates a constant period. The system was brighter on two nights than on the others. 
The preliminary elements and their permissible ranges are computed. 


I, INTRODUCTION 


HE variability of V Tuc (CPD—72°69, 185.1907, 
HV 3021, HD 5148) was discovered by Miss 
annon (Pickering 1907, 1908). Further observations 
ere made in 1917-18 (Dawson 1918) and the period of 
ie variable was found to be 0.8709 day. Dawson also 
und a length of minimum to be about 0°18 and the 
»pth 1.8 mag. Thearbitrary scale he used did not permit 
ymputation of the physical elements of the system. 
aposchkin (1947) included the system in his general 
irvey of eclipsing systems by photographic estimates. 
liss Cannon (Cannon and Pickering 1918) classified 
ie spectrum as A2. Improved light elements and a 
ary brief description of the light curve, based on the 
resent observations, have been given by Wood (1959). 


II. OBSERVATIONS 


The observations were made by the senior author 
| the Mount Stromlo Observatory of the Australian 
‘ational University, using the Reynolds 30-inch re- 
ector and a 1P21 multiplier photocell refrigerated 
ith dry ice. The filters have been described by 
ascoigne and Burr (1956). The spectral ranges cor- 
‘spond closely to the International P and V system. 
lhe linearity of the direct-current amplifier has been 
wefully checked in the laboratory of the Mount 
tromlo electronics staff; rougher checks on the stars 
towed that any variation was less than the error of 
bservation. 

The observations in Table I are expressed in helio- 
tntric Julian days, with phases from the elements 
iven by Wood (1959). 

| 


Pr. Min.=JD 243 6139.142+-0987091649 FZ. = (1) 


‘he comparison star was BD—72°70. Correction for 
ifferential extinction is no case exceeded 07002. 
lecause of the rapidity of the light changes during 
clipse, it was found desirable to take two readings in 
ach color on the variable star between successive ob- 
arvations of the comparison; in order to make the 
bservations as uniform as possible, the same practice 
‘as carried out on all parts of the light curve. 

Careful inspection on the best nights showed no 
ebulosity or faint stars in the region around the vari- 
ble included in the diaphragm used. However, difficulty 
fas encountered in attempting to take sky readings 


nearby in that the sky brightness did not repeat from 
region to region and was consistently brighter than the 
sky readings near the comparison star. A photograph 
taken by Dr. A. R. Hogg with the 74-inch Mount 
Stromlo reflector clearly shows nebulosity in some of 
the regions where sky readings were attempted. Hence, 
for purposes of reductions, the sky readings near the 
comparison star were interpolated to the time of the 
observations of the variable. Further justification for 
this treatment will be found in the discussion of the 
reliability of the solution. 


Ill.. EPOCH AND PERIOD 


The published times of minimum are listed in Table 
II, together with residuals from the new light elements 
given in (1). There is no evidence for change of period 
but the evidence is so meager that small changes could 
easily have gone undetected. 


IV. DISCUSSION 


When the reflected yellow observations between 
minima were plotted, they fell into two clearly distinct 
sets, separated by about 0.07 mag. A similar situation 
existed in the blue, except that the separation was only 
about 0.04 mag. and the division was not as clearly 
shown. The observations comprising the upper curve 
were all taken on two consecutive nights; the system 
was obviously brighter on these two nights than on the 
others and the excess brightness continued through 
secondary. 

Obviously, any attempt to make a rectification based 
on the average of the observations would produce a 
meaningless solution which would fit neither set satis- 
factorily. Therefore, four separate rectifications were 
carried out, i.e., one in each color for both the upper 
and lower curves. Rectification coefficients were com- 
puted for the upper and lower curves separately. Rec- 
tification of each followed the standard procedure 
(Russell and Merrill 1952). In the conventional equation 


IT=Ao+A1cosB+A. cos20+ B, sind+B, sind, (2) 


the following values were found from a least-squares 
solution : 


Upper blue curve, A =0.9522+-0.0011, 
Ai=—0.0362+0.0007, 


A2=—0.0286-0.0022, 
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TABLE Ia. Yellow observations of V Tucanae. 

JDyHel. JD Hel. JD Hel. ot 
2430000+ Phase m 2430000-++ Phase m 2430000-++ Phase m | 
6134.9839 0.1965 —0.309 6136.1977 0.5394 —0.297 6165 .9871 0.7176 —0.38) 

9881 0.2007 0.318 6137.9576 0.5574 0.399 9984 0.7289 0.38 
9971 0.2097 0.299 9698 0.5696 0.377 6166 .0038 0.7343 0. 
6135 .0008 0.2134 0.314 6138.0244 0.6242 0.388 .0306 0.7611 0. 
.0105 0.2231 0.310 .0291 0.6289 0.389 0349 . 0.7654 OF 
.0138 0.2264 0.309 0399 0.6397 0.379 .0465 0.7770 0.4 
.0318 0.2444 0.317 9943 0.7232 0.391 .0515 0.7820 0. 
0352 0.2478 0.320 9982 0.7271 0.390 .0636 0.7941 0.45}, 
0435 0.2561 0.317 6139 .0088 0.7377 0.349 .0676 0.7981 0.48) 
04814 0.2607 0.308 .0148 0.7437 0.419 .0792 0.8097 0.57" 
.0588 0.2714 0.308 .0250 0.7539 0.389 0839 0.8144 0.63) 
.0642 0.2768 0.357 .0293 0.7582 0.398 0877 0.8182 0.68 
0741 0.2867 0.304 0415 0.7704 0.407 .0989 0.8294 0.89) 
.0805 0.2931 030075: .0466 0.7755 0.410 1034 0.8339 1.02} 
.0924 0.3050 0.319 0968 0.8257 0.840 . 1084 0.8389 1.091)" 
.0976 0.3102 0.328 . 1002 0.8291 0.910 1281 0.8586 1.93} 
. 1092 0.3218 0.289 . 1082 0.8371 1.084 1316 0.8621 2.128 
.1145 0.3271 0.292 .1122 0.8411 Nol) .1385 0.8690 2.71) 
. 1396 0.3522 0.328 1154 0.8443 te o22 . 1433 0.0029 2.664) 
1435 0.3561 0.330 1253 0.8542 1.649 1477 0.0073 2.319 
1552 0.3678 0.360 1291 0.8580 1.837 . 1530 0.0126 1.96,|| 
1594 0.3720 0.374 . 1336 0.8625 2.144 1576 0.0172 1.73) 
.1703 0.3829 0.367 1410 0.8699 2.676 1701 0.0297 1.213} 
.1760 0.3886 0.380 1518 0.0098 2.079 ails! 0.0347 1.0932) 
1888 0.4014 0.395 1562 0.0142 1.868 1778 0.0374 1.003} 
1937 0.4063 0.405 .1601° 0.0181 1.681 1816 0.0412 0.855) 
2055 0.4181 0.411 . 1646 0.0226 1.498 1951 0.0547 0.67 
2107 0.4233 0.410 "1684 0.0264 1.374 . 1989 0.0585 0.62 
52255 0.4381 0.412 . 1808 0.0388 0.938 . 2038 0.0634 0.57. 
2297 0.4423 0.418 1857 0.0437 0.861 6169.0534 0.3002 0.39 
.9200 0.2617 0.305 . 1900 0.0480 0.775 .0587 0.3055 0.38 
.9247 0.2664 0.292 2035 0.0615 0.580 .0716 0.3184 0.387) 
.9374 0.2791 0.282 2079 0.0659 0.554 .0776 0.3244 0.383) 
9417 0.2834 0.313 2119 0.0699 0.545 .0910 0.3378 0.389), 
9541 0.2958 0.316 6165 .0332 0.6346 0-.378 .0965 0.3433 0.384) 
9591 0.3008 0.315 .0499 0.6513 0.366 6170.9490 0.4540 0.470) 
.9691 0.3108 0.335 0599 0.6613 0.361 .9619 0.4669 0.446) 
9744 0.3161 0.324 .0707 0.6721 0.370 .9730 0.4780 0.417) 
9872 0.3289 0.322 .0757 0.6771 0.372 9847 0.4897 0.417) 
.9930 0.3347 0.312 0868 0.6882 0.367 .9966 0.5016 0.399)) 
6136.0396 0.3813 0.368 .0912 0.6926 0.353 6171.0077 0.5127 0.379) 
0434 0.3851 0.366 1203 0.7217 0.379 .0201 0.5251 0.375} 
.0539 0.3956 0.393 1256 0.7270 0.369 .0385 0.5435 0.367}, 
0583 0.4000 0.401 1412 0.7426 0.385 .0515 0.5565 0.370) 
.0714 0.4131 0.390 1535 0.7549 0.407 .0666 0.5716 0.384} 
.0777 0.4194 0.409 . 1580 0.7594 0.414 .0808 0.5858 0.365} 
.0906 0.4323 0.430 eh 22 0.7736 0.426 6187 .9604 0.0471 0.751) 
.0949 0.4366 0.434 ~ .1780 0.7794 0.438 -9657 0.0524 0.677 
1061 0.4478 0.427 *.1922 0.7936 0.473 .9767 0.0634 0.549) 
.1110 0.4527 0.426 1982 0.7996 0.500 .9809 0.0676 0.519) 
. 1238 0.4655 0.390 2040 0.8054 0.562 .9926 0.0793 0.426 | 
. 1300 0.4717 0.378 9382 0.6687 0.363 .9966 0.0833 0.414) 
1585 0.5002 0.324 9494 0.6799 0.354 6188 .0072 0.0939 0.409 | 
. 1635 0.5052 0.331 - .9540 0.6845 0.363 6221.9803 0.1012 0.440 | 
ESV 0.5174 0.323 9641 0.6946 0.371 9890 0.1099 0.400 | 
. 1806 0.5223 0.319 -9688 0.6993 0.372 6222 .0023 0.1232 0.376 | 
1920 0.5337 —0.309 9827 0.7132 —0.386 -0111 0.1320 —0.398 
- SS 


Lower blue curve, 


Lower yellow curve, 


Upper yellow curve, 


Ao=0.9739-+0.0003, 
A,=—0.0105+0.0014, 
A2=—0.0274+0.0006, 


Ao=0.9614-£0.004, 
A1= —0,0091-+0.0029, 
A»=—0.026-+0.0000, 


Ay=0,9332-£0.0003, 


A,=—0.0455+0.0034, 
A»=—0.0405+0.0006. 


For each of these rectifications, 
attempted, following the procedure of Russell an 
Merrill (1952) and using the Merrill nomograms ae 
tables (Merrill 1950, 1953). Because of the shalloy 
depth of secondary, the solution was made using th 
x function derived from the shape of primary minim 
and the usual relationship found from the rectifi 
depths of the minima. i‘ 


solutions wer 


When the observations were rectified by the coeff 
cients found in the upper curve in each color, no solutior 
was found. The depth lines and the shape lines did no 


itersect. Solutions were found using observations 
ctified by the coefficients found from the lower curves. 
his suggests that on the nights when the system was 
nusually bright, factors other than those covered by 
le conventional theory were present. 

If we assume that light was added to the system by 
temporary shell or stream to form the upper curve, 
e can then carry through solutions after correcting 
r this added light. (The fact that the increased light 
as greater in the yellow than in the blue argues 
yainst "its being a thermal effect connected with 
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TABLE Ib. Blue observations of V Tucanae. 
JD Hel. JD Hel. 

Phase m 2430000-+- Phase m 2430000+- Phase m 
0.1976 —0.107 6136.1987 0.5404 —0.094 6165 .9850 0.7155 —0.134 
0.2017 0.103 6137.9589 0.5587 0.180 -9887 0.7192 0.146 
0.2107 0.098 .9710 0.5708 0.171 6166.0005 0.7310 0.154 
0.2141 0.104 6138 .0253 0.6251 0.144 .0074 0.7379 0.149 
0.2239 0.107 .0302 0.6300 0.142 .0326 0.7631 0.163 
0.2271 0.103 .0409 0.6407 0.136 .0367 0.7672 0.168 
0.2452 0.097 -9951 0.7240 0.158 .0489 0.7794 0.169 
0.2485 0.110 .9989 0.7278 0.157 .0533 0.7838 0.170 
0.2569 0.093 6139 .0099 0.7388 0.140 .0657 0.7962 0.223 
0.2614 0.087 .0161 0.7450 0.180 .0696 0.8001 0.251 
0.2724 0.089 .0258 0.7547 0.158 .0818 0.8123 0.378 
0.2778 0.132 .0306 0.7595 0.167 .0857 0.8162 0.423 
0.2878 0.095 .0428 0.7717 0.173 .0899 0.8204 0.527 
0.2938 0.098 - .0476 0.7765 0.160 . 1008 0.8313 0.724 
0.3067 0.123 .0976 0.8265 0.639 . 1060 0.8365 0.867 
0.3114 0.134 . 1009 0.8298 0.687 .1101 0.8406 0.973 
0.3229 0.104 . 1090 0.8379 0.879 .1295 0.8600 1.903 
0.3281 0.106 oleae 0.8416 1.023 .1339 0.8644 2.382 
0.3530 0.124 .1167 0.8456 eo .1409 - 0.0005 2.943 
0.3570 0.128 .1260 0.8549 1.539 . 1449 0.0045 2.590 
0.3685 0.155 .1298 0.8587 1.797 . 1508 0.0104 2.048 
0.3730 0.154 . 1346 0.8635 2.187 n552 0.0148 1.769 
0.3839 0.158 .1416 0.8705 2.940 .1598 0.0194 1.483 
0.3896 0.159 .1527 0.0107 1.981 sey? 0.0328 0.935 
0.4024 0.168 .1568 0.0148 1.709 .1796 0.0392 0.782 
0.4073 0.188 . 1610 0.0190 1.483 . 1842 0.0438 0.674 
0.4191 0.188 .1655 0.0235 1.293 1971 0.0567 0.451 
0.4243 0.197 . 1693 0.0273 1.109 .2031 0.0609 0.371 
0.4392 0.203 .1818 0.0398 0.736 . 2065 0.0661 0.322 
0.4429 0.198 . 1868 0.0448 0.668 6169.0561 0.3029 0.148 
0.2626 0.092 . 1909 0.0489 0.549 .0610 0.3078 0.154 
0.2674 0.092 . 2047 0.0627 0.415 .0745 0.3213 0.141 
0.2800 0.084 . 2088 0.0668 0.320 - 0800 0.3268 0.158 
0.2844 0.093 .2131 0.0711 0.276 -0936 0.3404 0.142 
0.2968 0.117 6165 .0353 0.6367 0.144 .0988 0.3456 0.145 
0.3016 0.107 -0519 0.6533 0.115 6170.9527 0.4577 0.204 
0.3118 0.126 .0625 0.6639 0.108 -9645 0.4695 0.219 
0.3171 0.108 .0737 0.6751 0.128 .9754 0.4804 0.173 
0.3302 0.115 .0774 0.6788 0.126 .9869 0.4919 0.180 
0.3357 0.115 .0890 0.6904 0.120 -9990 0.5040 0.171 
0.3821 0.141 -0935 0.6949 0.119 6171.0110 0.5160 0.133 
0.3859 0.159 .1236 0.7250 0.147 .0221 0.5271 0.137 
0.3965 0.175 .1294 0.7308 0.140 .0412 0.5462 0.127 
0.4010 0.167 .1438 0.7452 0.147 .0546 0.5596 0.136 
0.4143 0.174 ES) 0.7569 0.158 .0712 0.5762 0.121 
0.4208 0.180 .1610 0.7624 0.156 .0832 0.5882 0.129 
0.4331 0.194 .1744 0.7758 0.180 6187.9636 0.0503 0.513 
0.4378 0.202 .1811 0.7825 0.178 .9677 0.0544 0.440 
0.4489 0.196 p .1952 0.7966 0.228 .9787 0.0654 0.327 
0.4536 0.192 2011 0.8025 0.280 .9830 0.0697 0.287 
0.4666 0.160 . 2063 0.8077 0.325 .9946 0.0813 0.196 
0.4727 0.165 .9410 0.6715 0.129 -9987 0.0854 0.191 
0.5010 0.129 .9519 0.6824 0.126 6221.9857 0.1066 0.159 
0.5064 0.122 .9559 0.6864 0.118 -9916 0.1125 0.157 
0.5186 0.099 .9664 0.6969 0.138 6222 .0053 0.1262 0.159 
0.5235 0.112 .9719 0.7024 —0.138 .0081 0.1290 0.176 
0.5352 —0.111 


intrinsic changes in either component of the variable 
in the comparison star.) The assumption of added light 
equal to from 0.06 to 0.07 of the total light in the yellow 
and 0,02 to 0.04 of the total in the blue permits us to 
find solutions in which the elements fall well within 
the range of those given by the observations rectified 
from the lower curves. 

In actual practice, uncertainties in both the depths 
and shapes of the minima, as given by the scatter and 
distribution of the observations, were determined and 
a number of solutions were carried through for assumed 
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TABLE II. Observed minima.* 


WOOD AND Gy EL MeGLUSKE Y; 


O 

E JD 2400000-++ O-C Reference 
—22 363 16662.851 +0.014 Dawson (1918) 
—21 851 17108.743 —0.003 Dawson (1918) 
—21 480 17431.833 —0.023 Dawson (1918) 
—16 784 21521.678 —0.002 Dawson (1918) 
—16 783 21522.551 0.000 Dawson (1918) 
—16777 21527.779 +0.003 Dawson (1918) 
—16775 21529.518 0.000 Dawson (1918) 
—16 767 21536.483 —0.002 Dawson (1918) 
—16 761 21541.712 +0.001 Dawson (1918) 
—16 760 21542 .582 0.000 Dawson (1918) 
—16 752 21549 .540 —0.009 Dawson (1918) 
—16 690 21603 .547 +0.001 Dawson (1918) 
—16 684 21608 .765 —0.006 Dawson (1918) 
—16 644 21643 .609 +0.001 Dawson (1918) 
—16 629 21656.671 —0.001 Dawson (1918) 
—16 628 21657 .543 0.000 Dawson (1918) 
—16 451 ~21811.648 —(.047 Dawson (1918) 
—16 450 » 21812.564 —0.002 Dawson (1918) 
— 7828 29321.598 —0.010 Gaposchkin (1947) 

0 36139. 142 0.000 Wood 

Sie 31 36166. 140 0.000 Wood 


a The first three minima are as given by Dawson from the times listed 
by Pickering. The original publication included an error in one date which 
was corrected by Dawson. The first minimum is the mean of two times 
on the same night at which the star was listed as ‘‘faint.”’ 


degrees of limb darkening, «=0.4, 0.6, and 0.8. The 
possibilities were then reduced by considering only 
those which give common elements (k, 7, 7,). Table III 
lists the elements from the solution we consider the 
most satisfactory. Table IV is more realistic and lists 
the permissible ranges. 

The weakest part of the whole solution depends 
on the observations in secondary minimum. Observa- 
tions at the bottom of secondary were made on three 
nights; two of these were the two nights when the 


TR 
TaBLE IIT. x=0.4 (assumed). V mag. at max. 10"7) 
A; (Y) 2™40 A, (B) 2™85 
Az (Y) om11 Ao (B) 0=10 
ag’ 0.975 : 
ao™ 0.979 
IE ROD) 0.941 Ly (B) 0.979 
t 88°23 
rq 0.298 
Ts 0.286 


night alone. The computed curve thus appears to ¢ 
a very poor fit in secondary to the observations aj 
whole. Fortunately, two factors make this less seril® 
than it otherwise would be. The first is that the obj 
vations on the two nights when the star was brigh> 
run consistently above the curve throughout —_ 
by the same amount as they do outside eclipse. 
additions of an amount of light from the out-of-ecli 
observations will give a satisfactory fit to these. 7 
second is that if secondary is appreciably shallower tl] 
the value used no solution can be found; if it is apr 
ciably deeper, it will move the solution nearer td 
total eclipse. Since the sharp bottom of primary sho} 
that this is a partial eclipse and since ap (ratio of tj 
light lost at the bottom of minimum to that whi 
would be lost at internal tangency) is 0.98, little cal 
adjustment in this direction is possible. The same li 
of argument indicates that little, if any, of the nebulosi} 
in the general region was close enough to the varial} 
to be included in the diaphragm and that the meth) 
of correcting for the sky background was the prog 


i] 


i 


i 


i] 


Fic, 1. Blue light curve of 
Tuc. 


phase 


O6 


Taste IV. Limiting solutions. 


+ 
x=0.4 x=0.6 
| 0.96 to 0.98 0.95 to 0.97 
#| 0.30to 0.29 0.30 to 0.29 
4 88°3 to 89°3 88°3 to 89°0 
| 0.0 for lower curves 0.0 for lower curves 
» 0.04 for upper blue curve 0.02 for upper blue curve 
' 0.07 for upper yellow curve 0.06 for upper yellow curve 
= | 
i} 


ie. In this case it is the depth of primary which would 
"; most seriously affected by systematic errors in 
eatment of sky background. The result of this would 
‘sain be to alter the depth function, and the objections 
:“ cited apply to any significant change in this. 

" In view of the uncertainties of rectification, it seems 
") be unwise to carry the discussion further in detail at 
‘resent. This is especially true, since Dr. A. R. Hogg 
‘ready has a number of three-color photoelectric ob- 
‘rvations and with J. S. Sommerville is continuing 
bservation of the system. The large difference of 
‘iminosity between the components suggests a rela- 
lvely large mass ratio. This is borne out by the fact 
‘hat only a small ellipticity effect is shown by the ob- 
‘ervations between eclipse, even though the system is a 
lose one. This again indicates only minor distortion of 
he larger component and this in turn requires a large 
hass ratio. It is possible to estimate this ratio from the 
ight curve alone (Wood 1960), but we feel that even 
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an approximate numerical value should await a curve 
obtained when intrinsic variations were not present. 
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Mean densities of 52 visual binaries and 82 eclipsing variables have been calculated without employing 
parallax values by using a method originally proposed by Opik. The visual double stars seem on the average 
to be 10 times denser than photometric binaries, and their mean density is normally similar to the sun’s 
density. Very likely, however, the difference between densities in the two types of binary systems is the 


result of observational selection. 
INTRODUCTION 


HE mean density p of a visual binary system can 
be known even though the parallax is unknown. 


Tn fact, in the determination of p from the masses and ~~ 


volumes of the binary stars, the parallax is automati- 
cally eliminated. This fact was brought out as long ago 
as 1916 by Opik, who derived a relation giving the 
mean densities expressed in terms of the elements of the 
orbit, of the mass ratio, and of the surface brightness. 
Opik obtained this last ee depending on the absolute 
temperature, from the spectral type; for the mass ratio 
he adopted some hypothetical values empirically de- 
pending on the difference of magnitude of the com- 
ponents. Using these methods he computed the mean 
density for the principal component of each of 39 
binaries. The results were necessarily only rough 
approximations, because the data were more or less 
seriously affected by the inadequacy of the above con- 
versions as well as by errors. 

Bernewitz (1921) and Thompson (1929) later pub- 
lished similar investigations concerning the density of 
binary stars. They attempted to derive carefully the 
relation between surface brightnesses and temperatures 
and the values of the mass ratios. The agreement 
between their results is not good. The principal cause 
of such deviations is to be sought in the uncertainty 
inherent in the derivation of the temperature scale from 
the spectral types. Also, for computing the mean density 
of the principal component, the combined spectral type 
of the system was considered, which is not always that 
of the single stars. 


DETERMINATION OF THE DENSITY 


Giannuzzi (1958) has examined the problem again 
in order to give a formula that is easier to apply for the 
mean density p of a double star. 

From Kepler’s third law, taking for units the solar 
mass and the year, and with the parallax p and semi-axis 
major a, both in seconds of arc, 


Mat po= a°/p P?. (1) 
The volume measured in terms of the sun’s volume is 
VitVo= (Ri/Ro)?+ (Ro/Ro)3. (2) 


Moreover, the radius R of a star is related to the ab- 
solute bolometric magnitude M and to the effective 
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temperature T by the Stefan formula for blackboj 
radiation : : 


Mo—M =S5\og(R/Ro)+10 log(T/To). 


From the absolute magnitude M and apparent maj 
nitude m, in the standard formula, q 


—M=—m—5—5 logp | 
=5 log(R/Ro)+10log(T/To)—Mo, (| 


from which 


—,) 


Ro tbe 
—--(—) 10-10 .2Mo—-0.2m_ ¢ 
Ro p\To 


Introducing this result in (2) we have 


10-3+0.610 100-61 10—0-6m2 
VitV2= aah 
P U(Ts/To)* (after 
whence the mean density p is | 
mate _ Co 1Q-0-671 {Q-0-6m2 4-1 : 
| . +——_| (i 
ts ae Pp? (T1/To)$§ (T2/To)§ 


where C= 10?-°-%0=1.45, if we assume Stebbins an 
Kron’s (1957) value of Mee +4.73. For To we wi 
assume Unséld’s (1955) value of 5780°K. 


TEMPERATURES AND ACCURACY OF DATA 


‘In (7), p is a function of the sixth power of 7 
therefore a small error in T greatly affects p. But thi 
difficulty is not unique to our method, since in the usua 
computation of stellar density the radius (indispensabl 
in calculating the volume) depends on the second powe 
of T; and hence T is present in the calculation of p t 
the same sixth power. Therefore it is necessary to hay 
values of the temperature as exact as possible. 

The stellar spectra from which temperatures wer 
derived usually refer to the combined light of th 
system ; colorimetric observations of single component 
are very few. ; 

The modern colorimetric measures of visual doubl 
stars have been made by Wallenquist (1950) (about 2! 
systems), by Muller (1952) (51 binaries), and b: 
Hopmann (1953-1956), who has observed the singl 
components of about 700 binaries, but until now ha 
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| p.e. in T p.e. in p 
} +850° +0.50p 
+340 +0.34p 
| +120 +0.12p 


jublished the data for only 400 of them. Giannuzzi also 
ook some part in Hopmann’s observations at Vienna 
h 1955. The values for T from Hopmann’s list, in view 
if its homogeneity and much greater size, have been 
ised in the present paper. 

|| Accurate colorimetric observations must be corrected 
or differential atmospheric extinction, for changes in 
‘he color of the photometric lamp from night to night, 
tnd for the observer’s nightly personal equation. 
(Therefore Hopmann (1953) corrected each measure by 
in appropriate amount calculated from the theory of 
itmospheric extinction, and always compared the color 
bf each star with that of many stars of known color, 
situated in the same region of the sky, and assumed as 
colorimetric fundamental stars. He then reduced his 
measures to Becker’s scale of radiation temperatures, 
and gave for each component of each double star a 
value of c/T,, where c has the value 14 320° and 7, is 
the radiation temperature. The average mean error of 
a c/T, measure, according to Hopmann, is +0.147; 
accordingly the percentage error in 7, is higher for 
hotter stars. 

_ From radiation temperatures we have derived 
efiective temperatures T according to Kuiper’s (1938) 
‘conversions. According to Becker (1952), such a scale 
‘introduces an error of about +1000° for the highest 
and of about +100° for the lowest temperatures; that 
lis, the uncertainty of the scale is no more than 5% of 
its value. 

_ In formula (7), putting 71= 7» and taking the loga- 
rithmic derivative, we find dp/p=6 dT/T. This formula 
is not quite correct, since the bolometric magnitudes 
also depend on 7, but we think that the discrepancy in 
‘temperature is not likely to alter the bolometric cor- 
rections appreciably. The probable errors for different 
values of 7 are given in Table I. 

In order to determine the distribution of double stars 
according to their mean densities, we take logp in 
steps of 0.50. Our table of errors suggests that even 
under the worst conditions any binary system will have 
been placed in the interval to which it belongs or in 
that beside it. Moreover, only nine systems of those 
studied by us have temperatures higher than 9000° 
and hence the highest probable error; the temperatures 
of other systems, ranging from 3500° to 9000°, are 
sensibly more correct. Thus the general distribution of 
double stars according to their densities should be very 
nearly correct. 

By adopting Unsdéld’s value for the sun’s effective 
temperature, we have set the zero point of our system 
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such that dwarf stars of type G2 have effective tem- 
peratures of 5780°K. Should the zero point change, the 
T values change accordingly. In particular, again 
according to Unsdld (1955), if from a system of T 
values defined on 7% as base value one passes to a 
system of 7’ values defined on 7,’ as base value, one 
has the difference 


CaP) Paty tr 


Thus a change of -+-100° in the value of the sun’s 
effective temperature would alter all values by little 
more than 1.5%, so that the densities should vary by 
+0.09p. 


MAGNITUDES AND ORBITAL ELEMENTS 


To determine the apparent bolometric magnitudes 
we have used Kuiper’s (1938) paper, where, for stars 
hotter than 7000°, tables of bolometric corrections are 
constructed using theoretical spectral energy curves, and 
for later types, bolometric corrections are derived from 
radiometric observations. Our work has pursued the 
following plan: 


(a) Obtain the photovisual magnitude m in the 
Harvard Photometric System (Harvard Mimeo. Ser. 
III, 1, 2, 1938) which usually gives the combined mag- 
nitude of a binary system. 

(b) To find the magnitude difference Am between the 
separate components, we consulted the following 
catalogues: Kooreman, Bull. Astron. Insi. Neth. 10, 
No. 375; Muller, Ann. Strasborg 5, Heft 3; Miczaika, 
Veroff. Heidelberg-Kénigstuhl 15, No. 6; Detre, Astron. 
Nach. 273, No. 253; Hopmann, Mitt. Wien 6, No. 8. 

We have not used Wallenquist’s (Uppsala) recent 
photometric measures since his observations show 
systematic differences with the observations of other 
catalogues. When several determinations were available 
for a stellar system, we averaged Am, since Hopmann 
has shown that the above sets of measures have no 
systematic difference between them. For a few stars, 
for which no magnitude differences are available in 
these catalogues, we adopted the values of Harvard 
Photometry or those given in Baize’s catalogue of 
orbits. 

(c) Determine m; and mz by the relations 


mMm=Mm+x, Me=m+Am, 
x= 2.5 logio(A-+2.512-4”). (8) 


(d) Reduce m; and mz expressed on the Harvard 
Photovisual scale to International Photovisual mag- 
nitudes on the system of the North Polar Sequence 
1922, by the formula in the Harvard Mimeograms 
referred to above: 


Mrpy= MyPvt+0.10+0.06C ;, 


where C; is the color index, known from the tem- 
perature. 
(e) Finally, apply the bolometric corrections. 
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TaBLE II, Visual double stars. 
Int. pv 
Name my m2 Ti is p logp References 
588 7.82 9.66 6160° 4700° 1.00 0.00 Rabe. 1948, Astron. Nach. 276. 
755 6.45 6.89 4260 3950 0.014 —1.85 van den Bos. 1938, Union Obs, Circ. No. 100. 
0.013 —1.89 Baize. 1946, Astron. J. 51. 
0.015 —1.82 Rabe. 1951, Astron. Nach. 280. 
0.015 —1.82 Rabe. 1955, I.A.U. Circ. 6. 
0.014 —1.85 Muller. 1956, 7.4.U. Circ. 10. 
1538 6.88 6.93 6430 7200 0.96 —0.02 Volet. 1933, J. Obs. 16. 
1.64 0.21 Palacios. 1954, Urania 37. 
1615 4.40 5.35 15400 10100 0.47 —0.32 ‘ Rabe. 1954, I.A.U. Circ. 2. 
2034 8.5 Ont 5450 5150 1.30 0.11 Pensado. 1954, Bol. Madrid 4. 
2122 PP RSHE Fees lke) 6010 4700 0.78 —0.11 Rabe. 1954, 7.4.U. Cie. 2. 
2416 15 8 6120 6300 2.44 0.39 Dommanget: 1954, I.A.U. Circ. 1. 
3264 5.91 8.44 9200 6080 1.42 0.15 Kuiper. 1937, Astrophys. J. 86. 
3711 5.94 6.69 9180 9180 0.66 —0.18 van den Bos. 1933, Union. Obs. Circ. No. 90. 
1716 0.06 Voronov. 1933, Tashkent. Obs. Circ. No. 10. 
1.46 0.16  Baize, 1953; J. Obs. 36. 
5400AB 5506! 46.20 7270 8500 0.10 —1.00 Broche. 1956, 7.4.U. Cire. 8. ee 
5871 ULE Theee? 4470 5250 0.17 —0.77 Van Biesbroeck. 1925, Monthly Notices Roy. Astron. Soc. §\\— 
Op22 —0.66 Karmel. 1940, Astron. Nach. 269. Hl 
6650AB 5.62) 6:01 6080 5250 0.42 —0.37 Makemson. 1933, Astron. J. 42. 
0.34 —0.47 Gasteyer. 1954, Astron. J. 59. 
7307 6.62 6.85 9700 9700 Dod 0.74 Giintzel-Lingner. 1952, Astron. Nach. 281. 
1.97 0.29 Arend. 1954, Bull. Astron. Belg. TV. 
7704 7.26 7.42 9000 7700 0.78 —0.11 Wierzbinski. 1955, J.A.U. Circ. 5. 
7724 2.51 3.84 4300 5060 1.5xX10-£ —3.82 Giintzel-Lingner. 1956, 7.A.U. Circ. 8. 
8032 9.36 9.89 5800 4100 4.35 0.64 Baize, 1948, J. Obs. 31. 
259 0.41 Baize. 1954, [.4.U. Cire. 4. 
8119 4.49 4.96 6580 5130 1.98 0.29 Rakowiecki. 1938, Wiadom. Matem. Warsaw 46. 
1.84 0.26 Rakowiecki. 1938, Wiadom. Matem. Warsaw 46. 
8630 Ji Oo ee 5000 5590 0.10 —1.00 Strand. 1937, Leiden. Ann. XVIII. 
0.10 —1.00 Wolf. 1949, Astron. Nach. 277. 
0.10 —1.00 Wolf. 1949, Astron. Nach. 277. 
8695 Sel Seah ke, 3950 7860 1.2X10°% —2.92 Schmeidler. 1939, Astron. Nach. 268. 
1.4K10-3 —2.85 Rosino, 1944, Mem. Soc. Astron. Ital. 15. 
8949 7.60. 8.47 3850 3950 0.09 —1.04 Luplau-Janssen. 1931, Astron. Nach. 241. 
0.13 —0.89 Florsch. 1954, 7.A.U. Cire. 3. 
eal —0.96 Heintz. 1955, [.A.U. Cire. 7. 
8974 5.04 7.50 7000 3750 0.017 —1.77 Wierzbinski, 1956, 7.A.U. Cire. 8. 
9031 (pert SEO 3550 3120 0.043 —1.37 Rabe. 1927, Astron. Nach. 231. 
9182 thes) bees) 6300 6700 2.39 0.38 Von Bezold. 1938, Astron. Nach. 267. 
2.49 0.39 Hopmann. 1945, Verdff. Univ. Sternw. Leipzig 8. 
9343 4.60 4.80 11600 12200 0.13 —0.89 van den Bos. 1937, Union Obs. Circ. No. 98. 
0.13 —0.89 Wierzbinski. 1954, J.A.U. Circ..4. 
9380 T0353) ieee: 5710 4140 Ond2 —0.92 Pensado. 1954, Univ. Madrid Pub. 13. 
0.13 —0.89 Wierzbinski. 1956, J.A.U. Circ. 8, 
9413 4.75 6.95 5630 3400 0.39 —0.41 Strand. 1937, Leiden Ann. XVIII. 
0.40 —0.40 Kiimmritz. 1956, 7.A.U. Circ. 10. 
9578 D2 aixoe 6060 6060 1.2t 0.08 Danjon. 1938, J. Obs. 21. 
1.38 0.14 Danjon. 1938, J. Obs. 21. 
1.10 0.04 = Aller. 1939, Astron. Nach. 268. 
0.93 —0.03 Muller. 1952, Bull. Astron. 16. 
‘ 0.93 —0.03 Van Biesbroeck. 1955, Yerkes Publs. 8. 
9617 5.66 6.03 6160 6300 1.42 0.15 Silbernagel. 1927, Astron. Nach. 234. 
1.12 0.05 Danjon. 1938, Bull. Astron. 11. 
9626BC 7.03 7.66 4930 6300 0.32 —0.49 Baize. 1952, J. Obs. 35. 
9716 7.41, 7-49 5120- 5080 1.50 0.17 Pitman. 1934, Publs. Astron. Soc. Pacific 46. 
1.86 0.27 Stephens. 1939, Publs. A. A. S. 9. 
9979 5.78 6.84 8000 6120 3.40 0.53 Strand, 1937, Leiden Ann. XVIII. 
3.09 0.49 Giannuzzi. 1951, Contrib. Sct. Oss. M. Mario 164, 
3.42 0.53 Rabe. 1954, 7.A.U. Circ. 4. 
10075 EES 4950 5000 Sil 0.51 Finsen. 1933, Union Obs. Circ. No. 90. © 
2.90 0.46 Siegrist. 1954, Univ. Madrid Publ. 11. 
10157 3.05 5:60 6140 5080 0.26 —0.58 Silbernagel. 1928, Astron. Nach. 233. 
0.26 —0.58 Baize. 1949, J. Obs. 32. 
10235 6.84 8.70 6200 4140 0.14 —0.85 Rabe. 1926, Astron. Nach. 231. 
11046 4.12 5.93 6500 4000 2.24 0.35 Strand. 1937, Leiden Ann XVIII. 
2.24 0.35 Geffers. 1954, Veréff. Bonn 39. 
11186 8.33 8.63 9550 10080 4.56 0.66  Luyten. 1934, Astron. J. 43. 
Salle) 0.57 Luyten. 1934, Astron. J. 43. 
3.65 0.56 Wilson. 1935, Astron. J. 44. 
11635AB Sle One, 8800° 10 540 0.18 —0.74 Giintzel-Lingner. 1956, [.A.U. Cire. 8. 
11635CD Beil sie 5ul 10320 10620 ike ZAb 0.08 Giintzel-Lingner. 1956, J.A.U. Cire. 8. 
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TABLE II (continued) 


References 


Int. pv 

a Jame my, Mes Ti T2 p logp 
871 Sl Tear 5700 =4150° 0.63 —0.20 
| 0.70 —0.15 
447 Geode Sco! 6280 5080 0.62 —0.21 
_ B89 8.63 8.84 5320 5480 11.4 1.05 
13.0 Det ts 
12.6 1.10 
13.6 Silane) 
12.4 1.09 
169 9.70 10.24 6310 6280 4,33 0.63 
3.28 0.51 
296 Avis) 6011 10360 11710 0.038 —1.42 
360 6.41 7.24 6390 6630 0.21 —0.68 
499 5.89 6.40 6340 ~ 5050 0.073 —1.14 
636 5.41 6.18 3500 3300 0.66 —0.18 
0.66 —0.18 
-}270 4.81 6.14 6600 6000 0.66 —0.18 
(971 4.42 4.63 7000 6050 0.18 —0.75 
0.15 —0.82 
| 0.10 —1.00 
5988 5.90 7.29 9080 5890 0.47 —0.33 
5428 6.08 7.43 7400 4150 0.020 —1.69 
5665 8.48 9.12 4970 3900 0.46 —0.34 
‘7149 6.36 4.49 5810 6020 0.16 —0.80 


Rechtenstamm., 1940, Astron. Nach. 268. 
Baize. 1943, J. Obs. 26. 

Finsen. 1937, Union Obs. Circ. No. 98. 
Strand. 1937, Bull. Astron. Inst. Neth. 8. 
Baize. 1945, l’Astron. 59. 

Rabe. 1948, Astron. Nach. 276. 

Cester. 1950, Oss. Trieste P. 227. 

Baize. 1954, 7.A.U. Cire. 4. 

Wierzbinski. 1951, Acta Astron. 4. 

Baize. 1953, J. Obs. 36. 

Rabe. 1948, Astron. Nach. 276. 

Baize. 1943, J. Obs. 26. 

van den Bos. 1933, Union Obs. Circ. No. 90. 
Flechter. 1933, Monthly Notices Roy. Astron. Soc. 92. 
Zagar. 1934, Atti R. Inst. Veneto di Sc. 
Komendantoff. 1936, Astron. Nach. 260. 
Strand. 1941, Astron. J. 49. 

Giannuzzi. 1954, Contrib. Sci. Oss. M. Mario 204. 
Rabe. 1954, J.A.U. Circ. 2. 

Van Biesbroeck. 1937, Yerkes Publs. 5. 
Giintzel-Lingner. 1956, /.A.U. Circ. 9. 
Tannenbaum. 1938, Union Obs. Circ. No. 99. 
Franz. 1955, 1.A.U. Circ. 6. 


There are several sources of uncertainty in the 
lerived apparent photovisual magnitudes. First is the 
verage probable error +0.04 in the combined mag- 
jitude taken from Harvard Photometry (1938). 
lecondly, errors affecting the individual magnitudes m 
nd m2, according to (8), depend on the amount of Am 
nd are largest when Am=0.00. In the worst case 
\im=(.00-+0.10 (0.10 being the highest average p.e. 
f the photometric catalogues examined), whence 
'=0.75+0.06 and the resulting errors are m+0.07, 
ro 0.12. 

Another source of error is introduced with the 
olometric corrections. According to Kuiper, this should 
ot exceed about 0™10 p.e., whence the total probable 
tror of a bolometric magnitude is not more than 
-0.15 (0.12 for photovisual magnitudes and +£0.10 
or bolometric corrections). Since this error has been 
omputed in the most unfavorable case, we can assume 
hat on the average the p.e. of a bolometric magnitude 
; sensibly smaller and about +£0.10. In formula (7), 
utting 7;=T2, mi=m2, and taking the logarithmic 
erivative, we find 


dp/p= +0.6 dm/Mod= 1.38 dm. 


fonsequently as a result of an incorrect knowledge of 
lagnitudes, p should be affected on the average by 
e.+0.149 and in any case no more than +0.20p. 
‘his considers the accidental errors in m only. But 
losterhoff (1945) has found evidence of considerable 
rrors in the zero point for different areas. A good 
stimate of the systematic p.e. in m may be +0.13. 
his affects p by da=+0.18p p.e. The zero point for 
olometric magnitudes is that value for which Mo 
--+4.73. This value has been deduced from the sun’s 
hotovisual magnitude (Stebbins and Kron 1957), 
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Fig. 1. Comparison of visual (full line) with eclipsing (dotted 
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line) binaries in percentage distribution of logp. 
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Taste III. Eclipsing variables. 
Name logp Name logp Name logp 
Detached systems 
V 805 Aql —0.69 o Aql —0.95 TT Aur —0.76 
WW Aur —0.57 AR Aur —0.39 8B Aur —0.76 
SZ Cam —1.66 GL Car —0.95 YZ Cas —0.69 
AR Cas —1.38 SZ Cen —0.62 AH Cep —1.12 
a CrB —0.88 AI Cru —1.09 Y Cyg —1.07 
V 477 Cyg —0.10 YY Gem 1.43 RX Her —0.60 
TX Her —0.36 UV Leo —0.13 U Oph —0.83 
V 451 Oph —0.69 VV Ori , 1.03 AG Per —0.71 
& Phe —0.73 RS Sgr —1.33 WY Ser, —0.65 ; 
Y Cam —1.30 S Cnc —1.71 RS CVn —1.56 \ 
TV Cas —1.08 Z Her —1.00 AR Lac —1.01 
Semi-detached systems ~ 
RT And —0.19 TW And —1.37 ZX And —1.02 
R CMa —0.57 CV Car —2.51 RZ Cas —0.61 
U Cep —1.26 U CrB —1.39 W Del —1.53 
Z Dra —0.65 TW Dra —1.34 RW Gem —1.40 
SZ Her —0.55 u Her —1,21 RX Hya —0.97 
T LMi —1.26 6 Lib —1.36 RW Mon —1.04 
TY Peg —1.18 RT Per —0.41 RY Per —1.77 
B Per —1.21 V Pup —1.13 U Sge —1.42 
V 356 Sgr —2.21 V 505 Sgr —0.81 uw’ Sco —1.04 
RW Tau —1.02 d Tau —1.62 X Tri —0.63 
TX UMa —1.25 W UMi —1.28 Z Vul —1.31 
RS Vul —1.63 
Contact systems 

AB And 0.03 S Ant —0.55 i Boo 0.19 
TX Cnc —0.11 VW Cep 0.12 TW Cet 0.08 
RZ Com —0.01 YY Eri 0.07 SW Lac 0.07 
V 502 Oph —0.28 ER Ori —0.18 U Peg —0.06 
RZ Tau —0.17 W UMa 0.01 AH Vir —0.18 


mo=—26.7340.03. We think the zero point for 
bolometric magnitudes should change no more than 
+0.10. Accordingly, from this cause the densities 
might vary by an amount +1.33dMop=+0.13p. 

This discussion shows that the combined effects of 
different kinds of errors in the magnitudes are as 
serious as those of the errors in T. 

We have neglected interstellar absorption since these 
stars lie within about 100 parsecs. Also, the selective 
absorption producing stellar reddening leads to esti- 
mate temperatures a little lower than the true value, 
and this effect should tend to compensate a diminution 
of m in (7). 

To obtain the orbital elements P and a, we have 
consulted Baize, Bull. Astron. 10, 273, 1937 oad J. Obs. 
39;11;).1950, and for the orbits later computed, the 
annual reports of the Monthly Notices Roy. Astron. Soc. 
and the Circulaire d’Information of I.A.U. Commission 
26. 

We have not considered the error in p due to an 
incorrect knowledge of orbital terms, since these affect 
p very slightly, as one can see in Table I, where the 
values of density are given for each binary and for 
different orbits. 


CONCLUSIONS 


In Table II, we have collected the results for 52 
double stars having satisfactory data. The columns 


magnitudes m, and mz on the International scale) 
(3) effective temperatures JT; and T2; and (4) mea| 
density in units of solar density. | 

For comparison, eclipsing binaries are considered fo| 
which the mean density is known by means of phote 
metric and spectroscopic observations only. Adoptin, 
the data of the Kopal and Shapley catalogue (1956) 
we have calculated logp for 82 systems given in Tabl: 
Iil. ji 

Figure 1 shows the resulting percentage distribution 
with log of the visual binaries (full line) and of th 
eclipsing binaries (dotted line). For those visua 
binaries having more than one set of orbital data anc 
hence more than one determination of p, we have, 0 
course, counted each system only once, using the mor 
recent orbit in those three cases where the logp value 
were not included in the same interval. 

It is clear from the figure that the visual double star 
of this selection are on the average eight times dense 
than photometric binaries, also that their mean densit; 
is normally little different from the sun’s density 
However most eclipsing binary stars of Kopal’s cata 
logue have spectral types A and B, while the visua 
doubles which we have considered are predominentl: 
of spectral type later than dFO. 

Unfortunately, not enough systems are available t 
compare between the two groups of binaries the dis 


»sion in density among stars of the same spectral 
1s. But we would in fact expect that B and A stars 
Juld have about 0.1 the mean density of G and K 
ss (see, e.g., Ledoux and Walraven 1958). Thus it is 
jjte probable that the differences between the histo- 
as of star density for the present sample of the two 
hes of binary systems are due entirely to the result of 
) ervational selection. 


j 
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The computation is based on 243 measures of position covering the period from 141 days before to 356 
days after perihelion. Planetary perturbations were taken into account but remained small. The hyperbolic 
excess comes out 48 times its probable error. The “original” orbit, however, turns out to be a very large ee 


ellipse. 


OMET 1947 k& (1948 I) was discovered on Sep- 
tember 24 by M. J. Bester at the Boyden Station, 
Bloemfontein, South Africa on a sky-patrol plate taken 
with a 3-inch Ross camera. It appeared as a small, 
diffuse, centrally condensed object of magnitude 11 
moving southwestward in Eridanus. 

The first accurate positions were obtained September 
28 at the Lick Observatory. On October 9 the mag- 
nitude was estimated as 9.8 and a broad tail was noted 
in 350° which extended on October 14 to a distance of 
10’ from the well-defined nucleus. By the end of 
October the comet was out of reach of northern ob- 
servers and went as far south as —55° on November 22. 
In December the comet moved far in the western sky, 
Simonov calling the magnitude 9 on December 22. The 
last observation before conjunction with the sun was 
by Johnson at Johannesburg on 1948 January 12, when 
he estimated the brightness as 7th mag. and noted a 
tail of one degree in length. 

After perihelion passage on 1948 February 16, the 
comet was first observed in the morning sky by Hirose 
on March 3, when he called the magnitude 7. By the 
middle of March it could be seen with the naked eye. 
It appeared then as large round coma centrally con- 
densed with neither nucleus nor tail. Being conveniently 
located, the comet was abundantly observed in March 
and April, when the declination went as far north as 
73° (April 23). A minimum distance of 0.74 a.u. from 
the earth was reached on April 6, at which time the 
daily geocentric motion was nearly 3°. In May the 
brightness decreased rapidly. On May 28, although the 
magnitude had dropped to 11, a tail 20’ long was again 
visible and there was a sharp nucleus in a coma 1’ in 
diameter. Observations now become scarce. A last 
evening measure was obtained at low altitude by Van 
Biesbroeck on August 1. After the comet reached con- 
junction with the sun, Jeffers secured a pair of plates 
on 1948 November 11, when he estimated the mag- 
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“at 17.5. The measures therefore cover an arc of 
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nitude of the nearly stellar object as 17.7. On 4 
February 6, Jeffers secured the last observations w, 
the Crossley reflector, the magnitude being estima] 


from 141 days before to 356 days after perihelion. Fr} 
short arcs, preliminary parabolic elements w) 
deduced by Cunningham and Maxwell. A more accur} 
orbit, given in the third column of Table I, was obtair} 
by Cunningham from measures extending to Octobei 

The third orbit was used to establish an accurj 
ephemeris for comparison with all the measures. Ta’ 
II gives the list of residuals of the individual obs 
vations. The dates are corrected for aberration tii 
and reduced from UT to ET. Published parallax facté 
were corrected when necessary. The references to t 
original measures are found in the Astronomisch 
Jahresbericht 1947-49. A plot of the residuals agair 
time was useful in illustrating observations showil 
abnormal discrepancies and also the scattering of t] 
measures by the individual observers. By correspon 
ence with the observers most discrepancies . we 
eliminated. A few outstanding ones could not be a 
counted for. These are put in parentheses i in the lis 
meaning that they were not included in making ti 
normal places. Table III gives the list of the observe 
and the weights from 3 to 1 generally attributed to the 
from the inspection of the scattering of the measur 
on the plot mentioned above. One numerous series | 
measures, namely that by Chis (1949) at Ch 
(Rumania), showed such wide scattering that it w 
thought better to omit it altogether rather than | 
include it with low weight. The graphica! plot of tl 
residuals showed the large uncertainty of the position 
A well-defined nucleus was mentioned in the ear. 
measures in October 1947; however, when the com 
reappeared after perihelion it had a large coma 8’ : 
10’ in diameter (Krumpholz, March 15 and April 1 
centrally condensed but without nucleus. Besides, tl 


TABLE I. Orbits computed. 


Cunningham (1947a) 
T=1948 Feb. 17.727 UT 


Maxwell (1947) 
1948 fae 17.195 UT 
349°2 


Cunningham (1947b) 
1948 Feb. 16.43302 UT 


w= 348°50’) 350° 24100 
2=270 2 '1947.0 270 1G 1947.0 270 . 69135 }1947.0 
i=140 12 140 21 140. 56326 
g=0.7677 0.7596 0.7478089 
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f) 10.92977 
14.30930 
14.89298 


17.87679 
18.39304 
19.04058 
19.05202 
§ 19.47478 
 20.45962 
) 21.29532 
2.06005 
~ 23.09170 
23.34160 
23.79649 


} 24.01380 


W., 27.81526 
. 29.79989 
ty. 1.02085 
5.02748 
5.78934 
8.05261 
9.80221 
12.02740 
13.01841 
13.72427 
15.04078 
16.05354 
17.78984 
18 .04317 
22.02018 
| 23.73518 
Ww. 27.72754 


2.08185 
2.73408 
5.10573 
6.74373 
13.07210 
16.76741 
 18.03330 
te. 22.46467 


1948 
4.75577 


n. 12.75374 


tb. 25.10846 
. 3.80290 


3.80915" 


3.81401 
5.45078 
5.81267 
5.82778 
7.20737 
7.20983 
7.21136 
8.79362 
8.79779 
8.80126 
9.12669 
9.84436 
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TABLE II. Residuals. 


On 


oe (ar ae ee pe 


atest (eli: eedieny 


TA TT) TP AE LSI Sica ee ee Sc ot bal a af te 


Aa A6 

0804 — O%4 
0.01 — 0.8 
0.02 + 0.1 
0.03 + 2.0 
0.04 +. 0.5 
0.33 a hl 
0.01 + 0.8 
0.06 aaa eA: 
0.25 ap ale® 
0.06 + 3.2 
0.23 + 3.7 
O525 5.9 
0.54 ap was 
0.51 + 1.2 
0.13 + 6.2 
0.02 On) 
0.14 + 5.3 
0.17 + 2.3 
0.82 corona 
0.99 + 5§.2 
0.44 + 4.0 
0.24 + 4.3 
0.57 + 6.8 
0.63 + 9.1 
0.60 (— 5.3) 
0.30 + 12.3 
0.56 + 8.4 
0.37 + 8.1 
0.21 Sie whet) 
0.65 + 8.6 
1.85 + 7.6 
0.77 + 12.9 
0.40 + 3.4 
0.97 + 15.8 
1.01 + 13.7 
0.43 + 4.5 
1.43 “EY s2 
0.89 + 11.3 
1.83 + 9.0 
0.70 + 9.5 
0.52 a 9.0 
0.70 + 10.6 
0.19 + 11.8 
1225 + 8.5 
1.05 + 15.0 
1.14 + 13.4 
0.76 + 13.0 
111 + 14.1 
0.62 + 14.3 
1.03 + 17.2 
0.26 + 5.9 
0.59 — 15.8 
0.44 — 17.9 
0.36 — 17.6 
0.66 — 25.2 
0.95 — 25.1 
0.90 — 29.7 
0.15 — 34.6 
0.21 — 27.0 
0.80 SURV) 
0.68 — 40.3 
0.64 — 51.6 
0.87 — 39.5 
0.64 — 44.8 
0.62 — 50.0 


Wt Obs. ET Aa Ad Wt Obs 
Mar. 10.20181 — 0:90 — 44°6 3 PB 

3 a Mar. 11.52330 SA EAGT 3 I 

5 IL 11.52606 8134 — 49.0 3 I 
3 IL 12.14322 Bri 31 Ege 3 HP 
3 RW 12.19713 eG 5A. 3 a PL 
5 ni 1219892 Sis LETTS 3 PL 
> I 1220093 — 0.96 meee 3 PB 
3 VB 12.20258 Reis 50.5 3 BA 
5 nl 14. 12826 = 1.68 62.8 3 HP 
3 VB 14.18471 es — 64.8 3 PL 
5 BJ 14.18615 2 its0 265.2 3 PL 
14.21363 isk E1638 3 PB 
2 By 14.76172 S76 — 65.5 3 HT 
3 Si. 14.76589 er Le50 — 68.4 3 HT 
1 D 14. 76936 — 1.60 60.3 3 HT 
3 SA 15.13970 — 1.66 S68 3 HS 
3 G 15. 14892 — 1.88 Sen 1 SI 
3 G 15.14995 S78 705 3 HP 
3 VB 15.17238 S137 — 69.2 3 KV 

i D 16.10786 Sri ga erg Z O 
1 D 16.18592 a11.62 ance 3 BA 
3 VB 17.20419 Se iv6s eS eese7 3 PB 
By 18.79855 24.98 eS0729 2 Mi 
1 D 19.17668 = 2.08 2580755 3 BA 
21.20742 953 —108.4 3 PB 

2 BJ Mar. 22.14241 e570 —122.8 1 SI 
4 1 Mar. 24.14690 — 2.98 ~142.5 3 HS 
5 24.14749 =53.80 —142.6 3 PT 
3 a 24. 16190 22318 138.7 3 KV 
i 4.17114 335 —140.1 3 PB 
= i 2508619 = 3107 —146.3 3 PL 
i hs 25 .09013 = 2.95 aeicat 3 PL 
; ic 25.11624 2.66 wisss 3 HS 
: i 25. 15964 = 53.04 151.9 3 PT 
‘ 26.11205 32 ~162.1 3 HS 
; o 2707922 233.39 ary 3 PL 
; i 27 08274 3644 ~172.9 3 PL 
; 5 27.14615 en —169.4 3 KV 
i a 27. 15925 45.04 2175.0 3 HS 
5 i 2716684 T3538 21720 3 PT 
i i 2809393 3156 24531 3 PL 
2809746 3.65 —182.6 3 PL 
; 28. 10661 Sr 3 — 184.0 3 HS 
2 a 29. 12652 2253.00 ios 4 3 HS 
i 29. 14409 es —194.2 3 KV 
z i 29.81113 — 4.06 —203.0 3 HT 

‘ 5 29.97528 £2393 —206.7 2 O 
; a 29 98440 S470) —204.8 3 HP 
Ns 3003608 ee sed —204.6 3 HS 
3 ce 30.72990 =A —210.4 2 Mi 
31 .03372 235495 —214.4 3 Ma 
31.03574 — 4.07 ~214.7 3 HP 
3106988 — 4.00 2319.7 a HS 

1 J Apr. 1.05564 — 4.19 —~ 228.6 2 O 
1.14318 Ta 226.4 3 KV 

1 JJ Apr. 1.98891 — 4.02 —236.6 2 Vo 
2 J Apr. 2.96208 25207 —247.7 2 KK 
3 HT 2.96258 AS —~248.3 3 HP 
3 HT 3.01784 — 4.18 —247.3 3 Ma 
3 HT 4.90053 are —267.3 3 Ma 
3 VB 4.94469 = 4.37 —~260.8 3 HP 
3 HT 4.97687 = AIG ~257.0 2 KK 
3 HT 5.69663 2549 ~269.7 2 Mi 

3 PB 5.84196 S80 273.6 2 O 
3 PL 5.93365 — 4.48 —262.5 3 HP 
3 PL 5.97068 aay «3 —268.4 2 Ma 
3 HT 6.07096 en —267.9 3 PL 
3 HT 6.07483 — 4.00 —270.9 3 PL 
3 HT 7.01540 3% $9 —281.5 3 Ma 
3 AA 7.11493 os as —284.5 3 BA 
2 Mi 7.97112 — 3.66 —286.9 3 HP 
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TABLE II (continued) 
ET Aa A6 Wt Obs. ET Aa A6 Wt 
Apr. 8.05960 — 4:07 —283"4 3 KV May 1.89865 +4115 — 48%5 3 
8.61894 = OD == 288 3 HT 2.92925 +39.90 = 54058 3 
8.62310 =216 —289.4 3 HT 3.83143 (+33.44) 31.9 3 
8.62538 57 = 2887 3 HT 3.92392 +38.08 = SO 3 
8.62866 Ne oe —285.4 3 HT 3.96547 +37.43 — 39.4 3 
8.71561 Bye —293.5 2 Mi 5.00353 +36.40 ee 3 
8.91718 a5 52) = 772 3 KV 5.02993 +35.94 Sie 3 
9.69537 — 2.46 2703.8 3 HT 5.67480 +34.17 pe oe 3 
9.85860 = OTT —290.0 p Vs 5.68244 +34.31 aie 3 
9.94282 =. 78 —295.8 3 HP 5.83658 +35.55 =e 3 
11.01087 — 0.69 —298.4 3 Ma 5.97035 +35.90 = ORs 3 
11.01101 — 0.44 Bey ae 3 HS 6.94571 +33.83 = 35.0 3 
11.69330 4+ 1.55 = 801.8 2 Mi 6.95160 —+33.80 = 25h 3 
11.90552 + 0.50 312.2 2 Vs 7.89433 “433.01 = ae 3 
11.99996 + 0.96 —296.6 3 Ma 7.90022 +31.40 =" 49.9 3 
12.08619 41250 = S02 1 SI 7.96166 +32.53 = Os 1 
12.11704 1:20 —300.4 3 BA 8.82807 +31.76 Be 3 
12.87244 4D: Gl —284.6 gram e | 8.84311 +29.94 a SN | 2 
Apr. 12.96975 + 3.36 —307.4 3 PB 8.87626 +31.04 st 3 
8.94136 239.25 OF 3 
Apr. 13.93649 + 5.89 —295.2 3 HP 9.80057 +30.75 — 18.6 3 
13.95434 + 4.13 — 304.2 2 VS 10. 71341 +2994 = line 2 
13.99226 9 BF OIF 1 SI 11.86459 +30.56 =) ida 2 
14.92505 + 9.13 —292.9 3 PB May 12.85468 +29.50 lice 2 
15.04658 4+ 9.47 —292.0 3 HP 
15.88555 +12.43 —291.6 3 HS May 16.93383 +24.07 10.9 2 
15.95727 +14.39 —298.9 2 Vs 17.97231 +24.13 = 4148 2 
17.04995 421.23 —278.0 3 HP 19.87922 +23.01 ASS 3 
17.90949 1:95 58 —286.1 3 HS 25.87103 +21.03 = 14.3 3 
17.93095 +23.40 —293.7 3 Ma 2588940 +20.29 — 15.5 3 
18.04654 +28.16 187.8 2 Vo 27 89888 +19.86 ede 3 
19.12045 430.93 —247.4 3 AU 2791913 +19.88 = 4g 3 
1990276 +34.77 —241.6 3 KV 28. 15697 +19.91 =) 154 3 
20.02321 +37.03 = 940.9 2 KK 30. 86807 +19.35 Bee ie 3 
21.82094 +42.91 —199.3 3 KV 31.82558 +19.40 = 27 3 
Apr. 21.83647 445.18 (—263.7) 2 Vo 31.86325 +18.39 = 10 3 
June 1.90298 +19.10 (+ 16.8) y. 
Apr. 22.82885 +46.60 —182.0 3 KV 2.50127 +18.18 — 19.6 2 
22.98851 449.57 —176.9 3 HS 3.88078 +16.61 =A OF 1 
2383466 451.54 —192.5 3 AA 4.87060 +17.68 (= $8.0) 1 
24.78009 +50.09 + 147,5 2 Vs 6.56530 +17.49 pri) 2 
24. 80347 +52.14 —136.6 2 Vs 9.16191 A725 21.0 3 
24. 82994 +39.06 1378 3 KV June 11.88322 416.47 = 29% 3 
Apr. 24.98660 +48 .65 Sees 2 BL 
y June 24.86164 +14.89 = WO 3 
Apr. 27.94068 +48 .04 ee 3 AU July 3.17988 +14.23 Seyi 3 
2795193 +48.11 S576 3 AU -15.18710 +13.49 32.6 3 
28.61441 446.47 = 79.3 2 Mi 15.19474 +13.46 re | 3 
28. 87467 +4745 SEIS ne SI 24.12121 +12.69 — 36.6 3 
28.93579 +4717 — 69.9 3 HP Aug. 1.09880 12.73 = NIONS 3 
poe +45.42 — 66.6 3 AA 
- 79020 144,84 — 67.8 3 Ma Nov. 11.46974 +13.97 — 60.8 3 
29 .81504 +43 .42 — 70.6 3 KV “41. ; cae 
99183384 duslog |i. = ology tn es amen eg i 0 ee 
30. 83680 +47.02 obi: 1 KK 1 
30. 88221 4.42.35 2.750 3 Ma i 
30.93365 +4266 — 54.4 3 PB Feb. 6.42341 +23.08 — 46.0 3 
May 1.82438 +46.43 — 49.4 3 Ma Feb. 6.45258 


fast motion of the comet at that time and its circum- 
polar position added to the difficulty of the measures, 
which show a surprisingly large scatter for such a 
bright object. In May 1948 there appeared a better 
defined nucleus, but this advantage was soon lost by 
the rapid decrease in brightness. 

The blank lines in Table II indicate how the 243 
observations have been grouped in 17 normal places. 
For each group the mean was obtained by a least- 
squares solution including a term in ?. Table IV gives 


the list of normal places and their weights in con 
2, 3, and 4. 

The noninclusion of a few bracketed large residue 
made slight differences in the weights of the two co 
dinates. It was thought best to disregard those mine 
differences. 

Next the perturbations by all planets except Mere 
and Pluto were computed at 20-day intervals by Encke 
method, using 1948 February 9.0 as osculation dat 
The comet did not come near any of the major planet 


f minimum distance to Jupiter was 4.5 a.u. on 1948 
ch 9. Minimum distance to Venus was 0.7 a.u. 
9/948 January 8 and to the earth 0.7 a.u. on April 6, 
| these approaches were of short duration. Hence 
4) mall corrections for perturbation listed in Table IV 
plumns 5 and 6. From these we find in columns 7 
8 the final residuals from Cunningham’s second 
bolic orbit. 
was though that these residuals were small enough 
} se the same parabolic orbit in computing the dif- 
# tial coefficients of the equations. For that purpose 
tcke’s (1929; p. 330) form was used. Table V lists 
» 34 equations of condition. 
or the least-squares solution of these equations I 
= indebted to the National Bureau of Standards at 
Ider, Colorado. Through the good offices of Dr. F. E. 
© ch and W. B. Jones the work was done on an IBM 
{ computer “‘in less than a minute [sic ].” The cor- 
ions to the elements and their mean errors are 
nin Table VI. It will be seen that the two elements 
nd were mainly responsible for the deviation of 
| representation. 

able VII shows the final representation of the normal 
tes. The mean error of weight 1 comes out as large 
1794, which corresponds to 1712 for a good obser- 
| jon of weight 3. This is not surprising in view of the 
‘sical appearance of the comet when the large 
. jority of the measures were obtained. The hyperbolic 
sess of the eccentricity is small but well established 
ize it is 48 times as large as its mean error. 

it is of interest to find out what the “original” 
‘nents are and in particular what the initial value of 

eccentricity was before the comet entered into the 
'er part of the planetary system. 
|| was fortunate in securing the interest of B. G. 
irsden at the Yale Observatory in this question. 
uding the perturbations by the earth, Jupiter, and 
turn, he wrote a program and computed the back- 
rd Bad forward values of e and 1/a for an interval of 
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TABLE III. Observers. 


~~ 


Abbreviations Observer 
AA Adamopoulos 
AU Arend 
BA Boyer 
BC Bobone 
BJ Bruwer 
BL Brown 
D Dujinin 
G Giclas 
HP Hurnik 
HS Hoffmeister 
HE Hirose 
I Inglis 
JJ Johnson 
JL Jeffers 
KK Konopleva 
KV Krumpholz 
Ma Martynov 
Mi Mitani 
O Oterma 
PB Polit 
ID, Pels 
12M Prétre 
RU Rigaux 
RW Reuning 
SA Schmidt 
SC Simonov 
SI Schneider 
SL Shane 
VB Van Biesbroeck 
Vo Vorsat 
Vs Vsekhsvyatsky 


Location 


Weight 


Athens 

Uccle 

Algiers 
Cordoba 
Johannesburg 
London 
Santiago 
Flagstaff 
Poznan 
Sonneberg 
Tokyo 

Lick 
Johannesburg 
Lick 

Kiev 

Vienna 
Kazan 

Kyoto 

Turku 


‘Barcelona 


Leiden 
Toulouse 
Uccle 
Washington 
Algiers 
Canberra 
Innsbriick 
Lick 
Yerkes 
Kiev 

Kiev 
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12.5 years on an IBM 650 computer 
end.” The values, referred to the center of gravity of 


the solar system, came out as follows: 


“over the week- 


Elements € 1/a 
Osculating (1948 Feb. 9) 1.0003071 —0.0004105 
+0.0000064 +0.0000085 
Original (1935 July 17) 0.999935 +0.000087 
Future (1960 Sept. 23) 1.000237 —0.000317 


Once more we find that the apparent hyperbolic 
excess of the eccentricity changes into an original elliptic 
one corresponding to an orbit of very large size and 


TaBLeE IV. Normal places. 


) O-C Correction for perturb. Corrected O-C 
| ET Aa cosé Ad Wt Aa cosd Aé Aa cosd Aé 
| = 
| 1947 Oct. 6.0 -— 172 + 174 29 +1°6 179 + 074 — O%5 
21.0 — 3.9 + 4.8 26 +1.6 —1.9 — 2.3 i eae 
Nov. 12.0 — 6.9 + 8.8 25 1.2 —2.2 == taal + 6.6 
Dec, 11.0 — 5.6 + 11.9 14 +0.1 —0.7 — §.5 + 11.2 
1948 Jan. 4.756 — 8.2 + 14.2 1 0.0 0.0 — 8.2 + 14.2 
12.754 — 12.4 + 17.2 1 0.0 0.0 — 12.4 + 17.2 
Mar. 6.0 — 8.3 — 29.5 46 0.0 +0.1 — 8.3 — 29.4 
16.0 — 25.5 — 78.8 69 —0.3 0.0 — 25.8 — 78.8 
29.0 — 55.7 —194.1 86 —0.7 —0.1 — 56.4 —194.2 
Apr. ) 7.0 — 40.4 —283.1 91 —1.0 —0.6 — 41.4 — 283.7 
17.0 + 43.4 —278.1 41 —1.4 —1.4 + 42.0 —279.5 
; 24.0 +-222.3 —153.6 18 —1.6 —0.7 +220.7 —154.3 
May 5.0 +-221.1 30.7 99 —1.8 +0.6 +219.3 — 30.1 
30.0 +175.2 — 17.6 45 —1.7 +1.8 +173.5 — 15.8 
July 16.0 . +153.6 — 31.5 18 —1.2 +2.4 +152.4 — 29.1 
1948 Noy. 11.485 +168 .4 — 60.0 6 oe +4.9 +170.6 — 55.1 
1949 Feb. 6.438 +219.2 — 46.2 6 -+-5.0 +6.3 +224.2 — 39.9 
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TABLE V. Equations of condition. 


Date de dq dT dw di dQ O-C 
Ina a 
1947 Oct. 6.0 —0.57096 —0.74970 + 50.97 +1.18984 —0.40114 —0.81325=+4+ 0.4 4 
21.0 —0.50972 —0.97588 + 363.71 +1.32937 —0.49702 —0.95880 = — 2.3 | 
Nov. 12.0 —0.16330 —1.38149 +1744.07 +1.29548 —0.47941 —1.14836 = — 5.7 i 
Deets G +0.16610 —1.03334 +2779.61 +0.57989 —0.05715 —0.77049=— 5.5 | 
1948 Jan. 4.756 +0.11862 —0.70947 +2748 .85 +0.36716 +0.07224 —0.55178 = — 8.2. 4 
12.754 +0.09826 —0. 63584 +2788 .16 +0. 34666 +0.08823 —0.51457 = — 12.4 ¥# 
Mar. 6.0 —0.09029 —0.18717 +3090. 58 +0.38498 +0.02803 —0.30608=— 8.3 4 
16.0 —0.17347 —0.20919 +3201.44 +0.41617 —0.05568 —0.23318 = — 25.8 4 
29.0 —0.33987 —0.21592 +3532.95 \ +-0.50627 —0.26199 —0.13655 = — 56.4 |} 
Apr. 7.0 —0.43255 +0.03275 +3737 .81 +0. 68091 —0.48541 —0.20188 = — 41.4 i} 
17.0 —0.26830 +1.03042 +3648. 41 +1.17349 —0.73384 —0.78299 = + 42.0, 5 
24.0 +0.12738 +1.67348 +2413.65 +1.29535 .——0.58644 —1.23155 = +220.7') 
May 5.0 +0.39924 +1.40924 + 832.06 +0.90912 —0. 20964 —1.11187 = +219.3 | 
30.0 +0.44106 +0.94567 + 329.85 +0. 63067 —0.03929 —0.87097 = +173.5 io 
July 16.0 +0.42811 +0.70829 + 306.34 +0.55966 —0.06586 —0.77047 = +152.4 | 
Novy. 11.485 +0.45352 +0.67100 + 351.97 +0.64678 —0.26908 —0.89176 = +170.6 | 
1949 Feb. 6.438 +0.85692 +0.76173 + 260.87 +0.83778 —0.10521 —1.18787 = 4224.2 | 
In 6 | 
1947 Oct. 6.0 —0.29153 —1.11378 +1026.75 +0.22554 +0.34358 —1.06807 = — 0.5 i 
21.0 —0.12771 —1.27399 +1590. 28 +0.25422 +0.41068 —1.12555=+ 2.9 | 
Nov. 12.0 +0.17886 —0.85773 +1897 .81 —0.15223 +0.61199 —0.66287 = + 6.6 
Dec. 11.0 +0.04235 —0.01256 + 246.30 —0.45071 +0.71900 —0.01114 = + 11.2 
1948 Jan. 4.756 —0.02694 +0.02727 — 356.05 —0.20515 +0.58466 +0.03802 = + 14.2 
12.754 —0.02913 +0.02017 —_ 491.92 —0.13512 +0.53923 +0.05050 = + 17.2 | 
Mar. 6.0 +0.04998 —0.47797 — 2634.55 +0.07808 +0.06220 +0.42136 = — 29.4 © 
16.0 +0.09086 —0.90637 — 3621.52 —0.05336 —0.09556 +0.68819 = — 78.8 
29.0 +0. 10874 —1.85436 — 4987.96 —0.57792 —0.21102 +1.28750 = —194.2 | 
Apr. 7.0 —0.02337 —2.37974 — 4493.73 —1.04069 —0.09306 +1.62186 = —283.7 | 
17.0 —0.37053 —1.94818 —1366.11 —0.78711 —0.18857 +1.34227 = —279.5 © { } 
24.0 —0.44454 —0.80454 + 853.60 +0.06295 —0.64459 +0.57570 = —154.3 
May 5.0 —0.23004 +0.01631 + 968.41 +0.52697 —0.89051 +0.00497 = — 30.1 
30.0 —0.05586 +0.02888 + 189.29 +0.35981 — 0.84334 —0.01255 = — 15.8 
July 16.0 —0.00592 —0.08865 — 105.20 +0.12492 —0.74783 +0.08400 = — 29.1 
Nov. 11.485 —0.04175 —0.19841 — 155.79 —0.11082 —0.90071 +0.24159 = — 55.1 | 
1949 Feb. 6.438 —0.19362 —0.05620 + 27.99 —0.08203 —1.18614 -++0.11525 = — 39.9 ae 
TABLE VI. Corrections to elements and their mean errors. 4 } 
i] 
Initial elements (1950.0) Correction Mean errors | Final orbit (1950.0) 
e=1 +0.0003071 +0.0000064_ 1.0003071 y 
q=0. 7478089 a.u. +0.0003268 0.0000075 0. 7481357 a.u. q 
T=1948 Feb. 16.43334 ET —0.01014 0.000009 1948 Feb. 16.42320 ET 
= 270°73277 +-0°00742 ~ 0200069 270274019 EE: 
z= 140.56330 +0.00719 0.00057 140.57049 i 
w=350.24039 —0.02305 0.00088 350. 21734 { 
Taste VII. Final residuals. that the future perturbations will eject the comet fro 
the solar system. 
Date Aa cosé Aé Weights 
1947 Oct. 6.0 —1"3 —0"1 29 BRIGHTNESS OF THE COMET 
21.0 —2.8 Spoil 26 
ne hae nae ee a Numerous attempts have been made to represet 
1948 Jan. 4.756 49.5 44.0 1 the brightness of the comet as a function of its distant 
“A oe Saleh +6.5 1 from the sun. From photovisual observations coverin 
Mar. ee oe AeA a the period 1948 March 14 to May 6, Bacchus (1948 
29.0 906 ett 86 finds that the magnitudes reduced to unit distance fo 
Apr. ie se ce s 91 low 6.75+6.35 log r. This refers to the part of the oo 
24.0 5G _0'5 a included in a field of 18’ diameter, which therefore 
May 5.0 —3.1 1.5 99 cludes part of the tail. Ahnert (1948) has estimated 
Tay ae BER SAG 43 total magnitude visually during the interval March ‘ 
Nov. 11.485 40.9 a atG 6 to May 1 by extrafocal comparison with stars. He repr 
1949 Feb. 6.438 +3.0 —1.1 6 sents his numerous observations by 6.85+10.21 log 


and notes a marked bump in the light curve early i 
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TaBLe VIII. Comparison of results. 


1947 Oct. 4.09 


1947 Oct. 9.36 1948 Nov. 11.51 1949 Feb. 6.46 


: 
mtric distance 


2.375 2.306 3.972 4.889 
| tric distance 1.776 1.661 4.180 4.083 
ed magnitude Vis. 10"0 Vis. 928 Pay, ey Phi £75 
(Johnson) (Van Biesbroeck) (Jeffers) (Jeffers) 
‘ude reduced to unit distance from the earth 8.8 8.7 14.6 14.4 
‘is formula 9.1 9.1 10.6 iat il 
‘formula 10.7 10.6 1320 14.0 
» and Bertaud formula Seal 13.0 16.3 17.6 
 Baldet and Bertaud (1948) represent their photo- ACKNOWLEDGMENTS 


ic magnitudes covering the same period by 7.74 
4 logr. Bakharev (1949) computes the photo- 
¢ parameters representing the magnitudes ob- 
fi by five different observers. The constant term 
sfrom 6.29 to 6.76 and the coefficient of log 7 from 
) 10.5. He also expresses these magnitudes in the 
iative form suggested by Lewin. 

we compare what these formulas give for the 
st and latest observations (Table VIII) we see that 
pmet appears to have been brighter before perihelion 
after. 

) single formula represents the whole array of 
litude observations. 


In the laborious discussion of the observations I was 
assisted by W. Caldwell. Financial support from the 
National Science Foundation is hereby gratefully 
acknowledged. 
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The effects of solar radiation pressure on the motion of an artificial satellite are obtained. The results 
include the effects produced by the intermittant acceleration which results from the eclipsing of the satellite 
by the earth. Vectorial methods have been utilized to obtain the nonlinear equations describing the motion 
and the method of Kryloff-Bogoliuboff has been applied in their solution. 


INTRODUCTION 


EVERAL authors have considered the effects of 

solar radiation pressure on the motion of an artificial 
satellite (Parkinson, Jones, and Shapiro 1960; Musen, 
Bryant, and Bailie 1960; Musen 1960). With the 
advent of satellites of large area-to-mass ratio such as 
the Echo I balloon satellite, the magnitude of the 
radiation pressure effect is substantial and must be 
taken into account in the analysis of the tracking data. 
The central objective of this paper is to present the 
effects, including that of the earth’s shadow, which 
radiation pressure has on the motion of a satellite. 


KRYLOFF-BOGOLIUBOFF METHOD 


In an article entitled ‘‘Low acceleration takeoff from 
a satellite orbit,” Lass and Lorell (1961) have applied 
the method of Kryloff-Bogoliuboff to low thrust in the 
following cases: constant radial, constant transverse, 
and intermittant. The basic reference for the Kryloff- 
Bogoliuboff method is the book Introduction To Non- 
Linear Mechanics (Kiev, 1937), Chaps. 10-12, by 
N. M. Kryloff and N. N. Bogoliuboff. 

The method may be simply described as follows: 
Given a nonlinear equation 


dx/dt= pf (x, sint), 
where wv is a small number, one obtains upon integrating 
wt lr) —20) pe pte 
= i f(a, sinr)dr. 
Qa Qa J, 


Since pu is small, x may be considered constant during 
the integration from r=/ to r=f+2z. Since sint is 
periodic, the above equation becomes 


x (t+ 2) — x(t) 
Qn 


be Qa 
x if f(a, sinr)dr=pF (2). 
2m Yo 


Replacing the left side here by dx/dt, since the slope of 
the secant line is approximately that of the tangent, 
gives 


a ~ iG eon 


This result may be arrived at simply by averaging the 
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original equation over one cycle on the motion Wy 
respect to time. 


__CHANGE*IN ENERGY PRODUCED BY 
RADIATION PRESSURE 


Of fundamental importance in the examination! 
atmospheric drag on a satellite is the ability to sepia 
the energy loss (or gain) produced by radiation pr 
from the loss due to drag. The lifetime of the sa 
will be substantially affected by the changes in e 
resulting from radiation pressure. 

The effect of radiation pressure, neglecting the 
of the earth’s shadow, will produce only short- 
terms in the development for the semimajor axis of 
orbit (Musen 1960). However, the inclusion of © 
shadow results in quite a different phenomenon. 

Using the energy integral and Brown’s 6/dt oper. 
(see, for example, Musen 1954), we obtain 


da/dt=2a?v-F/p. 


r 


% 


The velocity vector v may be represented in terms | 
unit vector P directed from the center of the eartl 
perigee and a vector Q at right angles to P such t 
P, Q, R, with R the unit vector in the direction of 
Gace momentum vector, form a right-handed sy 
The vector F is the dace acceleration, @ i 
semimajor axis, and u=GWM, the gravitational const 


multiplied by the mass of the earth. We have # 
v={Q[ua(1—e?) ]} cosE—P (ua)! sinE}/r, a 
F=— Fo’, a 


where e is the eccentricity, E the eccentric anon. 
F the magnitude of disturbing acceleration, u? the t 
vector in direction of the sun, and r= cae cosE) 


+cosz sinw cosQ-+ cosw sinQ 


—cos? sinw sinQ-+cosw sn 
L +sini sinw 


—cos? cosw sinQ—sinw cosQ 
+cosi cosw cosQ—sinw sinQ -, 
q +sini cosw 


(+sinz sinQ 
—sini cosQ -. 
+cosz 


| 


| we obtain for (1) 


20°F 
= {Q-u.a(1—e&)? cosE—P-u'a sine}. (2) 
t (ua)? 


widering the slowly varying quantities to be constant 
© g one revolution and using 


d/dE= (r/an) (d/dt), 


| 
| 
z 


) tain 
| 20 Be 
4 — {o-wa—e) f coskdE 
n(ua)? Ey. 
E2 
— Pew f sinbak |, (3) 
a eel 
Fr E2 
\-—| Q-uw(1—e)? sinE+P-w? cost | : 
hac p 


ve E, and E» are the eccentric anomaly at exit and 
ince of shadow, respectively, and are obtained from 
- quation 


eae rp"), (4) 
1 terms of £, 


/ (cosE—e)+Q-u0(1—e)? sinE 


=—[(1—e cosE)?—(p/a)?}*, (4’) 


ve p is the radius of the earth. The solution of this 
ition is readily obtained by the iterative technique 
‘the entrance and exit roots are easily identifiable. 
‘e may write the expression for Aa as follows: 


|= (de /y) singe 2) 
XE(Q-w)2(1—2) + (Pew) 

| Xsin[3 (E:+- Es) +], 

arc tan[Q-u?(1—e)!/—P-w°], 


(3’) 


nination of this expression shows (1) if there is no 
owing, 
sin} (H»— £,)=0, 
lting in Aa=0, and (2) if e=0, ¢ becomes 
—32(E1+ £2) 


a consideration of Eq. (4), and again 


Aa=0. 
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COMPLETE SYSTEM OF EQUATIONS FOR THE 
OSCULATING ELEMENTS 


The equations for the elements may be obtained 
using (Strémgren 1929) 


da/dt=2av-F/p, 
de/dt=[(a/u)(1—e) ]?(T cosE+N), 


dt rw 
a ual) } cos(w+f), 
_ dQ rw : 
Ee re sin(w+f), 
day dQ 


1 
+cosi—=—{ Ta? sinE—[a(1—e?) }#L}, 
dt di pe 


where T=F-RXPr°, VN=F-Q, W=F-R, f=true anom- 
aly, L=F-P, and r°=unit vector in direction of r. The 
complete procedure will be carried out for only one 
equation. The others are obtained by exactly the same 
procedure. 


de/dt=(_(a/p)(1—e?) }}(T cosE+N), 
where 
T=— Fu: RX r= — Fu: RX r/r, 
N=—Fu°-Q—FQ-u"(a/r)(1—e cosE£). 
= — (Fa/r)_Q-u°(cosE—e)— P-u°(1—e?)? sinE], 
r= Pa(cosE—e)+Qa(1—e?)? sinE. 
We obtain 
de/dt=—[(a/u)(1—e?) }}F{Q-u?(a/r) 
< ($—2e cosE+# cos2E) 
— P-u°(a/r)(1—e?)? sinE cosE}, 


elttr)—e) _ _F [[e-0] 
Le 


Ti of t 


X{Q-u"(a/r) ($—2e cosE+} cos2E£) 
— P-u°(a/r) (1—e*) sinE cosE}dd, 
where 7 is the period of the satellite. 
Since F is small we may consider slowly varying 
quantities constant during the integration from / to 
(+7. Changing the independent variable for the inte- 


gration and replacing the secant by the tangent on the 
left side, we obtain 


de 1fa i 


$—2e cosE+% cos2E)dE 


E2 
—Pew(—e) f sinE cosh. 
Ey 
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The limits are EZ, to E» since F is zero outside this range. dQ aFR-u? | 3e a | 
: sin oe sini—= — ——E+(1+e*) sinE 
Hence, the resulting equation is di Papua(l melt , 
de a ae . 
=--—|--)] —* sinze | sinw q 
whe ii 4 
3 DiS nal ot e 
x<P| o-w(“E- 2¢ sinE{--E—— sink) lela eT cos2E | cos 
2 4 4 
AG pes Bb, dw dQ ak il 
P-ui(i—e) COsle ~ +008 = ——— | O-ui(¢ cosH—— cost) 
Ee ae dt Qarpte 4 
The complete set of equations is 3E : 1g 
“ : = P-w'(1=2)!(——c sink-— sind 
da wFP Ba z wee 2 4 
SS -uw9(1—e?)? sinE P-we E : 
dt art eam Nat Pig \ The integration of this system of equations must 
be carried out using a large-scale computer. 
de [a(1—e?) |? E 3E ee: 
Fe =— sy P| Q-u (— 2e me sind ) LoeNOuEee ee 
Pe! ra es The author wishes to express his gratitude to. 
be ‘u’(1—e’)} ened z| Ann Bailie for her work in verifying the developm 
2 EY 
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